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SERIES EDITOR’S PREFACE 


THERE is at present a growing interest in exact results in statistical 
mechanics and I am, therefore, very pleased to incorporate in the 
International Series of Monographs in Natural Philosophy the present 
volume by Dr. Bogolyubov in which he presents a detailed study of 
the use of model Hamiltonians which allow one to obtain solutions 
which are exact in the thermodynamic limit as the volume becomes 
infinite. The author has himself made many contributions to this 
field and his account bears the mark of an authority; it will, no doubt, 
be welcomed by workers in the field of statistical physics. 


Oxford D. TER HAAR 


PREFACE 


IN THIS book methods are proposed for solving certain problems 
in statistical physics which contain four-fermion interaction. 

It has been possible, by means of “approximating (trial) Hamil- 
tonians”, to distinguish a whole class of exactly soluble model systems. 

An essential difference between the two types of problem with 
positive and negative four-fermion interaction is discovered and 
examined. The determination of exact solutions for the free energies, 
single-time and many-time correlation functions, T-products and 
Green’s functions is treated for each type of problem. 

The more general problem for which the Hamiltonian contains 
some terms with positive and others with negative four-fermion 
interaction is also investigated. On the basis of analysing and general- 
izing the results of Chapters 1 to 3, it becomes possible to formulate 
and develop a new principle, the minimax principle, for problems in 
statistical physics (Chapter 4). 

I take this opportunity to express my gratitude to Dr. P. J. Shepherd 
for his comments and thorough checking of the manuscript. I also 
express my sincere thanks to Academician N. N. Bogolyubov for 
valuable remarks and to Dr. D. ter Haar for useful discussions. 

The material of this book was the subject of theoretical seminars 
and lectures at Moscow State University. I consider it my pleasant 
duty to thank Professors A. N. Tavkhelidze (Corresponding Mem- 
ber of the USSR Academy of Sciences), A. A. Sokolov, I. M. Ternov, 
I. L. Bazarov and N. A. Chernikov, Dr. V. D. Kukin, and also the 
participants in the theoretical seminars at the V. A. Steklov Mathe- 
matics Institute and the Theoretical Physics Laboratory of Joint 


ix 


Institute for Nuclear Research at Dubna for their interest and en- 
couragement, and for discussions. Some of the results of this work 
where discussed in a theoretical physics seminar at the E. T. H., 
Zurich. The author is indebted to Prof. Dr. R. Jost, Prof. Dr. K. 
Hepp and Prof. A. M. Jaffe for discussions. 


Moscow N. N. BOGOLYUBOV, Jr. 


INTRODUCTION 


§ I. GENERAL REMARKS 


In many-body theory most problems of physical interest are rather 
complicated and usually insoluble. Model systems permitting a mathe- 
matical treatment of these problems are therefore acquiring consider- 
able interest. 

Unfortunately, however, in concrete problems in many-body theory 
there is usually no adequate correspondence between a real system and 
its mathematical model; one must be content with a model whose 
properties differ substantially from those of the real system and in 
solving problems one must use approximate methods lacking the 
necessary mathematical rigour. 

Of considerable interest in this connection is the study of those 
few models which have some resemblance to real systems yet admit 
exact solution; fundamental properties of many-body systems can be 
established in this way. Systems of non-interacting particles can be 
taken as examples of systems which can be solved exactly. Although, 
of course, this model seems rather trivial, it is used as a starting point 
in most problems in many-body theory. In the theory of metals one 
can often leave the mutual interaction of the valence electrons out of 
consideration. In the shell model of nuclei in its simplest form one 
can explain many general properties of nuclear spectra without 
introducing interaction between the particles into the treatment. 

One of the most important problems in statistical physics is the 
study of exactly soluble cases. The fact is that this study makes an 
essential contribution to our understanding of the.extremely complex 
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problems of statistical physics and, in particular, serves as a basis 
for the approximate methods used in this field. 

Up to the present time, the class of exactly soluble dynamical model 
systems has consisted mainly of one- and two-dimensional systems. 

In this study we shall concentrate on the treatment of certain model 
systems of a general type which can be solved exactly, e.g. models 
with four-fermion pair interaction which have as their origin the 
BCS model problems and are applicable in the theory of superconduc- 
tivity; the determination of asymptotically exact solutions for these 
models has been investigated by Bogolyubov, Zubarev and Tserkov- 
nikov (1957, 1960; see also Bogolyubov, Tolmachev and Shirkov, 
1958). 

In these papers an approximation procedure was formulated in 
which the ideas of a method based on the introduction of “approxi- 
mating (trial) Hamiltonians” were propounded and reasons were 
given for believing that the solution obtained was asymptotically 
exact on passing to the usual statistical mechanical limit V + 0. 

Bogolyubov, Zubarev and Tserkovnikov (1957, 1960) considered 
the case of a model with the Hamiltonian 


H = Hot Ain. 
Hy = > (E(p)— 1) ap.aps, (1) 


1 ; 
Hin, = — »y Vv FP, p')a*p, 112 p, +112 Fp’, 412 9-p, -112» 
P 


where a, 44/2. and a+ 44/2 are Fermi operators and V is the volume 
of the system. The kernel (p, p’) is assumed to be a real bounded 


function which effectively vanishes outside a certain region of mo- 
menta. The summation in Hj, is performed over the momenta p and 


t See Lee and Yang, 1952; Yang and Yang, 1966; Lieb and Liniger, 1963; 
McGuire, 1964; Lieb, 1963; Yang, 1967, 1968, 1970; and Sutherland, 1968. 

t Although the term “approximating Hamiltonian” is the literal translation 
of the Russian and is the form found in most of the English translations of the 
author’s and other Russian articles, the term “trial Hamiltonian”, which accords 
better with British and American usage, will be used in this book. (Translator.) 
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p’ belonging to the energy layer: 
Er—w < E(p) < Ert+o. 


It was shown that for such a Hamiltonian it is possible to construct the 
free energy asymptotically exactly (as V +o). The idea of this 
method consisted in the introduction of a so-called “trial Hamiltonian” 
H,(C), which is a quadratic form of Fermi operators and depends on 
certain arbitrary constants C. 

Such a Hamiltonian is easily diagonalized and it is then not difficult 
to calculate the free energy corresponding to it. In this work reasons 
were given for assuming that the free energy Fo(C) becomes equal to 
F as V +0, This result was obtained by means of perturbation 
theory. The derivation was based on the fact that in the perturbation 
theory series each of the terms used to calculate the correction to this 
solution is asymptotically small when V + oo. However, the question 
of the convergence of the perturbation theory series was not investi- 
gated and a more rigorous basis was required for the results of this 
work. In view of this, Bogolyubov, Zubarev and Tserkovnikov (1960) 
treated the problem without applying perturbation theory. 

The system studied was of the BCS type and had the form 


=D Napa—zy LIALS apatyayaytrd, — @) 


where 
a= 4 y w(f) (a_yay+ aja*,); v=0. 


Sf =(, 0), —f = (—p, —09), o is the spin index, taking the values + and 
— 4, p is the momentum, which takes the usual series of quasi-discrete 
values p, = 27n,/L, where, for a fixed value of L(L? = V), n, runs 
over all integers, T(f) = (p?/2m)—, wis the chemical potential, and 
a, and a+ are operators satisfying the commutation relations of 
Fermi-Dirac statistics. Jf, f’) and w(f) are real functions with the 


properties 
FSS) = FLA) =-H SS), 
w(—f) =—w(f). 


4 A METHOD FOR STUDYING MODEL HAMILTONIANS 


For example, 


IS, f’) = 2p, p’) [5(o-0')— 6(o+0')], 
Fp, P’) = FP’, P) = H-p, P’), 


where (o—o’) is the Kronecker symbol. 

In this work a chain of coupled equations for the Green’s functions 
was studied and it was shown that the Green’s function for the 
“integrable problem” with Hamiltonian Hp satisfies the whole chain 
of equations for the exact Hamiltonian (2) with an error of order 
(1/V). From a purely mathematical point of view, however, it is 
obvious that reasonings of this type are not entirely final. Never- 
theless, these papers are a substantial contribution to the study and 
determination of the asymptotically exact solution. 

We remark that arigorous proof of the asymptotic exactness of the 
results of Bogolyubov, Zubarev and Tserkovnikov (1957, 1960) pre- 
sents substantial mathematical difficulties. 

Such a problem was studied by Bogolyubov (1960 a, b) in a purely 
mathematical formulation in the case of zero temperature. He treated 
the model system characterized by the Hamiltonian (2) with the 
assumption that the kernel Z(f, f’) has the simple factorizable form: 


HAS, f') = AP) AP’). (3) 


For the supplementary conditions it was assumed that the functions 
A(f) and T(f) satisfy the following general conditions: 


t From a physical point of view, it should be noted that the simple factorization 
of the kernel (3) describes so-called s-state interactions; when one needs to take 
account of higher states (e.g. p- and d-waves), the kernel must be written in the form 


TLL) = ecb 4eF, fa > 0. 


Some model systems of this type have been treated in the zero-temperature case, 
i.e. for the ground state of the system; this treatment, however, succeeded only 
in determining the ground-state energy asymptotically exactly as V + oo (Tareeva, 
1963, 1964). 
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(1) The functions A(f) and T(f) are real and have the symmetry 
properties : 


M—f)=-Af), T(-f) = Th), 


@) yee} fe ine 
; ) 
(3) VY [Ar| <const for V+, 


(4) tim >a Ss 
(V+) VF VIAf)x+Tf) 
for sufficiently small positive x. 

This model problem (2) with the kernel (3) and the conditions (4) 
was investigated fully for zero temperature. 

It was found that this problem can be solved asymptotically exactly 
as the volume of the system tends to infinity, i.e. the ground-state 
energy and the Green’s functions and correlation functions charac- 
terizing the dynamical behaviour of the system can be determined. 

We remark that in statistical physics we are always interested in the 
limiting values (as V > -) of the dynamical quantities and functions 
being studied. In this work, Bogolyubov (1960 a, b) also showed, for 
application to the above case, how quasi-averages should be intro- 
duced; namely, source-terms were added to the Hamiltonian (2), pro- 
portional to pair-creation operators and pair-annihilation operators. 
After passing to the limit V -- ©, the intensity of the sources was set 
to zero and a definition of quasi-averages was thereby constructed. 

Of considerable interest is the analogous problem for arbitrary 
temperatures, i.c. 9 ~ 0 or 6 = 0. In this case, however, the most 
simple generalization of this problem is found to be impossible. 
The technique applied in the investigation of this problem was inher- 
ently limited to treating the ground state. Papers (Bogolyubov, Jr., 
1965, 1966a) have been devoted to this complicated problem for model 
Hamiltonians of the form (2) (with attraction) in the case of non-zero 
temperatures. 


> 1, 
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In these papers a detailed study was made of a generalized model 
system (with attraction close to the Fermi surface), characterized by 
the Hamiltonian 

H=T-W Y gJJi. (5) 
l<c<s 
This system goes over to the usual BCS model system studied in the 
theory of superconductivity, if we take as the operators T and J, the 
quadratic Fermi-operator forms 


T=ST fata Jum ap A fiapaty 6) 


Then for the Hamiltonian (5) the kernel Z(, f’) will be 
KAS) = oe Bl ) Af’). (7) 


In the above papers the case with all g, = 1 was studied. It is clear 
that with the appropriate choice of functions A,(f) we could obtain 
the systems studied by Tareeva (1963, 1964) at zero temperatures, 
and also more general systems' in which “pairs” interact not only in 
s-states, but also in p-, d-states, etc. 

For the model problem (5) under consideration, we have constructed 
a trial Hamiltonian 


H° = T-2V BAC SE+CI)+VW YY glCe?. (8) 
1 s 1 


Sax sass 


The complex constants C,(1 =< « «<s) occurring here are determined 
from the condition for the absolute minimum of the function 


Sul(Ci, sees Cs) = 5 In Tr e~ #0 (9) 


t See Pitaevskil, 1959; Privorotskil, 1963; Anderson and Morel, 1960, 1961; 
Brueckner and Gammel, 1958; Anderson, 1962; Balian, Nosanow and Wertheimer, 
1962; Brueckner, Soda, Anderson and Morel, 1960; Cooper, Mills and Sessler, 
1959; Klein, 1962; Markowitz and Kadanoff, 1963; Thouless, 1960; Gor’kov and 
Galitskil, 1961; Yang, 1962; and Bogolyubov, Jr., 1966b. 
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in the region of complex quantities (C,, ..., C,). On the basis of the 
results of papers by the author (Bogolyubov, Jr., 1965, 1966a), we 
have developed a technique by means of which, with sufficiently 
general conditions imposed on the operators T and J,, it is possible 
to find an upper bound for 4 = fyo—fy, the difference in the free 
energies per unit volume. It follows from this bound that this difference 
tends to zero as V +o. The results obtained for finite s have been 
generalized further to the case when s = . The results found were 
valid for both cases, 6 > 0 and 6 = 0. 

As we have shown, for such systems it is possible to construct or 
find an expression for the free energy, which, since the technique has 
been developed for arbitrary temperatures, can be used in the treat- 
ment of phase transitions. 

Furthermore, using the example of the model system (5) with 
s = 1, we (Bogolyubov, Jr., 1966a) have indicated a way of obtaining 
asymptotically exact estimates for the simplest, binary, correlation 
functions. 

In a series of papers (Bogolyubov, Jr., 1966c, 1967, 1968 a, b, c, d, 
e, f, g, 1969, 1970 a, b, c, d, e, f), which form the basis of this book, 
we have developed a new method which permits us to study the systems 
(5) with kernels of the type (7) at any temperature, in cases when the 
parameters g, are positive, negative or of different signs. We have 
succeeded in deriving asymptotically exact expressions not only for 
the free energies but also for the many-time correlation functions and 
Green’s functions which completely characterize the dynamica! be- 
haviour of the system. 

In the process we have made clear the necessity of introducing 
a new definition of quasi-averages, since the definition accepted until 
now has proved to be inadequate in the cases we have considered. 

In order to unify all the model problems we have studied and to 
simplify appreciably our account of this material, the plan of the 
book has been arranged as follows. 

All the model problems being treated are reduced to a general 
form, characterized by the Hamiltonian 

Pr=I,4+M (10) 
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where 
r,= ra y {A*(f)a_rart+ A(f)atat,}+ const, (11) 


H, =—V Y.GAJ.—C,) Jt -C), (12) 
T=PT(Najay, J. = wv) Aa fapacy, (13) 


and the function A(/) satisfies the conditions I referred to below. 
The notation is the same as in the treatment of the model system (2); 
a, and af are Fermi operators, V is the volume of the system, and 
the summation over f is performed over the quasi-discrete set oy. 
In the sums over «, the index « takes integer values. 

This model system is treated under certain supplementary condi- 
tions (conditions I) (Chapter 1, § 1), of which the main one is 


© |G.|(0.-C)UE-C2)) «ey, (14) 
a r 


where €y + 0 as V +o. The symbol (...), implies statistical aver- 
aging over the Hamiltonian I’. These conditions are chosen in such a 
way that they are satisfied by all the problems considered later, in 
Chapters 3 and 4. 

We note that the reduction of Hamiltonians to the form (10}(13) 
and the proof that condition (14) is fulfilled require considerable 
mathematical resources. 

For the technique of Chapters 1 and 2 to be applicable, it is necessary 
only that the conditions I be satisfied; then the results of Chapters | 
and 2 and the Theorems I-IV which are proved there are automatically 
extended to the cases of model systems with negative four-fermion 
interaction (cf. Chapter 3), model systems with positive four-fermion 
interaction, and also to the mixed case when the Hamiltonian contains 
both “attraction” and “repulsion” terms (cf. Chapter 4). 

As has been pointed out, all the dynamical problems studied in this 
work can be reduced to the model Hamiltonian (10)-(13). This 
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Hamiltonian consists of two parts, one (11) quadratic in the Fermi 
operators, and the other (12) quartic in the Fermi operators. 

We note that the quadratic form (11), which will be referred to 
below as the trial Hamiltonian, can be easily diagonalized by means 
of the well-known w—v transformation techniques, and for this system 
it turns out to be not difficult to calculate explicitly the corresponding 
expressions for the free energy, correlation functions, etc. 

The quartic form (12), which describes the interaction, is in some 
sense, because of condition (14), a small correction. 

It should be emphasized that an extremely complicated mathemati- 
cal technique was required to study the “effect” of this correction 
on the free energy, correlation functions and Green’s functions, and 
to show, further, that this effect disappears in the passage to the limit 
V+, 

The main feature of Chapters 3 and 4 is the reduction to the “uni- 
versal form” (10)}{13) of the model Hamiltonians dealt with in these 
chapters and an actual proof of condition (14) for the specific systems 
considered there. 

In Chapter 1 we examine earlier methods (Bogolyubov, Jr., 1967) 
and the difficulties associated with the introduction of quasi-averages 
for certain particular model systems with four-fermion negative inter- 
action. In this connection we formulate a new principle for treating 
model systems of the form (10)(13) under sufficiently general condi- 
tions (cf. Chapter 1, § 1). The model systems (10)-(13) represent a 
fairly broad class of model systems, including the following important 
particular cases: (a) systems with four-fermion negative interaction 
(e.g. the BCS model systems studied in the theory of superconduc- 
tivity), (b) systems with positive four-fermion interaction, and also 
(c) the mixed case for which the Hamiltonian contains both terms 
with “attraction” and terms with “repulsive” interaction. 

For this general class of model systems we demonstrate how single- 
time and many-time correlation averages can be determined asympto- 
tically exactly (as V + oo). In passing we point out the possibility 
of obtaining upper bounds which are uniform in 9. 

In Chapter 2, using the results and upper bounds obtained in 
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Chapter 1, we prove generalized asymptotic relations for many-time 
correlation averages (composed of products of field functions), T-pro- 
ducts and Green’s functions. 

Noting that the difference n,—n_,in the numbers of particles with 
opposite momenta is an integral of motion for the systems (10)-(13), 
we formulate here special selection rules for averages composed of 
products of Fermi creation and annihilation operators. 

On the basis of this treatment we formulate Theorem III and also 
the Theorem IV for the existence of generalized limits for averages 
taken over the trial Hamiltonian. 

In Chapter 3 we examine a concrete application of the results of 
the first two chapters to model systems with four-fermion negative 
interaction. We draw attention here to certain difficulties associated 
with the definition of quasi-averages for model systems. 

We develop a new method of introducing the auxiliary system and 
defining quasi-averages for a broad class of model systems containing 
negative interaction. 

In §§ 1-3, we consider the question of the asymptotic expression 
for the free energy as VV > o. 

A model system, characterized by the Hamiltonian (5), is studied. 
we treat the case when the operators T and J, have the form (6), and 
all the g, = 1. 

In§ 1 of Chapter 3 we study certain properties of the free energies 
derived on the basis of the model Hamiltonian (5) and the trial 
Hamiltonian (8) (with all g, = 1). The results concerning the closeness 
of the free energies constructed on the basis of the Hamiltonians (5) 
and (8) are briefly formulated. The bounds obtained are found to be 
uniform in the temperature 6. We note, furthermore, that the existence 
of an asymptotic expression for the free energy constructed on the 
basis of the Hamiltonian (5) does not follow from the inequalities 
for the difference in the free energies constructed on the basis of the 
model Hamiltonian (5) and the trial Hamiltonian (8). 

In this connection, in Chapter 3 we formulate those conditions for 
which we can give a correct proof of the existence of a limit as the 
volume tends to infinity for the free energy constructed on the basis 
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of the model system (5), i.e. 
lim f(H). 
(V — oe) 


In§ 7 of Chapter 3 we consider the question of the definition of the 
quasi-averages. Let @ be some operator of the form for which the 
limit theorems are formulated in Chapters 1 and 2, e.g. a product 
of Fermi amplitudes or field functions. 

It has been proposed that the quasi-average < (> of such an op- 
erator for the above Hamiltonian be defined as the limit 


(» +0) \(V — 0) 
of the usual averages (7), taken over the Hamiltonian I” obtained 
from (5) by the addition of “source-terms” 


Pr =H-VY (,Ji+vtJ,). (16) 


In§ 7 of Chapter 3 we draw attention to the difficulties associated 
with the definition (15). In this definition there is, for example, no 
indication of the region in which the parameters » must lie or of how 
they must tend to zero in order to ensure convergence in the definition 
(15). It is shown in § 7 of Chapter 3 that, even in the simplest cases 
of binary averages for a kernel of the form (3), if » tends to zero 
arbitrarily, the limit in the expression (15) may not exist, since this 
expression is a finite statistical average multiplied by the factor »/|»|; 
the presence of this factor means that as » tends to zero the expression 
(15) does not tend to any limit. The limit does exist, for example, when 
» tends to zero in such a way that the ratio »/|»| remains constant. 
In the general case of a non-factorizing kernel (i.e. when the number 
of terms s > 1), with the passage to the limit » - 0 the situation is 
found to be even more complicated. Apart from gauge invariance 
(due to the gauge group), other groups of transformations, e.g. the 
rotation group, may also be manifested. 

To avoid this type of difficulty, in § 8 of Chapter 3 we propose that 
» be taken proportional to C with a constant coefficient of proportio- 
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nality. Accordingly, it is shown in § 8 that the auxiliary Hamiltonian 
with sources takes the form 


P= H+ Sigut(J.—C.) (JZ —Cd, (17) 


in which the t, are positive parameters (« = 1, 2, ..., s), and H is 
given by the expression (5). We emphasize that C, denotes the point 
at which the absolute minimum is attained of the function (9) taken 
after the passage to the limit V > , so that C, does not depend on 
V. We see that with the above choice (17) of the auxiliary Hamiltonian, 
no difficulties arise in defining the quasi-averages. These results are 
summarized by the following theorem: 


THEOREM 3.1V. If the conditions of Theorem 3.1\1 are fulfilled and 
the Hamiltonian is represented by expression (17) in which 


0<17,< 1 (a = 1,2, ..., 5), (18) 


then the following inequalities are valid: 


0<f(T)—f(H) < 5y+ by + 0 (19) 
as V+, 
Ledde-C)UE-C)r= Ht 20 0) 
; as Vo, 
where Tt, is the smallest of the quantities t,, ..., T,. 
On the other hand, 
Pr = H(C)-2V » 81 —7.) (Ja—C,) (JE —Ca), (21) 


where H(C) is the trial Hamiltonian, defined by expression (8). By 
virtue of the above result for the smallness of the “correction term” 
(20) quartic in the Fermi operators, all the conditions I of Chapter 
1 will be fulfilled and consequently all the results of Chapters 1 and 
2 are automatically extended to the above case. 
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To summarize, we shall write down theorems on the existence of 
limits for the quasi-averages 


< U ~H= lim (Wr = lim (Wc); (22) 
(V—> ce) (V — 09) 


in which t can take any values from the range (18). 

We draw attention to a surprising feature of the asymptotic relations 
(22). In spite of the fact that the Hamiltonian (21) depends on the 
parameters t,, the expression in the left-hand side of the asymptotic 
relations (22), which is constructed on the basis of this Hamiltonian, 
is found to be independent of the parameters t, lying in the range 
(18). 

Using the new definition (22) of the quasi-averages and taking into 
account that H(C) is independent of the parameters 1, lying in the 
range (18), we can write the old definition in the form 


lim ( lim (Wr) = lim (Macey. (23) 
(«+ 0) \(V +) (V > 0) 

The expression in the round brackets in the left-hand side of (23) 
always has as its limit the right-hand side of (23), provided that the 
parameters t, lie in the range (18). Consequently, when all the t, tend 
to zero while remaining positive, the limit as 7 — 0 is trivial. 

The main argument of § 7 of Chapter 3 is the proof of the inequality 
(20), which is based on the inequality (19). It is not difficult to see that 
the asymptotic relation for the free energies (as V + -) also follows 
from the inequality (20). 

We remark that for our method it is not necessary to confine our- 
selves to the case g, > 0. In Chapter 4 we examine the situation 
when the constants g, have different signs, and in this case it is also 
possible to reduce the Hamiltonian to the form treated in Chapters | 
and 2. 

In Chapter 4 we consider model systems with positive and negative 
interaction components. In § 1 we investigate the properties of the 
free energy for a system with positive four-fermion interaction. We 
prove the existence of an absolute maximum for the free energy 
constructed on the basis of the corresponding trial Hamiltonian. 
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It is shown that the solutions of the equations (formulae (4.9)) for 
the point C at which the absolute maximum of the free-energy function 
occurs are unique. 

In § 2 we examine the question of whether the analogue of Theorem 
3.1 is valid for the case of a model system with positive interaction 
components when the operators T and J have arbitrary unspecified 
form. 

It is shown that, in the general case when the operator forms of T 
and J (for interaction of the repulsive type) are unspecified but satisfy 
the same general conditions as for Theorem 3.I, the free energy con- 
structed on the basis of the model Hamiltonian does not tend, as 
V — co, to the free energy constructed on the basis of the correspond- 
ing trial Hamiltonian. 

In § 3 of Chapter 4 we investigate the problem of calculating the 
free energy and correlation averages for a model system with positive 
four-fermion interaction, corresponding to repulsion of fermions. 
The essential difference between this and the case of systems with 
negative four-fermion interaction (cf. Chapter 3) is that the free energy 
constructed on the basis of a trial Hamiltonian regarded as a function 
of the variables ...C,... has an absolute maximum for a certain 
fixed set of quantities ...C, ..., and the maximum point is unique. 
As a result we obtain upper bounds for the free energy and quasi- 
averages and from these there follows an unexpected result: the 
contribution to the Hamiltonian of the system from a four-fermion 
positive pair interaction is found to be negligibly small for large 
systems (V - ) (cf. § 3). 

On the basis of this treatment, in § 4 we also investigate a model 
system of a more general form (4.30). This study is essentially auxiliary 
for § 6, in which we formulate the minimax principle. As a result, we 
obtain bounds for the free energies and quasi-averages. 

In§ 5 we prove that the solutions of the equations (4.47) are unique 
and also show that these solutions C(U) are continuous functions and 
possess continuous partial derivatives with respect to the variables 
...U,... and ... Us .... In this proof we make use of a theorem 
concerning implicit functions (since all the conditions for its applica- 
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bility are fulfilled). The results obtained are summarized by Theorem 
41. 

In § 6 we investigate Hamiltonians with coupling constants of 
different signs and formulate the minimax principle. We construct 
an auxiliary system and set up the asymptotic relations for the free 
energies and quasi-averages. An expression is found for the asymptoti- 
cally exact calculation of the free energy using the minimax principle 
(formula (4.73)). 

We emphasize that in deriving a proof for the determination of an 
asymptotically exact solution for a model system (5) in which the 
g, take different signs, we start from a specially constructed trial 
Hamiltonian H(C, S) which is a quadratic form in the Fermi operators 
and depends on the groups of complex variables C,, ..., C,, and 
S,, ..., S,. These variables, which appear initially as independent of 
each other, are then related by means of the minimax technique, 
namely, from the condition for the absolute maximum of the free 
energy with respect to the variables C,, ..., C,, for fixed S,, ..., S, 
we obtain a system of equations from which a unique solution C = C(S) 
can be determined. The solution found is then put into the trial Ha- 
miltonian. Using this trial Hamiltonian, we construct the free energy 
and in the final upper-bound inequalities take the absolute minimum 


min f{H(C(S), S)} 
(S) 


of the free energy with respect to the variables S,, ..., S,. 

It is found that taking the maximum and minimum in the reverse 
order leads to an incorrect result. 

We also consider the question of the definition of the quasi-averages 
and the construction of a Hamiltonian with sources. As a result we 
obtain upper bounds, uniform with respect to the temperature, for 
the corresponding quasi-averages. It is shown that this problem re- 
duces to the Hamiltonian (10)-(13) and, since all the supplementary 
conditions I are fulfilled, we can apply the limit theorems of Chapters 
1 and 2 to this case and obtain the same results for the quasi-averages 
as in Chapter 3. 
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In this book we have solved the complicated problem of constructing 
upper bounds for many-time averages for the case of a fairly broad 
class of model systems with four-fermion interaction. The methods 
proposed in this book for solving this problem will undoubtedly 
find application not only for the model systems, associated with the 
theory of superconductivity, considered here. The theoretical methods 
developed in Chapters 1 and 2 are applicable to a much broader class 
of model systems from statistical physics and the theory of elementary 
particles. 


§ 2. REMARKS ON QUASI-AVERAGES't 


The expression “quasi-average” will be used frequently in this book; 
we therefore include this introductory section to clarify, without 
the required mathematical rigour, the physical meaning of the concept 
of a “quasi-average” and of the connection between quasi-averages 
and the question of the degeneracy of a state of statistical equilibrium. 
This section is of an introductory nature and has no direct connection 
with the main text of Chapters 1-4. 

Here, using a number of simple model examples borrowed from 
statistical physics, we shall elucidate the physical meaning of the 
quasi-averages introduced in the papers of Bogolyubov (1960 a, b). 

First we shall examine the idea of the degeneracy of a state of 
statistical equilibrium. We remark that the concept of degeneracy 
is well known in quantum mechanics (Heisenberg, 1966). 

In treating problems on the determination of eigenfunctions in 
quantum mechanics, we know that the usual form of perturbation 
theory, developed for problems without degeneracy, cannot be applied 
to problems with degeneracy; for such problems it must be modified. 

Because of the existence of additive conservation laws, there are 
always cases of degeneracy in problems in statistical physics. 

It might seem at first sight, however, that in these problems the 
degeneracy has no effect and need not be taken into account in 


t See Bogolyubov, 1960b, 1963. (See also Selected Papers of N. N. Bogolyubov, 
Vol. 3, p. 174, Naukovaya Dumka, Kiev, 1971 (in Russian).) 


INTRODUCTION 17 


practice. In fact, in problems in quantum mechanics, a linear manifold 
of eigenfunctions can correspond to one energy eigenvalue; in this 
case, the eigenfunctions contain undetermined constants. 

In statistical physics, however, the average value of any dynamical 
variable @ is always uniquely defined by 


Tr U e~ He 
(U) = a =H 


Consequently, the Green’s functions and T-products constructed 
from the usual averages must also be uniquely defined. 

It may appear from this that in studying a state of statistical 
equilibrium, e.g. by means of a diagram technique, we need not take 
the presence of degeneracy into account. In fact, the situation is 
more complicated. 

To form an intuitive idea of the nature of the difficulties arising 
here, we shall consider the case of an ideal isotropic ferromagnet. 

We shall take the Heisenberg model system 


H=-% 4, I fi-F-(Sp,° Sp); (24) 


Is 


where the f are spatial points corresponding to the lattice sites (lying 
in the volume V), S, is the spin operator of an electron at site f, and 
I(f1—f2) is the exchange integral. For definiteness we put I(fi—f2) = 0 
for all f; and fo. 


We note that, for the system (24), each component of the total spin 
SS ey (a = x, y, Z) (25) 
is an integral of motion. 
We have also 


S89 —S*S* = iS!, S¥S*¥—S#S” = iS*, 


S#S*—S*S* = iS’. (6) 


Taking these commutation relations into account, we write 


iTr (S? e~#/*) = Tr {(S*S’ — S¥S*) e~ #9}, 
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But since S* commutes with H, we obtain 
Tr (S”S* e~ 4/8) = Tr (SY e~ 4S) = Tr (S*S” e~#/9) 


and, therefore, 
Tr (S? e—#*) = 0. 


Repeating these arguments for the components S” and S*, we find 
Tr (S* e—#*) = 0, (a = x, y, 2). (27) 
We introduce the magnetization vector referred to unit volume: 


eo 1 oe AL 
M a hep Lr eS 

We have 
Tr (Mle e- #/*) = 0, 


and, consequently, 


 9—H/0 
Ue) = tin =0, («=%y%2), (28) 
i.e. the mean magnetic moment of the system is zero, which corresponds 
to the isotropy of this dynamical system with respect to the spin- 
rotation group. We emphasize that the expression (28) is valid for all 
temperatures 6, and in particular for temperatures below the Curie 
point. It is known, however, that at temperatures below the Curie 
point the magnetic moment of the system is non-zero, although its 
direction can be chosen arbitrarily. In this sense, a state of statistical 
equilibrium in this system is degenerate. 

We now switch on an external magnetic field ve (v > 0, e? = 1), thus 
replacing the Hamiltonian (24) by 


H,. = H+v(e-MV. (29) 


Bearing in mind a characteristic property of isotropic ferromagnets 
at temperatures below the Curie point, we see that 


(M*\ = e*M,, (a = x,y,z) (30) 
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where M, will tend to a finite (non-zero) limit when the intensity » 
of the external magnetic field tends to zero. 

In this case we thus have an “instability” in the usual averages; on 
addition of the term »(e-JOV with infinitesimally small »*t to the 
Hamiltonian (24), the average (") experiences a finite non-zero 
increment: 

em; m = lim M,. 
v0 

For the system with the Hamiltonian (24) we introduce the concept 
of a “quasi-average”. 

We take the dynamical variable of to be a linear combination of 
the products 

A = S(t)... SH(t). (31) 


We shall define the quasi-average < of > of this quantity as 
<A> = lim (Ave, 


(v +0) 
where (c#),, is the usual average of o# for the Hamiltonian H,,. 
Thus, the presence of degeneracy is directly reflected in the quasi- 
averages by their dependence on the arbitrary unit vector e. 
It is not difficult to see that the quasi-average is related to the 
usual average (c#) by 


(A) = J <A> ae. (32) 


It is clear that quasi-averages are more convenient for the descrip- 
tion of the above case of degeneracy of a state of statistical equilibrium 
and also more “physical” than the usual averages. The latter are the 
quasi-averages averaged over all directions e. 

The usual averages 


(Sil) .. SHC) (33) 
must be invariant with respect to the spin-rotation group. 


+ When speaking of infinitesimally small », we mean that » tends to zero after 
the statistical mechanical passage to the limit V -> oo. 
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The corresponding quasi-averages 
< Sita)... SH(t)> (34) 


will only possess the property of covariance; on rotation of the spins, 
the vector e must be subjected to the same rotation, so that the expres- 
sion (34) is unchanged. Thus, the selection rules for the usual averages, 
resulting from the invariance of the latter with respect to the spin- 
rotation group, will not apply to the quasi-averages. It is seen that 
the unit vector e, i.e. the direction of the magnetization vector, 
characterizes the degeneracy of this state of statistical equilibrium. 
To remove the degeneracy, the direction of e must be fixed. We shall 
take the z-axis for this direction. In this case, all the quasi-averages 
become well-defined numbers. It is with these averages that the theory 
of ferromagnetism is concerned. In other words, to remove the degene- 
racy of a state of statistical equilibrium with respect to the spin-rota- 
tion group, we must include in the Hamiltonian A an invariant 
extra term »/,V with infinitesimally small y. In the above example 
of an isotropic ferromagnet, the perturbation term is the energy of 
interaction of the system with an external magnetic field. 

We turn now to consider a second example containing degeneracy, 
taken from the theory of the crystalline state. We shall take a dynamical 
system of spinless particles with two-particle interaction, characterized 
by a Hamiltonian in the second-quantized representation 


H FP 4\ gt 
= (Smi-2) 84% 
1 ’ , , \ 

+55 DY 95,95 ,4p,4p( Pi — Pi) 8(P1+P2+Pi—P2), (35) 

2V vy, per Pin D2 
in which d(p) is the Kronecker 6 and v(p) is the Fourier transform 
of the potential energy of interaction G(r) of a pair of particles. We 
shall further assume that this interaction is such that our dynamical 
system is crystalline at fairly low temperatures 0 < 6,,. 

We shall consider the particle-number density e(r), which must 
clearly be a periodic function of r with the period of the crystal lattice. 
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It would seem natural to assume that g(r) is equal to the usual average 
value of the density operator Y*+(r) P(r): 


or) = HPO) = FLT Catan deler. 36) 


This, however, is not so. 
In fact, in this case a conservation law for the total momentum is 


valid: 
K= vk ata,, HK—KH =0. 
k 


Taking into account the translational invariance of averages of the 
form 


(MH) = (P(t, r+ &) Pt r+ €)), (37) 
where 
P(t, r) = a Y a(t) ett, (38) 


it is not difficult to arrive at the selection rules: 


(afaz,g=90 if gO, 
whence 


(Ptr) Pir)) = 5 py (ata,) = se = const. (39) 


Thus, the usual average of the density operator cannot be equal 
to the periodic function g(r). It is clear that this position has come 
about because of the degeneracy of the state of statistical equilibrium 
under consideration. 

A crystalline lattice can, as a whole, be positioned arbitrarily in 
space. In particular, our Hamiltonian possesses translational invari- 
ance and, therefore, the lattice can always be given an arbitrary trans- 
lation. There is nothing to distinguish any particular spatial position 
of the crystal lattice, and when we take the usual average, we average 
at the same time over all possible positions of the lattice. To remove 
the degeneracy, we include in the Hamiltonian (35) a term of the form 


vf Un) Yt) Pid; »>0, »+0 (40) 
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corresponding to an infinitesimal external field vU(r). For U(r) we 
take a periodic function of r with the periodicity of the lattice, so that 
the external field vU removes the degeneracy, fixing the position of the 
crystal in space. 

Since we naturally consider only physical, stable systems, it is clear 
that the inclusion of an infinitesimal external field can only produce 
an infinitesimal change in the physical properties of the system being 
studied. 

By taking the usual average of the density operator + (r) V(r) for the 
Hamiltonian H, with infinitesimal », we in fact obtain an average 
for the system with the original Hamiltonian H, but without the 
additional averaging over the spatial positions of the crystal lattice 
(as a whole), since its position is now fixed. 

In this way we obtain the observed spatial particle-distribution 
density g(r). We shall define the quasi-averages formally by 


<e.. Prt, r;) Pent Wt, rs) > 
= lim (... W(t, nj)... Mts, rs)... ay. 


(» + 0) 


Then, as we have already remarked, 
<P*(r) P(r)> = oft). 


Bearing in mind that 


1 


<¥t(r) P(r)> = y > {z <ahaurg> beter 
q 


we see that the quasi-averages 
(ata,  (k' # k) (41) 


cannot all be equal to zero. 

It is clear that the selection rules arising from the law of conservation 
of the total momentum are not obeyed by the quasi-averages we have 
introduced. 

We note that the quasi-averages depend in general on a number 
of arbitrary parameters, e.g. on the arbitrary vector &. In fact, by 
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replacing U(r) by the permissible function U(r+&), it is not difficult to 
see that the quasi-averages (41) are replaced by 


<apay > eltk-k-8, 


The quasi-averages become uniquely defined when the function U(r) 
is fixed. 

We have considered cases of degeneracy of a state of statistical 
equilibrium associated with a conservation law for the total spin 
vector and the total momentum vector. In both examples, the degener- 
acy could be removed and quasi-averages introduced by including a 
suitable infinitesimal external field. 

We give one further example, when the degeneracy is associated 
with a conservation law for the total number of particles. We take 
the well-known example of the Bose-Einstein condensation of an ideal 
gas. We write the Hamiltonian of the system in the form 


2 
H =—jagayt+ Pee (s~*) ata,, e>0. (42) 

The degeneracy will be lifted if we add the infinitesimal term 
Waytat) /V. 

In the above examples, the degeneracy was associated with the 
existence of additive conservation laws, or, in other words, with the 
presence of invariance with respect to the corresponding groups of 
transformations. We remark that not all the conservation laws apply- 
ing in a given system bring about degeneracy. In the last example 
(the Bose gas), the degeneracy was associated only with the law of 
conservation of the number of particles. In the corresponding quasi- 
averages, only those selection rules which resulted from this law were 
violated. In the second example (the crystal lattice), the degeneracy 
was associated only with the law of conservation of momentum. 
Selection rules due to, for example, the particle-number conservation 
law were not violated. One could extend the list of examples by con- 
sidering cases of degeneracy associated with other groups of trans- 
formations, or with several groups of transformations simultaneously. 
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We shall not dwell on this here. We shall pass to more general con- 
siderations and introduce the appropriate definitions. 

We shall take a certain macroscopic system with Hamiltonian H. 
To H we add some infinitesimal terms corresponding to the external 
fields or sources violating the additive conservation laws and thereby 
obtain another Hamiltonian H, (vy > 0). Then, if all the averages 


(A), ASS 505 wrt, r;) ates Pte, Vs) Pres (43) 


receive only infinitesimal increments, the state of statistical equilibrium 
will be said to be non-degenerate. Conversely, if some of the averages 
(43) receive finite increments in going from H to the infinitesimally 
close Hamiltonian H,, the state of statistical equilibrium will be said 
to be degenerate. We confine ourselves to treating stable systems, since 
only they have physical meaning. 

For cases of degeneracy, it is appropriate to introduce quasi- 
averages in place of the usual averages: 

<A> = lim (A)z,. 
(v +0) 

As we have already seen from the above examples the, selection 
rules resulting from the additive conservation laws are not necessarily 
all obeyed for the quasi-averages. We emphasize that by defining 
quasi-averages we are constructing a two-limit technique: first we put 
the volume equal to infinity, and then we put » equal to zero. 


CHAPTER | 


PROOF OF THE ASYMPTOTIC 
RELATIONS FOR THE MANY-TIME 
CORRELATION FUNCTIONS 


§ 1. GENERAL TREATMENT OF THE PROBLEM. SOME 
PRELIMINARY RESULTS AND FORMULATION OF THE 
PROBLEM 


We begin here with the problem of the asymptotic calculation of 
quasi-averages and prove a number of preliminary theorems. 
For a natural approach to the proper formulations of these we recall 
a series of results which we established earlier for a Hamiltonian of the 
form 
H = T-2VgJJ+ (1.1) 
in which 


2 ee 
T= LTA) afays Tf) = 50 —b 
s l + 
T= ap DMN apat, (1.2) 
a,and a; are Fermi amplitudes, m, u and g are positive constants, and 
f= (p,°) 


where @ is the spin index, which takes the values +4. The function 


Af) = Ap, 2) 


25 
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occurring in (1.2) is real and continuous in the spherical layer 


2 
Pale 
2m 
(where -/ is a certain positive constant) and equal to zero outside it; 


it also possesses the property of antisymmetry 
A-f)=-Af~); —f =(-p, —9). 


We note finally that in the summation “over f” the components p, 
(a = 1, 2, 3) of the vector p take the values 27n,/L, while n, runs 
over all integers (— <0, oo). L3 = V, where V is the volume of the 
system, which below will tend to . 

To introduce the quasi-averages we add to the Hamiltonian H terms 
with “pair sources”, e.g. 


—wW(J+I*) 


where » is a positive constant. 
Thus, the Hamiltonian under consideration will be 


P = T—WegJJt—wW(J+JI*). (1.3) 


The quasi-averages for the Hamiltonian 7 are introduced as the limits, 
as V + oo, of the usual averages for the Hamiltonian J’, with the 
sequential passage to the limit 

<...>y = lim lim (...)r. 


y—->0 V-> oo 


(v>0) 


We have shown (Bogolyubov, Jr., 1967)* that the simplest binary-type 
correlation averages 


(af (Na(t))r, (af (Nat (t))r, (a_(Daf(t))r 
are asymptotically (V — oo) close to the corresponding averages taken 
t We note that this technique has also been found to be useful in the study of 


exactly soluble quasi-spin models (Bogolyubov, Jr., 1968h; Thirring and Wehrl, 
1967; and Thirring, 1968). 
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for the “trial Hamiltonian” :* 
I({C) = T—2Wg(CI* +C*J —C*C)— WJ + J*) 


Eps YN megs Pe * 
_T 2V8 {(c+ 2) 7 +(€ +32) save C. 


(There is, however, another approach, viz. to deal with infinite volume 
from the start. Such a situation was studied in a paper by Petrina 
(1970).) The quantity C occurring in °C) is determined from the 
condition that the free energy has an absolute minimum 


S{EAC)} = min 


in the whole complex C-plane. 
Since [.(C) is a quadratic form in the Fermi amplitudes, this 
Hamiltonian can be diagonalized by means of a u-v transformation: 








ar = Uf ap— Uf az , (1.4) 
where «and a are new Fermi amplitudes, and 
T(f) 
ge. 
u(f) = NE Kf) 
A z 
Vif) pee 1 TA 
v2 | Wp) (C+ 32) Ef)” 
E(f) = | rare? C+ 





In the new Fermi amplitudes the Hamiltonian will take the form 
PAC) = [eC°C~ 5p PED -TI\ VE EDaiay (5) 


Tin this formula we should have written 2 VeC*C- i, where fis the unit operator; 
however, since this will not lead to misunderstandings anywhere, below we shall 
not write out the unit operator explicitly. 
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so that the free energy per unit volume calculated on the basis of this 
Hamiltonian will be 


AIPAC) = 2C"C~ 5p YABLN-T) 
(1.6) 
- : y In (1+ e-/*), 


It is clear from this that the absolute minimum of f{I°,(C)} as a function 
of the complex variable is found at a real value of C, and at a value 
such that 


v 


This minimizing value of C depends in general on V and »: 
C=CY, »). 


In our above-mentioned paper (Bogolyubov, Jr., 1967) it was shown 
that 
CV, ») + C(r) 


as V-+co  (v(fixed) > 0) (1.7) 
and Cl») Gear C0). 


Here C = C(») realizes the absolute minimum of the asymptotic 
expression 


f(F'(C)} = lim FUI(C) 
ae, 
= 2C'C- x5, | aflEf)-T(f)} 
— pap | inter eem 


As throughout the book, the “integral over f” denotes integration 
over p and summation over a: 


fa. =¥ fal...) 
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The value C(0)=O0 is chosen as the number giving the absolute 
minimum of the function 


folHAC) = 28C~r3 5 | MEVI-TD) 


- ene | df \n (1 +e-(/08), 
in which 
E(f) = VT*f)+49°C?2(f). 


As we have noted already, we have proved in our work that the 
difference of the binary averages constructed on the basis of the model 
(L) and trial (1’,) Hamiltonians tends to zero as V + oo for any fixed 
value of »y > 0. 

On the other hand, in view of (1.4) and (1.5), these averages are 
calculated easily for I". 

For example, we have 


e-E(N0 
(a (0) a,(t))r, = uf) BNCrr 1+ eke 


| 


+ v( fle ENE-9 Tpe- ke 5 


(1.8) 
(af (a2 ,(@))r, = WIP) {ene 
: ia 1+e-F(e 


_e-ene-y___ | sd 
1+e-V6 


As can be seen, the right-hand sides here are defined for all p, and 
not only for the quasi-discrete values 








_ (2am, 2nn2g ang 
‘Cra 
which are the only values for which the amplitudes a, and af, and 
thereby the left-hand sides of the expressions (1.8), are defined. 

Moreover, the right-hand sides of (1.8) as functions of f depend on 
V only through the quantity C = C(V, »). 
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The passage to the limit V ~ -o therefore reduces, because of (1.7), 
to replacing C(V, v) by C(v) in these functions. The subsequent passage 
to the limit » > 0 (» > 0) corresponds to replacement of C(») by 
C(0). 

Thus, if we put » = 0, C = C(O) in the functions of f occurring 
in the right-hand sides of (1.8), these functions will represent the 
corresponding quasi-averages for the Hamiltonian H. 

We note that the most complicated of our proofs was that establish- 
ing the relations 


oe Ce ea (1.9) 


as V +0 


for binary expressions of the type indicated above. To prove them, we 
had to show first that 


(J-C(V, »)) (J+ -C(V, »)))r +0 (1.10) 
as V +2 


and then establish the asymptotic relations (1.9a). We must emphasize 
that the problem of investigating situations with more complicated 
averages was not solved in the papers cited (Bogolyubov, Jr., 1966 a, c, 
1967). 

Moreover, the proof of the properties (1.10) and (1.7) was based 
on specific features of the Hamiltonian (1.1) and could not be extended 
to model Hamiltonians of more general form. 

In our work (Bogolyubov, Jr., 1966 a, c, 1967, 1968e), we have 
constructed a new method which enables us to extend the above- 
mentioned results to the case of many-time averages of Fermi ampli- 
tudes or field functions, and, furthermore, for model Hamiltonians 
of more complicated structure. 

If we wished to apply this method to the investigation of the Hamil- 
tonian I’ it would be appropriate to start from the representation: 


P =T,,(C(»))—2gV(J—C(v)) (Jt —C()) 
=T-+ y A(f) {afat + a_pay}+2gVC? 


—2gV(J—C(»)) (Jt¥—C(v)), (1.11) 
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where 


A(f) = 2gA(f) (co+32). (1.12) 


In this case, in view of (1.7) and (1.10) 
((J—C(v)) (Jt -C(»)))r < ev > 0. (1.13) 


However, we shail not study only this Hamiltonian here. 

As will be shown under much wider conditions in the following 
chapters, model Hamiltonians with properly chosen source terms can 
also be reduced to a form similar to (1.11). 

We turn now to consider the situation when 


Fr=T,4+M,, (1.14) 

where 
lr, = y T(f) afay— > y {A*(f) a_ya,+ A(f) atat\+X, 
where 
p 
XKz=const, T(f)= om 7b 
and (1.15) 
A, =-V > G(F.—Czx) (Jf —C%), 


where 
__! tat 
J. = wW > Alf afar. 


Here the summation over f runs over the above-mentioned quasi- 
discrete set, which we shall call the set @,. We shall examine this 
model system under the following conditions, which we shall call 
conditions I. 


1. The functions 2,(f) and A(f) are defined and bounded in the 
whole space @ of the points f = (p, 0). 
2. The series 


Y 1Ge1 LaF)? = Pf), 


converges uniformly in @ and the function P(f) represented by it 
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satisfies the inequalities: 
P(f) < M, = const, 


7 P(f) < M2 = const. 
3. The inequality 
o 1G.|.-C)UL-CD) <er 
* r 
is satisfied, where 
ty>O as Voro, 


4. The function A(/) and the constants C, satisfy the inequalities 


7 LIANE =< Ma 
> |G.||C.|? <= Mc. 


5. The functions A,(f) and A(/) are antisymmetric with respect to 


reflection :* 
A(—f) =—4(f), A(—f) =-A(f). 


Here « takes integer values. (In general, the sums over « imply an 
infinite number of terms; in Chapters 3 and 4, however, we consider 
problems in which « takes a finite number of values.) 

With these conditions we prove a series of theorems on the asymp- 
totic closeness of averages taken over the Hamiltonians I’ and I’, 
respectively. 


+ This last condition is not, essentially, a restrictive one. Any sum of the form 
LF asfaty;  Y Ff) a_sa, 
S IS 


can always be reduced to the form: 
F(f)-F(-f) F(f)—F(-f. 
Le ee FOE a yay 


$3 


in which the coefficient function [F({) — F(—f)]/2 is already antisymmetric with 
respect to the reflection: f + —/, 
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In the following chapters we shall reduce the Hamiltonians studied 
to the form (1.14), (1.15) and, as soon as we have established the validity 
of the conditions I, we can then make use of certain theorems. We turn 
now to the proof of these theorems. 


§ 2. EQUATIONS OF MOTION AND AUXILIARY OPERATOR 
INEQUALITIES 


We turn first of all to the equations of motion for the model Hamil- 
tonian I’. We have 


pH = ar —Ta, = T(f) a,—A(f) at, 


(1.16) 
ce GAdf) at (Jt —C). 
Taking the conjugate and replacing f by —f, we find also 
dat, 
i ara => -T(f) at,—A*(f)a, 
(1.17) 


—2h GAS) (Ja—C,) Gy. 


As can be seen, the first two terms in the right-hand sides of these 
equations originate from the so-called trial Hamiltonian I’,. But this 
Hamiltonian can be diagonalized by means of a canonical u—v trans- 
formation 


ap = u(f) a — Af)aty 





(1.18) 
at,= uf) att o*(f) a, 
in which we put 
T(f) Af) y T(f) 
ies So ie. 
sas Vi By "P= Vayagy Vg)’ a9) 


E(f) = VT*Xf)+/ AYE. 
As a result, this Hamiltonian is reduced to the form 


= y E(f) afay+const. 
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Expressing the “new” Fermi amplitudes «, in terms of the “old” 
a,, we obtain, for example, 


ay = Uf) a+ (f) aty 


Hence, taking eqns. (1.16) and (1.17) into account, we shall have 


iG = E(f) at Ry, (1.20) 


where 


Ry = —u(f) RY — v(f) RP, 
RY = ¥ GAAS) at (Jt -C%), (1.21) 


RY = pS G.AU(f) (Ja—Cy) ay. 


We must also keep in mind the conjugate equation 


1 CF = Kp yap + Rp, (1.22) 


in which 
Rf =—u f) REO— of) RF, 
RPO = VGA) eC) -p, Ge 


RAM = Y GALS) afIE—CD). 


Since our aim is to compare different quantities for the Hamiltonians 
I and I, it is natural that somehow we shall have to estimate the 
asymptotic smallness of the “correction terms” R, and Ry ; for this 
it will be necessary first of all to establish a series of operator inequali- 
ties. 

We shall show, in particular, that 


(5 »A\) (5 47) «< (5 ix) (5 4,47), (1.24) 


where the x, are complex numbers and the A, are operators. We have, 
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in fact, 


iA; + At) = x, AAt 
(54) (574) = Rava 


=} yi {x;xp AAgt +.x4x7 A, A}*}. 


But 
(xpAj— xj Aj) (x, Aj — x, At) = 0, 
i.e. 
xpXx,A Af +.xjx7 A At = 4x7 A, Aj + xjxf AjAf . 
Consequently, 


(54s) (tat) <5 (xe Adi +1 AeA) 


> [xul? 3 AAF, 


1 
k 
which is what we wished to prove. 

We shall now apply the inequality (1.24) to find a bound for R+ Ry 
defined using (1.21) and (1.23). We have 


RERy = (Wf) +! of) |?) (RFORP + RFORP) 


= REORO + RORY), (1.25) 
and also 
R,Rt < RORZ®+ RORE®, (1.26) 
Then, putting 
G, 


x, af), 4p = VG) Ja Ca) ay 


“VIG 


into (1.24), we obtain 


ROR = (FIG. ANE) (1G. I.- Ce) a-atfJE-C2)). 





But if we have an operator inequality 


A=<xB, 
the inequality 
UAU+ =< UBU* 
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is also true. Since 
a_ypaty< 1, 
we can write 


(J.—C,) pat (JE—C) < (Iu-C,) JE—CD). 
Consequently, 
RFORP <= P(f)Y |Ga! Ja—Cu) JE —CH). (1.27) 
In an entirely analogous way, we find for R®: 


RJORP = P(f)Y [Gal af (Jt —CD) J.—Ca)ay- 


With the intention of using the situation when a} and a, are found 
to be “squeezed” between (J,—C,) and (J+—C*), we shall have 
to interchange these factors, since a7 commutes with J, and a, with 
J}. For this, we note that 


(Jt —C2) (J.—C,) a (J.— C,) (J3—C)+IEI.—-ISt 


= (G.-C) ace IAAP (a,a¢ — af-a,) 


<= (J.—C,) JE —Cz)+ ae AAP, 


since 
1 
JEDI Jt = ce alae + a_y + arat — af a,) 
=a y |A.(f) |? (arat — a if ay) = ape we? -» |A.(f) |? bs 


Hence follows 
af (Jt —Ct) (Ja—Ca)ay 


<a}(J,—C,) (Jt -—CDay +a — yA. Sf)? 
= (J.-C) af a(Jz — Ca) + ey 


=< Ue-C_) UE- C+ gp8 FIAUNP. 
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Thus, 
R}ORP < P(f) iz IGa| (Jn—Ca) JE — Cz) + rg 2 1p) 


<PI{E IG Geez —Co+ 9 


Therefore, taking (1.25) and (1.26) into account, we obtain 


BPR <2. Gul Ge ~C,) (Jt—Ct) + a (1.28) 


We also obtain, entirely analogously, 
RR} < 2P(f) {x ieig.-coge-a +” wh}: (1.29) 


Hence, on the basis of conditions I, we see, in particular, that 


(RR) )r< PN fer+ el |< 2M, feet), 
(1.30) 
(RR )r = 2P(f) {ev v}< 2M, {ert 7 |. 


§ 3. ADDITIONAL INEQUALITIES 


Before we can begin to establish bounds for the differences of the 
averages taken over the Hamiltonians [’ and [°,, we shall have to 
obtain further inequalities for expressions of the form 


(By -»- By, RE RyBE -»- Ber, (1.31) 


(Br, «+» By Rp RABE». Bes 
in which 


+ 
By, => Oy, or ay, . 


Because of (1.28) and (1.29), this problem reduces to the problem of 
obtaining bounds for 


Da = (Bp, «++ By (Ja—Ca) (IEF — CIB}, ... Bar 
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We shall move the operator (J,—C,) successively to the left. We find 
Dz = (Je-Ce)Bp ++» By (JE —CaBR, .. Bie 

+ ¥ BAUE- CORE... Bide. 


where the operator B, is obtained if the operator f, in the product 
By, «+» By, (1.33) 


BiJe—IeBr- (1.34) 


To calculate this difference, we shall examine two separate cases: 


(1.32) 


is replaced by 


+ 
By, = and By = ay,- 
We have 


1 
A Ja— Jury = a5 » Af ){ayaf at ,— af at pay} 


I 
= 57 y Af) {[2,,a7 +af a) at,—af[agaty+a * pay ]}. 


(1.35) 
But, in view of (1.18), 


af = uf) af —o°*(f)a_y, 
at, = u(f)atpto"(f) a. 


It is clear that in the sum (1.35) the only terms which remain are those 
in which f = f, or —f = f;. We obtain 


ss 





Op Ja — Taty, = u(fi) avy, 


and also 


off Ja— Janz, =~ AG) 


vf) a, 
We have, consequently, a bound for the norm{ of expression (1.34): 


|ByJa—IaBy,| =< ACD ’ 


+t Everywhere below, |...| denotes the usual norm of the relevant operators 
in the second-quantized representation in the Hilbert space 2. 
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Since the norms of all the operators 8, are not greater than unity, and 
B, is obtained from the product (1.33) by replacing By, by the expres- 
sion (1.34), we see that 


(B,Bt)p < LAC 


Making use now of the general inequality 


|(BO)r| <« \W(BBY\ rn (O*OYr, (1.36) 
we apply it to formula (1.32). We obtain 


Da V(Je-Ca) Br, «+ BBA «+ BRUE —CD)r V Ds 
+ $ MD yo, 


j=1 
Setes te ak Sarr a, | fee eG fox 
< VU CVE Cir VB.+ & IO! 5, 
For brevity we put 


(Ja C,.) (J as C)r = A,. 
We write 





7% 7, & IACI 
Dy <= V4. + “, 
JD. <= Vv. Ley 
and 


Dz < Aq+2V/ Ay La [Aad +(5 Lappy” 


Applying this inequality to the averages (1.31) and taking (1.28) and 
(1.29) into account, we obtain 


(By «++ BRP R Ba +. Biyr < 2 2 {TI0. | Ae 
repre fs piel (EF) 
< 2P(f) 24 2P(f) fer2 var E1o($, BOI) 
+5l6 (i (20)) ‘ 
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But 
s 2 s 
(X|A0nl) <s 3 DP. 
Hence, 
Aa 2 
21G.| (3 a vi PD = pan 
Thus, 


Br «++ BRP RB... BED r 
ann {it +(Verty vin) | (1.37) 


<2 ly +(Verty VM.) | 
and, analogously, 
Br «+ BaRe RF BA, «+ Bhd r 
= 2m, {Mts (var+3 vin) | 


Having established the inequalities (1.30), (1.37) and (1.38), we can 


now turn to the proof of the asymptotic closeness of averages for the 
Hamiltonians J’, and I’. 


(1.38) 


§ 4. BOUNDS FOR THE DIFFERENCE OF THE SINGLE-TIME 
AVERAGES 


First of all we study the question of “single-time averages” and 
obtain asymptotic estimates for differences of the form 


Br ++» Bar— (Bp --- Barre: (1.39) 


We remark that these differences need be examined only for even 
values of s. For odd s, both terms in (1.39) are identically zero. In 
fact, both the Hamiltonians I and J’, are invariant with respect to the 
canonical interchange of Fermi amplitudes: 


ay > —ay, aj +—at, 
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i.e., in the amplitudes «, 


far: + —at 
ap Oy» ai > Of. 


(By, --- Barr. = (— 1) By --- Badr 


from which our remark follows. In the determination of estimates 
for the above differences, the spectral representations’ for two-time 
averages are of fundamental importance: 


Therefore, 


+00 
(A(t) B(t))r = f J 4, p(@) eC-9 da, 
oe (1.40) 
(B(t) A(t))r = f Ja, p(o) eeXt—9) den, 


where A(f) and B(r) are certain operators in the Heisenberg represen- 
tation, defined as functions of time by the equations of motion 


i 40 = A)H-HA(); A(0) = A, 
i BO = B(x)H—HB(t); BO) = B. 


We remark, incidentally, that we shall use these spectral representa- 
tions only for a fixed volume V; in this case, the Hamiltonian I has a 
discrete spectrum. This is ensured by the conditions I.4 (cf. Chapter 
1,§ 1). 

We denote the complete set of its eigenfunctions by y, and the 
corresponding eigenvalues by £,. 

Then, as one can convince oneself, 


(A(t) B(t))r = >> (G+, APu) (Pus Bor) e~ Fv + Er Ew C—9), (1.41) 
where Q = Tr e—#/, and 
1 
Ja, 3(@) = roy » (%, Agu) (Gu, Bg,) e~ Eee 6(E,—E,—o). 


+ See Bogolyubov and Tyablikov, 1959; Bonch-Bruevich, 1956 a, b; Zubarev, 
1960; and Bogolyubov, Jr. and Sadovnikov, 1962, 1963. 
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Consequently, the integrals occurring in (1.40) are found to be discrete 
sums. 

It is not the detailed structure of the so-called spectral intensity 
function J, ,(@) which will be of special interest to us, but only a 
number of properties, which we shall use constantly, of such spectral 
functions. We have in mind the following general relations: 


J4,a+(@) = 0, 
+00 
f Ja,+(@) do = (AAt)r, (1.42) 


—°9o 


+00 
f Ja,a+(@) eda = (At A)r, 


f |hy(w) how) J, a(w)| doo 


O i 
< / f lAx@) 2 F4, a+(@) deo J [holo |? Jas, alo) deo 
1) 
for the continuous functions h(w) and ho(w), where (/) denotes any 
particular infinite integral. 
We shall give special attention to the spectral representations 
+00 


(af (D) a(t))r = J J, ie «(®) eimt—2) doy, 


+00 


(a(t) af (1))r = f J, a -) €29 giolt—1) doy, 


(1.43) 


+00 


(RIOR) = J Tas, alo) eet? deo, 


+00 


(R(t) RF (Or = f J Rf. ni) 2/9 giolt—t) dey, 


(1.44) 


Here, R, and R7 are the same as in eqns. (1.20) and (1.22). 


RH) =-1 FO_ Bap, 


jo) 





R(t) =i —E(f) a,(T). 
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It is not difficult to see from this that 


dea ta (1) 


Rt t e—iE(f)t =_ 
(1) at ’ 


, det®V)*op(x) 


iE(f)r — 
R,(z) e' om 
Therefore, on the basis of (1.44), 
+00 
de—EM tat(t)\ (de® tax) ‘ 
(o-E(S) (t—r) — ie ee Se pc ee er AD 
[eee a |g) 


0° 


We integrate this equality over ¢ and t in the intervals (t, t+) and 
(t, t+) respectively, where 4 is an arbitrary number. We obtain 


+00 
Je(o-EY)) 4_ 4 |? 


| Paral (0 Ey 


= (Le EN Det(t + A)—e EW) tat (t)} 


X (eB + Dee + A)— eM Fat). 


ei(e-EV)) t—1) day 


But the right-hand side of this equality can also be determined from 
the formula (1.43). As a result, we obtain 


| ei(w-E(f)) 4 4 |2 Recep onas 


+ 
J @) ————___—_ 
| Rt Rf ) (o— Ef)? 
+00 
= | J.4 (a) [ele- 2) 4— 1 [2 e-FN)C-9 den, 
Oy, Oy 


Since this equality holds for arbitrary values of (f—t), we see that 


Trp, 26) = Jug, (0-EDY. 
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Therefore, from (1.42), 
+00 
J Jap, (o) (@- EY)? deo = (RER,)r. 
+0o 
Tap, (to) e°%(0—E(f))P deo = (RyR})r 


Taking notice of the inequalities (1.30), we see from this that 


+00 
[2 oe @) (w—E(f))? dw < 2M, fer +> ah 
< (1.45) 


+00 


J? asa) e*9(o— E(f))? dw =< 2M, {ert ar} 


Following the basic idea of the method described in the paper by 
Bogolyubov, Jr. (1968e), we now evaluate the difference 


D = (af B)p— e- = Bat) p, (1.46) 


where B is a product of an even number of operators Brp arranged 
in “normal” order, i.e. with all the a to the left and all the «, to the 
right, 

B= 1 Of nea Oy Bee 


and establish the asymptotic smallness of this difference. 

We note, first, that we shall need to consider only the situation 
when no two fj are the same and no two fy are the same, since other- 
wise B is found to be identically zero. But in such a situation, if we 
take into account that the total number of operators 8, in B is odd, 
at least one of the following statements is true: 


1. In the operator product B there is one operator a; such that 
the index f; is not equal to any of the indices f of the other operators 
8, contained in B. 
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2. In the operator product B there is one operator a, such that 
the index f, is not equal to any of the indices f of the other operators 
B, contained in B. 


For example, let statement (1) be true. Then we can write that 
B=a,,W, 
where W is the product of the remaining 8, in B. 
W=... By... (1.47) 


Since the number of operators occurring in W is even and their 
indices f all differ from f, we see that a must commute with W: 


B=07,W = Way. 
Taking the conjugate, we also find 
Bt = Wray, = a,Wt. 
On the other hand, it follows immediately from (1.47) that 


Wt*W<1, WWt <1. 
Therefore, 


BYB = a,,W* Wa, <ayaj,» 
BBt = ay, WW ay, KOs Oy, . 
Hence, 
BtB+BBt <1. 


We can convince ourselves of the validity of this operator inequality 
in case (2) in a completely analogous way. Thus, 


(B+B+ BB*)p< | 


always. Hence it follows that 


7 Jp+, p(w) (1+e°!*) do << 1. (1.48) 


—0° 
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We turn now to the quantity (1.46). Using the above spectral represen- 
tations, we write 


+0°o 
D= J Tay, (@) (1—ele-FN9} deo (1.49) 
and 
+00 


ID = f 





Jus, a()| |1—elo-FN9| doo. (1.50) 


We note here that if 


wo>E(f), 
then 


|1—el@-EUV9| = elw-EWV0_ | < eat elo EN 


ae (o-E(f)) e/9 
gonial ara , 


since, by definition, 


E(f)=0 (1.51) 
always. Further, if 
wo < E(f), 
then 
| 1—ele-2U8] = 1 — e-{EU)-0)/9 ta : 
Thus, 


\1 — elo ES} | < eae (al +e!) 


always, in view of which we obtain from (1.50) 


+00 


\Di<4 { 


Tap, go)] |o—E(f)| (1+ 60%) deo. 





THE MANY-TIME CORRELATION. FUNCTIONS 47 


Hence, using the inequality (1.42), we obtain 


SS 
oma) J Taz, (00) (0 E(f))P (+8) deo 


«|| [en torasenan f Joe, plo) (1-+62!) deo. 


Bearing in mind the inequalities (1.45) and (1.48), we therefore obtain 


D1 <5 Vf miler + Me), (1.52) 


§ 5. REMARK (1) 


This bound was obtained by starting from the trivial property 
(1.51). Since 6 occurs in the denominator of (1.52), we cannot use 
this bound for treating the zero-temperature limit 9 ~ 0, even if 
ty > 0 (VV — ~) uniformly with 6 — 0. 

Let 

Ef)=y, (1.53) 


for the (f) under consideration, where y is a certain positive constant. 
Since 


Bf) = (F-a) +14r. 


we see that the inequality (1.53) is always satisfied outside the spherical 
layer: 





We turn now to the inequality (1.50). We note from what has been said 
above that 


[1 —ele- £0 | jon BDI el e((2/3y-EN)/0 
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for 
w = E(f)—2y. 
Also, 
|1—ele- AN) << | 


for the values 
wo = E(f)—Fy. 


Thus, from (1.50) we find 


IDi= | |Jup,(@)|10-E(P)| ev" do 5 e-0 





dw. 





PRE sere Juz, (oo Elf )}re0" deo 


E-938 4-0 © 
al J 3+, Presrerane]/ | Jy a, ©) dw | J 3+, p(@) do. 


But it follows from (1.45) that 
E-2y/3 


Me 
Je, + (@) do = Qy 5 2M(er+ aw) 


Therefore, we obtain 
Me 
1D ={5 ent BB r+ > VEB;| |/2m(er+2#). 


Hence, taking into account that 


(BBt +B*B); <1, 
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| ray tt er sae 
=< —__ p— 27/36 panne eee 
\Di<\/ he 2y +(5,) Y/2ma(er+ am) 


But it can be seen that 
2 2 
x2e-* <= | — 
(<) 
and therefore 


1 wy\? of 3\2_ /2\2/3\2 
Re 27/30 — (33): e ad er < (z) (37) . 


we find 








Thus, finally, 
3 
ID| < ms y/2+2(2) *Y o(or+ M2) (r+ wv) (1.54) 
we combine the inequalities (1.52) and (1.54) in the form 
M, 
|D| + 4? (1.55) 
where the constant X can always be taken in the form 
; 
K=5 VM. (1.56) 


In the region of (f) for which the inequality (1.53) is fulfilled, we can 
also put 


z 2M,4+2 Ss *M ‘a 1.57 
25 i+ (z) i , (1.57) 


This latter value of X clearly no longer depends in any way on the 
temperature 6. Thus, we can write 





at B)p— e-2WO( Bat) n| =X err. 


Hence follows the inequality 


| +e-20) at B)p— e- EUW (at B+ Bart) p| =< K [err ‘ 


50 A METHOD FOR STUDYING MODEL HAMILTONIANS 


and therefore 


ise e-EUNy/0 r a 
(af B)r— T+e-EG ye (af B+ Ba} Dr 


— ie 
<= xX err ie . 


From this inequality we shall now be able to obtain the required 
estimates for the difference (1.39). 


(1.58) 


§ 6. PROOF OF THE CLOSENESS OF AVERAGES CONSTRUC- 

TED ON THE BASIS OF MODEL AND TRIAL HAMILTONIANS 

FOR “NORMAL” ORDERING OF THE OPERATORS IN THE 
AVERAGES 


We shall first of all examine an operator product of the form 


+ + 
U = af, 6. Oh BH +. OH 
in which 
k>qz2=0;  k+q is aneven number. 


We assume, of course, that all the f; (like all the f;’) are different, since 
otherwise U will be identically zero. Now, as soon as the number of 
creation operators «* in a given product is greater than the number 
of annihilation operators «, there is certain to be amongst them an 
operator aj such that the index f; is not equal to any of the indices 
fi, »» +> fy. But then we can write 


U= ay, B, 


where B is the operator product of the remaining £,. The essential 
point is that the number of operators 8, in the product B is odd and 
that their indices all differ from the index f,. 

Therefore, the operator Xp must anticommute with the product B: 


oj, B+ Boj = 0. 
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Using the inequality (1.58), we obtain the inequality 
M: 
Mat... apa, ... ag) r| < 2 er, (1.59) 
k>q=0; k+q even. 
We now take the opposite case. Let 


+ + 
U, = Of oe. OF, MY vee Offs 
q>k=0; k+q even. 

We have 


Uf =a... Of cep, oe 
and for Ut we can again make use of the inequality (1.59). But 


|(Us)r| = |(U)#| = |(U*)r|, 


and therefore 





[aft ... aay... wr =< XY eve Mt, (1.60) 
q>k=0; k+gq even. 
We now consider the case when k = q: 
Uszaz ... aay ... ay. 
If amongst the indices f,, ..., f, in this product there is at least one 


index f, which is not equal to any of the indices fj, ..., f, we have, 
analogously, 


U=afB; oafB+Baj =0, 


and from the inequality (1.58) it follows again that 


‘Me 
| aeft ee ony, orp ats a: r| =X / EV ay (1.61) 
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(if any one of the indices f,, ..., fj, is not equal to any of the indices 


ieee 9? 


Since in the “new” Fermi amplitudes the trial Hamiltonian I’, has 


the form 
= » E(f) afayt+const, 


we see that the averages examined above, constructed on the basis 
of the Hamiltonian I’ in the cases (1.59), (1.60) and (1.61), will be 
equal to zero if they are calculated on the basis of the trial Hamilto- 
nian I’. 

Consequently, the inequalities (1.59), (1.60) and (1.61) can be 
rewritten in the form: 


Ke det dra = afer, (1.62) 


where the averages are of the products treated earlier. 
Thus, of all the products 


Br, --- By (seven) 


in which the £, are arranged in normal order (i.e. with all creation 
operators «; to the left and all annihilation operators a, to the right) 
we shall only have to consider products 


+ + 
Of, ++ HFM... OF 


in which each of the indices f,, ..., ff, is equal to one of the indices 
fis -+-» Sq. But such products can obviously be reduced to the form 


+ + 
Lay... Of Oy, 16. OH. 
We put 
Ox = af, .+- Ope, ... %% = Of B, (1.63) 


and note (all the indices f,, ..., f, are different!) that 


+ +t + + 
Of, Bet Bray, = Hf, ey, 0, + OO, 6 OH, 
++ + = 
HOR Hf 6 BAO, 6. OH, = On-1 
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Therefore, from the inequality (1.58) we shall have 





<x ori, (1.64) 





| \ 
| eae ‘14+ eFhye (Ox-1)r 


and, in particular, for 


+ 
Or1 = a, a, 








we find 
ten e-em = KV cyt ae (1.65) 
But, obviously, 
teary = (apap dra: 
Also, 
(On)r.— Tae (Ox-1)r, = 9. (1.66) 


Therefore, it follows from (1.64) that 


(Oidr—(Ordr.~ <rgtagw (Or-)r—Oe-1)n} 





M2 
=X eV ay * 


Hence, transforming, we have 


Onr— (Opn) < 5 (Ona) (Oe-)r, 
(1.67) 
+ XK eyt+ 3 2 


It follows from the inequality (1.65) that 





|(Or)r—(Or)r, | « X Yet at , (1.68) 
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and from the inequalities (1.67) and (1.68) that 


Orden KY eve hia 5+... (5) 
kr kr, | = vray ger ee 2 
cea 
=<2X evt-ap - 


Thus, to summarize the above statements, we have shown that for 
products 


By ++ Br, (1.69) 


in which the operators 8, are normally ordered, the inequalities 





\ (By, shee By r— By, Peers By, r,| < 2X Yor (1.70) 


are valid. 


§ 7. PROOF OF THE CLOSENESS OF THE AVERAGES FOR 
ARBITRARY ORDERING OF THE OPERATORS IN THE 
AVERAGES 


We shall now consider the case of products (1.69) in which the 
operators B, are arranged in arbitrary order. We note that in the 
system of indices f;, ..., f,, some of them can be equal. 

From the aggregate of indices f,, ..., f,, we separate out all the 
distinct indices and denote them by fj, ..., f,, $0 that each (one or 
several) of the indices f;, . . ., f, is equal to one of the indices f/, ..., fy. 

We shall represent the operator product under consideration in the 


form 
By --- Brs (1.71) 


where B,, denotes a product of operators «,, az, all with the same 
index f’. 

However, by virtue of the commutation properties of the Fermi 
operators, it is not difficult to see that B,, can always be reduced 
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to one of the following forms: 
By =tayp, Bp =tap, Bp =tafay, (1.72) 
Bp =+ayaf =+(1—apa,.), (1.73) 


leaving aside the trivial case when B,, is identically zero. 

We shall reduce the products (1.71) to normal order, commuting 
the operators « only to the left and the operators «,, only to the right. 
The number of terms will not “multiply” as a result of the commuta- 
tion operation, since all the indices fis bere) te are different, and B,, 
consists of a sum of two terms in the cases (1.73) only. 

Therefore, the number of products of normally ordered operators 
into which the original product breaks down will be 2%, where q is 
the number of B,, having the form (1.73). 

But, obviously, g = s/2. Therefore, by making use of the inequality 
(1.70) for each of the “normal products”, we shall obtain for the gene- 
ral case under consideration: 


\(By --- Badr—(By --- Bera] = 26941 vta (1.74) 


To summarize the arguments given above, we shall formulate a 
theorem: 

THEOREM I. Let the conditions 1 (Chapter 1, § 1) be fulfilled. Then 
the bound (1.74) will be valid for the difference of singletime averages 
constructed on the basis of the Hamiltonians I’ and I,,. 

It must be emphasized that in this bound we can always take the 
value (1.56) for X: 


> ann 
— rf / M3. 
But then the estimate (1.74) does not allow us, in general, to pass to 


the treatment of the case 0 — 0, even if 


ey ~ 0 as V +o 


uniformly with respect to the temperature 0 in the interval 0 < 6 < 4, 
where 4p is any fixed temperature. 
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On the other hand, from what we said earlier we can use a “uniform 
estimate” for X% by taking the value (1.57) of the constant, 


2 / 
K= oa f2mi+2(z) iat * 
2y e 


if all the values of f, fy, ...,f, lie in the region in which 
E(f) > y. 


REMARK (II) 


We note, in passing, that it would be incorrect to think that the 
difference (1.74) in fact always tends to zero (as V ~ o) uniformly 
with respect to 6 (@ > 0), even in the region where E(f) = 0; the 
absence of this uniformity is due entirely to our method of treatment, 
which permits excessive majorization. 

To understand intuitively the reason for the absence of uniformity 
with respect to 6 in the region where E(f) = 0, we take an extremely 
simple (one could even say trivial) example: 


P=f,4+l; 
where 
Py = E(f)afay,; a =|£.- . 
Y) m 
I, =¥ eafay; e>0. 
Here, : 


e>0 as V+o, 


It is clear here that binary averages calculated over the Hamiltonians 
I’ and I’, have the form 


1 
(afoy)r = 7 + ele + 6/0” 


1 
(af: ar, 4 + eb" 
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We take some value f = fo such that 


E( fo) = Poy = 0. 


Then 


ee ee oe I I 
Capacr, rr —Kaerhanp Ta = TH gel ~ 3° 


Putting here 
6=2e>0, 
we find 
1 


| 
(ajraers rn — (ajeers re = EE 


Thus, even in this simplest case, we see that the asymptotic relation 


(apar)r—(ataryr, + 0 


is not uniform with respect to 6 ~ 0 in the region E(f) = 


§ 8. ESTIMATES OF THE ASYMPTOTIC CLOSENESS OF THE 
MANY-TIME CORRELATION AVERAGES 


We turn now to consider “many-time averages”; we shall estimate 
differences of the form 


(Bplts) .-» Belts)r— (By (tr)... Bylts)) re: 


Here we use the following conventions: 
1. Operators 8,(f) occurring in the average 


Coos 


are determined by the equations of motion for the Hamiltonian I" 
with the initial conditions 


(0) = By. 


2. Operators 8,(t) occurring in the average 


(My 
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are determined by the equations of motion for the Hamiltonian I, 
with the same initial conditions. 

We note that operators 8,(t) for the trial Hamiltonian I’, can be 
represented in a form with the time-dependence given explicitly: 


Brlt) = Bremen, 


of) = Bf) if B= ay, 
of) = Ef) if B= ap. 


Assuming that the conditions I are fulfilled, we establish the validity 
of the inequality* 


(Brats). - By(ts))r— (Bata) --- ByXts)>r,] 


<7 Syne J+ 26 +1K Vert M2 
Vv >» JT G+ 4V ’ 


where 





(1.75) 


where 





ae oo -_ ___\2yu2 
i) = 2M, {at 4 (Ver+ uv vm) (1.76) 
(j= 1,2,3,...)- 


We note, first of all, that since neither of the Hamiltonians and I, 
depends explicitly on time, both the averages considered will be 
invariant with respect to time translations 


t+ ¢t4+t 


with arbitrary t. Taking t =—t,, we obtain a situation in which f, is 
replaced by 0. 

Therefore, the inequality (1.75) will be proved to be true if we prove 
it only for the particular case t, = 0. Taking into account the nature 
of the time-dependence of the operators 8,(¢) for the trial Hamiltonian 


t On transforming to the old Fermi operators, we prove an analogous inequality 
by formulating a theorem (cf. Theorem II, p. 64). 
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I’, we see that we must prove the inequality 


I(Ba(ts) -.- By,s(ts—1) ByO))r 


—e-ifo(fat ..- +o(fe-ea} (By Sans Byor.| 
<P {|ti-tal+ . aes (1.77) 


+ 20m41Y > ee 





But 
(Bp Ct) «-- By,s(ts—1) By,O))r 
— e—HoUnat oon +ofe- Mead (By saad Byr,| 
= [(By (ts) --- By-a(ts—1) By, O)rea+ --- +o) 
—(By --- Badr. 
= (Bits) «.. By alts—1) ByO)r— (Bp, «-- Bydr.|> 
where 


BAO) = Ben, (1.78) 
and, therefore, the inequality (1.77) is equivalent to the inequality 


| (Br(ta) .-- By (ts—1) By (O))r — (By, --- By ra! 
=P {|ti—te|+ ... +(8—2) |ts-2—ts-11+(s—1) | ts-1|} 


(1.79) 
[ M 
+QWMHK Y vt ae ‘ 


It is not difficult to see that this inequality (1.79) will be proved as soon 
as we have established that 


Bats) «- » By alts—1) By O)r— (BO) -- - Br,-«(0) B,(0))r| 
<= {\ti—telt+ ... +(8—2) |ts-2—ts-1|+(5—1) |ts-il}. (1.80) 
In fact, if the inequality (1.80) is true, the left-hand side of inequality 


(1.79) will be less than the left-hand side of inequality (1.80) plus the 
expression 


|(By(0) ... By,(0))r—(By, - - - By) rel 
= |(By --- Byr—(By «+» Brdra|> 
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which, by virtue of the inequality (1.74), will be 























< 20941 VY ey + ca : 
To establish the inequality (1.80), we note that 
KXBa(ts) «+» Br-a(ts—1) B,O)r—(Bp(0) - - - Br,-s(0) B,))r| 
bs {(Br (ty)... Bt Bryulten By,(0))r 
= = = 1.81 
—(Baltian «+» Byltian) Bra ltras) «-- By O)r} | ve 
y+. A ae 
a / aBy(t) ... Bl \ 
= d On SET IN ZONE ‘ 
& : ‘ dt a 
where ; 
oA; = By,.,(tj41) -- - By, (0). 
In view of this, 
Aj A; <l. (1.82) 
But, 
tytn 4 re 
| at - «++ By) ,) 
rT 
4 ty 4a 2 { 
d 
< J {iio HD 4,,) |, 
where 
tgs = Brad --- By (ot; 
and, from (1.82), it is clear that 
Therefore, 





_ apy (t 
Ce 


rd r 


e dB,(t) dBt(t 
«| Gino... By(t) dBe(t) 


ener ares into), 
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Hence, taking inequality (1.81) into account, we obtain 


Balt). Bya(ts—s) By(0))r — (Br,0) - - - By,-(0) B(O))r 


tg+1 


s— : a aa 
Ds | Cato stake Br ) Pat y 








oe \ dt Bo), dt|. 


(1.83) 


We turn now to the equations of motion (1.20) and (1.22) for the 
Hamiltonian I and represent them in the form 


ao (t) 


uy] 





+io( Br) = ZO), 


in which 
Zp=—iRy if Bp =ay, 


Zp=iRt if Bp =at. is 
We therefore obtain for the functions (1.78): 
a) = Z, (Heo, (1.85) 


Consequently, we have: 
Cixo see ae “as iio), 
= (BAO... ZAOZEO «.. BEO)r 
= 8; (0) ... Z,(0) ZO) ... Bi O))r 
= Bi + ZyZig +> Bir 
We now take our inequalities (1.37) and (1.38) into account. Since 
ZZ} = R-Rf or RPR,, 





we convince ourselves that 


Chat BO 


(q) 
<=Ny <ny 





a 





dt : Auto ), 
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Thus, from the inequality (1.83) we find: 
Bilt) «+- By BaO)r— BO) - -- Ba(0) B,())rI 
< FS Pily— tail. 
A 


The relation (1.80) is now proved, and the proof of inequality (1.75) 
is thereby completed. 
We see from inequality (1.75), in particular, that the convergence 


Balt) --- Belts)r— Balt) -.- By (ts), + 0 


as V +o 


is uniform with respect to ¢,, f,, ..., ¢, in the region in which their 
differences are bounded :* 


I4;—tail <S. 


(Here S is any fixed constant that is, of course, independent of V.) 

Up to now we have considered averages composed of products 
of operators «,(?’), a(t’). We turn now to examine products of the 
“original”, “old” Fermi operators a,, a/. 


+ To elucidate intuitively the reason for the absence of uniform convergence 
with respect to f), ta, ..., ¢, in the entire space of the points f,, ..., ¢,, we shall 
consider the simple example when 


r= ret, 
with 
ro= PEN) a a, 


P= Deda, e>0 («+0 as V+ oo). 


We obtain, in this case, 


1 = 
(agt(ty) afte) = Tarr eM +8) (4-4 
= 1 $4) 
(af (ty) afte), = CO e aa) 


It is clear from this that these averages can be close to each other only so long as 
é€|t,—tel is small, 
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We shall be concerned with averages of the type 


(Ap (ty)... Ap(ts))r, (Ap (ti)... Ap(ts)>r, 
where 
A(t) =a,(t) or af(t). 


Here, the Fermi operators A,(¢) occurring in the averages 


(eciedpy Goce ry 


are determined by the equations of motion for the Hamiltonians J‘ and 
I’, respectively. 
We note that, by virtue of the relations (1.18), 


a,(t) = u(f) a(t)— o( f) «t,(2), 
af (t) = u(f)af()— v*(f) «_,(0).- 
We write these canonical u-v transformations in the abbreviated form: 


A(t) = un(f) By, (2) + ual J) By, of). (1.86) 
Here, 
lun(f) |? + | ual f)|? = 1, 
and, therefore, 


Y lu(f)| = V2. (1.87) 
y=1,2 
Using the relations (1.86), we obtain 


(Ap(ti) ... Ap(ts))r— (Ap (ta) -.- Ap (ts), 
= i, Wy) {Banta «~~ Brodts))r— (Banta) -- + Brrdts)ra}- 
tA ,,£ 


1,2 
Hence, using (1.75) and (1.87), we can convince ourselves that 
\(Ap(t1) «.. Ap(ts))r— (Ag (a) «-« Aylts))ra| 


s-1 M2 \ +3 
< ("5 YL Flt tral +2 éyt Mah I] » ef! 


lye, 


<2 [nf ? ps Flt —t rl FWO1X Yorr4 Mal 
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Thus, we have proved the following theorem: 


THEOREM II. Under the conditions I, the following inequality holds: 
| Ay (t1) Brae Ap(ts))r— (Ag (ti) ara Ay(ts))r, | <= 69) 





or = =e 2 1/2 54 
= Ys+12 MMs. vy, (Verse Vit) Y slh—-taa| 
av V jm 
ag 
$HK Yor Me. a 


Here we can always take the value (1.56) for X. If 
E(fi) = y, «++, EPs) > 


we can take for X the expression (1.57). 


CHAPTER 2 


CONSTRUCTION OF A PROOF OF 

THE GENERALIZED ASYMPTOTIC 

RELATIONS FOR THE MANY-TIME 
CORRELATION AVERAGES 


IN THIS chapter we shall obtain generalized asymptotic relations for 
many-time correlation averages, T-products and Green’s functions, 
starting, essentially, from the results of Chapter 1 and the theorems 
which were proved there. The systematic derivations for finding the 
generalized asymptotic relations are summarized by theorems. 


§ 1. SELECTION RULES AND CALCULATION OF THE 
AVERAGES 
We shall now study the “selection rules” for the averages considered 
in Chapter 1, § 8, i.e. for 
(Ay (ti) ... Apts). 
As has already been noted, for odd s these averages are identically 
zero. 
We shall now formulate more stringent selection rules. We take the 
operators 
+ a) aa 
Mp, = Aff, Nf = A_fpA_fy, 
and note that their difference 
(1%,—n_f.) (2.1) 
65 
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commutes, for any value of fo, with all operators of the form 
ata, Afar, apaty. 


Therefore, the operators (2.1) will commute with the model (/°) and 
trial (”,) Hamiltonians, and will consequently be integrals of motion 
for them. 

We introduce the unitary operators 


Us, = n—"-1), Uj = eG —"-1), 
where ¢ is an arbitrary real number. By virtue of the above arguments, 


TU,, = U,,0 
and, therefore, 
e- TU, = Uy,e-79. 


Consequently, we can write 


Tr {Ap(ty) ... Ap (ts) 07 7U, UF 
(Ay (ty)... Ay (ter = SEASON Arena) 


_ Tr{Uf g(t)... Ag(ts) e7?U;,} 
~ Tr e-76 


_ Tr{Up Apts)... Aplts) Uj} 
~ Tr e-T@ 


= (Up Apts). Ap(ts) Up)r- 
But in view of the unitarity, we have, identically, 
CUP Apt)... Apts) Up>r 
= {Uf Ag tr) Up} {Uf Alte) Un} «- {Uf Arts) Upghr. 
Thus, we have proved that, for any value of fo: 


(Ap (ti)... Ap(ts))r 
= {UZ Apt) Up}... {UR Apt) Up})r (2.2) 
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We note now that 
UprapUy, =e" ay,,  — Uppa_pUp, = €"a_g, 
Ujat,U; = eats,  UfapU, =e “af, 
Uj,yU, =a, if fx fo, —fo, 
UjafU,=af if f# fo, —fo- 


(2.3) 


ge the operator n,—n_, and, consequently, also the operators 
/,, and U;, are integrals of motion for the Hamiltonian I’, the rela- 
ms (2.3) will remain valid if in them we replace the operators a, 
and a; respectively by the operator functions of time a,(t) and a; (?), 
which are determined by the equations of motion for the Hamiltonian 
r. 
Having noted this, we turn to the equality (2.2). We denote by 
N (fo) the number of operators equal to a, or at, amongst the 
operators As, ae ee and by N_(fo) the number of these operators 
equal to a7, or a_,. The right-hand side of (2.2) is then equal to 


ellN+-N-Uol? (Ay (ti)... Ay (ts))r 
and we obtain 
(1—elN+f)-N-} 9) (Ap (ty)... Ap (ts))r = 0 


Hence, in view of the arbitrariness of the phase g, it follows that 


(Ay,(t1) oer Ar(ts)) =0 
if-the following inequality holds for any index fo: 
Ni(fo)—N_(fo) # 0. (2.4) 


Two operators A, and Ay, from the operators A,, ..., Ay, will 
be said to form a pair if 


1. for Ay, = 4,, Ay, = aj, or @_y,; 
2. for 4,=a7, A, =a, or ary. 
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It can be seen that if the operators 4, and 4, form a pair, then 
I, = th 

Let there be among the operators 4,, ..., dy an operator Ay, for 
which it is impossible to select an operator A, forming a pair with it. 
In this case, we shall say that 4, is an unpaired operator. 

We shall show now that if among the operators 4,, ..., A, there 
is at least one unpaired operator, then 


(Ap (ti)... Ap(ts))r = 0. 


In fact, let A,, be such an unpaired operator. There are two possibili- 
ties: 


1. A, = qy,; 
2. A, = aj, . 
In the first case, none of the operators 4,, ..., Ay is equal to a; or 


a_y,. In view of this, 
N_(f)) = 0. 


On the other hand, 
NiGf)=1, and NiG)-N_(fj) = 1. 
In the second case, 
N-(f)=1, Ns) =0, and Ny(f)-N-(f) <—1. 


Thus, by virtue of (2.4), our assertion is proved. We see, consequently, 


that the average 
(Ag(t) «Ay (ts))r 


can be non-zero only if all the operators Ay, Ap, ..+, Ay, can be 
grouped in pairs. 

We have just considered averages for the Hamiltonian I’. From 
the arguments given, it is clear that the results obtained are also valid 
for averages taken over the Hamiltonian I,. We further note that, 
since the trial Hamiltonian I’, is quadratic in the Fermi operators, 
in calculating expressions 


(Ap (ts)... Ap(ts))r.> (2.5) 
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we can make use of the generalized Wick procedure established by 
Bloch and de Dominicis (1958, 1959). In accordance with this proce- 
dure, the average (2.5) is equal to a sum of products of binary averages: 


LT] (4) 4,4 oN = tl, fy =f. (2.6) 


The summation here is performed over all possible pairings in the 
product A, ... A,. The number of pairs in the sum (2.6) is equal, 
consequently, to the number of ways of combining the operators 
Ay, -++, Ay, in pairs. This number is obviously less than 


3... (s—3)(s—1), seven, s > 2. 


The binary averages occurring in.(2.6) are calculated directly by means 
of the transformations (1.18). 
We have the following expressions for the binary averages: 


A eB I-#) “1-9 
(af(t) a_f(t))r, = OH ecco = Tem 


1 /, T(f)\ e#ne-9 
(a) ()r, = ( -ay) Trek 


T(f)\_e-inenn 
+3 (+E )) Tree 


iE(/) (t—¥) 
(af (4) af(t))r, = (+H) Tae 


1 T(f)\ e-#WU-9 
+2 (IE) treme 

A ) efBV) (t—-2) eB) (t—1) 
(a} (at At))r, = — ah Trees ~ Tp e- ee 


(2.7) 


Thus, the average (2.5) is a polynomial in the expressions (2.7), in 
which we replace f by f,, and ¢ and t by ¢, and ¢,. 
This polynomial 


ae 1; ere fe th, aerigihs)s (2.8) 
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is defined in the whole space ¢ (i.e. for all points (fi, ..., f,)), since 
the functions 


Af), TH), EP) 


are defined everywhere in ¢{. 

If the function A(f) does not depend on V, the above polynomial 
is also independent of V. 

On the other hand, the equality 


(Ay(t1) poets Ag(ts))r, = PV dyes shes ti, ..., by). (2.9) 


holds, obviously, only for the quasi-discrete indices 





- 2 2an, 2nN2 2nn3 
St =(p,9), p=( L a aa a | 


i.e. only on the set ¢%, since only for these indices (f) do the actual 
operators 
Ay => a, at 


have any meaning. Thus, the dependence on V of the above averages 
over the Hamiltonian I, (in the case when A(f) is independent of V) 
is due entirely to the fact that for (fi, ..., f,) we take the points of the 
set 6. 

§ 2. GENERALIZED CONVERGENCE 


We introduce the concept of convergence for functions 
Fy(fi, .-..fs) (WV >) 


defined on the sets $9. The function F\(fi, ..., f,) will be said to 
converge uniformly (as V + -) to some function F(/i, ..., f,) defined 
in the whole ¢, if one can find a sequence of numbers 


dy >. 0, V —~> oo 
such that the inequality 
IFA oo lhd—F(f, on Ss) <= by 
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is fulfilled everywhere on the sets ¢%. We note now that, by virtue 
of (1.88) and (2.9), we have 


(Ay (t1) ... Ap (ts))r—DA, -- +5 fe5 tay »--5 ts)! 


MM: 1 Ve 
i 7 vi) Y Jlt4—t+ril (2.10) 
4Vv j= 





aero 


M. 
+ 2+1K Vert ae. 


everywhere on ¢%. 

Thus, in accordance with the definition of convergence introduced 
above, the following rider to Theorem II can be seen to be valid. 
If the conditions I are satisfied and, in addition, the function A(f) is 
independent of V, then the limit 


lim (Ag (ty)... Ag(ts))r = DA,» +05 fe3 fay eee ts) (2.10) 


Me ui Ver eps 





exists. 

Up to now we have considered averages composed of products 
of Fermi operators a,, a. We pass on now to the operator field func- 
tions 


polr, t) = —= PY ay, (er, 
ve (2.12) 
yo(r, t) = phe it (Nee. 


Here, the summation runs, as before, over the set of quasi-discrete p. 

We shall be concerned with averages composed of a product of 
operator field functions, constructed'on the basis of model and trial 
Hamiltonians, I’ and I, 


(Po (Fr 1; f 1) eee Po {F 's> ts))r, Ts? 


where s is an even number (for odd s the averages are identically 
zero), and where 


PalF, j= Wal? t) or ye (r, n). 
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By virtue of (2.12), we can write 


1 
Polt, t) = VV yy Ap, Ader) , 


where 

Apo=4pc, €=1, if p=y, 
and 

Apo=@pc, €=—l1, if p=yt. 
Thus, 


(Poi, £1) «-» Po,(ss ts))r, Fe 
(2.13) 


1 
~ or, y > (Anos (t1) «+» Apgo(ts)>r, relator rot teleer a), 


We shall consider the difference 


(Po(Fiy t1) «+» Pa,(Pss ts) —Po(F1, £1) » ++ Po, (Ps ts))re5 


and show that for V > © it tends to zero, though now in the gene- 
ralized sense customary in the theory of generalized functions. We 
recall here that in this theory, 


Fri, ...,%s) > F(ri, ..., Ps), 
F(ri, ...,%2) = lim Fy(r, ..., ts) 


(V — oo) 


(2.14) 


is called generalized convergence. We shall consider the class C(g, v) 
(where g and » are positive numbers) of continuous and continuously 
differentiable functions A(r, ..., r,) such that 


{Irilt ... +1rs|}/ |A(ri, ..., %s)| < const 
(j= 0,1, 2, ...,»), 


OUNt «+» +Qaeh 
{lril+ ... +1rsl}!| -=——~=— | <const 
Oxi! a ee re) 33, a, 


and 


(j =0,1,2,...,%5 git ... +qgss =90,1,...,43 oa = I, 2, 3) 


in the whole space E® of points (r1, ..., ¥,)- 
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Here, 
(1,1, X1,2, X13) = Fr. 


Then, if we can fix the positive numbers g and » in such a way that 
the relation 
fra, «5 ts) Fra, «+. ts) dry... drs 
> fr. .- ey Ps) F(ri, ..., rsddri ... drs 
(VV + ~) 

is valid for any function A(r, ..., r,) from the class C(qg, »), there 
is said to be generalized convergence (2.14). For our purposes, we 
shall confine ourselves to examining the class -2 = C(qi, »1) in which 


the numbers q: and »; are chosen in such a way that the Fourier 
transform of the function h(n, ..., r,) 


h(pi, -.-, Ps) = fac, wo ay Pselrenit +. +Pe°r) dry... drs 


is a continuous function of (pi, ..., p,) in the whole space E®, and 
satisfies the inequality 


1 eee (2.15) 
[I (ate 


X, = const > 0, 
My, = const > 0. 


We shall show that for any function A(ri, ..., r,) in this class 2 


far, +s Ps) Poli, t1) --. Po, (Ps, ts)r dri... drs 
-f W(ri, ..., r;) (Po,(F 15 th)... Po,(Ps, ts)r, dr, ... dr, +0 
as V+, 


Re-expressing this using (2.13), we shall prove that 
1 
pa, Alera, --- €Pa)(Anoilts) --- Apyos(ts))r 
Pir +++ De 


—(Apr(ti) -.. Ap,o,(ts)r,} ~ 0. 


(2.16) 
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Here the summation runs over the points (pi, ..., p,) of the quasi- 
discrete set E®. We shall thereby also prove that, in the generalized 
sense 

(Po,(F1, t1) --- Pols, ts)>r 


2. 
—(Po,(1, £1) --- Pots, ts))ra 0 as V+ GM 


§ 3. REMARK 


The question may arise as to why in considering averages of products 
of operator field functions we have had to confine ourselves to proving 
only the generalized convergence (2.17). 

The point is that, even at fixed finite volume V, the sums (2.13) 
representing the actual averages can, generally speaking, be divergent 
in the usual sense. 

We take as an example the expression 


(Polhi, £1) ys (ro, t2))r,- 
We have 


(Pelt 1, 41) pr (re, fe))r, = y, y (Bp ft) Bolle) rekon sro} , 


Hence, taking into account the formulae (2.7) and the selection rules 
established earlier, we obtain 


(polis 4) pe (Ya f))r, = ; 2, oty) 23, ota) rele 1-"») 


=a7E ied) Cul pe een 
Ww E(f)) Teen 


T(f)\ ew) 
+57 ( TE) yi pe eee 


For simplicity, we take the case when 


Af) = 
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identically. Then, 


Ef) = IT(f)| = ae | 
and we have 
(Walh1, th) we (re, ty))r, 
i —#\4- 
= 2 — seca i el(p+r—r2) 
V cprbmy<ny 1+ €! (PH2m)—x| (2.18) 
1 e~H(p?/2m)— 4) (4 —t) 


v ——,,——__- ellp-n—rs) 
V pramy>uy Lt em (e*2m)—a} 8 ellpsr—re) | 


It can be seen that for fixed V the first sum contains a finite number 
of terms and the second sum an infinite number. In this case, 


1 


1+ e-eam—)18 1, 


(Ip| + =), 


Consequently, the second sum is a divergent series (in the usual 
sense). 

It is not difficult to see, in passing, that the given series is convergent 
in the generalized sense. In fact, multiplying it by a function A(n, re) 
from the class -€ and integrating, we obtain the sum 


1 e—i{(p?/2m)—1) (4-2), 


V oben +e Rm 16 h(p, —p), 


which converges absolutely, since 


l 
pu ee eres oa 
Ve cprdmeny 1 tem (etm — 1/6 |A(p, —p)| 


1 dp 
eons Ly 
eos (p?+ Mr)? : 


x 


{(p2/2m) =") 


t(yv+O as Vo, 
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Thus, even in this very simple case, we see that the average 
(PolPy> 4) Ps Pa» t2))r, (2.19) 


is given by a sum (2.13) which is convergent only in the generalized 
sense; the expression (2.19) is, therefore, defined only as a generalized 
function of r, and re, even for finite V. 


§ 4. PROOF OF THE ASYMPTOTIC RELATIONS 


We turn now to the proof of the asymptotic relations (2.16). Using 
the notation 


(2.20) 
Dp,o,; ne eed (Ap,o,(t1) tae Apots))r —(A po, (ta) sae Apo fts))r, ’ 


we obtain, by virtue of (2.15), 


Kn z 1 
| Dh) lS Fae i x Downs ---svenl UT Caray: (2.21) 


We shall take into account here the above-mentioned selection rules 
and use them for Dayo: ...;p,,- Owing to these selection rules, the 
summation in inequality (2.21) need be carried out only over those 
Pi, P2, .--, p, for which the system of operators A,,,,,.--, Apc, Can 


be divided into pairs A,,,, 4,,,, such that 
Pro=p; if op =a% 
and Pr=—p, if of =—o%. 
We denote 
Pi), =, Phy = 2 -- +> Pity = Usi2- 


We have 
s/2 1 


s 1 
Nig — qr: 


Further, we denote the number of ways the operators A,,,, ---» Ap,o, 
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can be divided into pairs by /(s). Obviously, 
K(s) < 1-3 ... (s—3)(s—1). 
Then in the right-hand side of (2.21), the sum 


can be broken down into /(s) sums of the type 


ae Oe (2.22) 


Qty +++, Gort 


We note that the same terms can be found in different sums (2.22). 
To exclude these from certain of the sums (2.22), we must not perform 
the summation over all the (q,, . . ., 4.) of the quasi-discrete set under 
consideration; a number of “exclusions” of the type 


(qi ~ tq7) 
must be introduced.t 


+ We take, for example, an average of the form 
(pa (ri, 4) Ys (re, t2) Pola, ta) Polls, %))- 
In this case we shall have s = 4 in (2.21), and 
Dow; woes PQO (as At) Opal te) Apyo(ts) Ap, (tar 
~ rots) Aprolts) Apso(ta) Apro(ta))Fo- 
From the selection rules, Dp,o;,..; ps can be non-zero here only if 


Pi=Pa» P2=Ps 
or 


Pi=P3, pPe=Ps 
so that / = 2. As stipulated above, we denote: 


a= Hh, Po = a 


4 1 
() ER | ee en 
a Pid; vol TT (p?+M,) 


We have 


2 1 
= YI Deo: a0: 210; 201 Be G@i+™,)? 


1s Oa 
3 1 
Deoieareo: Tl ——_,, 
Ee &, | €193 @205 @29;5, ool (q?+ M,)* 
(01 Fee) 


since in the second sum on the right-hand side we must exclude the term q, = q2, 
which is already accounted for in the first sum. 
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Since all the terms in the sums (2.22) being considered are positive, 
removing these exclusions only increases the upper bound. Thus, 
using (1.88), we find from (2.21) that 


3/2 
|Dh)| < 9X, ls) lam ; (2.23) 


But, on the other hand, it is not difficult to see that 


1 1 dq 


Ve @+Mie = Gay | Grm,F tT 


4n [ R2dR 


— Gm | Germ thy 


where 
ty +0 as V + c, 
Therefore we can write 


ve Ge mae = 


where C(M) is a constant depending only on M. We have now from 
(2.23) 
| Dh)| = 6PK,l(s)C7(M,). 


Thus, we have proved that the inequality: 


JP Arr, 5 td {polis tr) «+ Pos(tss ts) r 
= (Polis ti) +++ Pofltss ts)rg} ar... drs| = |Dy(h)| 


< C(s,XKn, Mn) {(2a scalp |vev+e5 vi | ) i (2.24) 





x Yilu-t PE re: Ml 


is true for any function A(r,, ..., r,) of the class 2. Here, 
C(s, Kn, Mn) = 2+ VIS) KyCo(My) 
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is a constant, depending only on s, X%, and M,,. The relations (2.24) 
also demonstrate the validity of the asymptotic relation (2.17). 

We shall make one more remark with regard to the possibility of 
obtaining upper bounds which are uniform in the temperature interval 
(0 < 6 = 6) even when the condition E(f) = y is not fulfilled. 


§ 5. REMARK ON THE CONSTRUCTION OF UNIFORM 
BOUNDS 


As was remarked earlier, for X we can always take the expression 
(1.56): 


| ee 
Ks 6 a/ My. 
In the case when, in the whole space ¢, 


E(f) = y, (2.25) 


where y is a positive non-zero constant, we can put (cf. (1.57)) 


x=2, 
Y 


where 
3 2\2 12 

in the inequality (2.24). Thus, in this case, if ey + 0 (V > -) uni- 
formly in the temperature interval (0 < 6 « 6), our bound (2.24) 
will also be uniform in this interval. We shall show that this is so even 
if the condition (2.25) is not fulfilled. 

In this case, we shall slightly modify the inequality (1.88), which 
we used to give explicit form to the relations (2.21). We fix some 
y > 0. Then (1.88) is true for X = Q/y in the region in which 


IT(p1)| > y, ---, |T(ps)| =v, 


since it follows from the inequality |T(p)| = y that E(/) = y. 
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In the region where at least one of the |T(p,)| satisfies the inequality 
|T(p;)| < y, we make use of the trivial bound 


|(Ay,(t1) ina Agits))r— (Ag,(t1) area A, (ts))r,| <2. 


Thus, we can always write 
| Doyo: seed 5 P40, =< ED+ yaa g 2 Wert] + ie 7 6,{T (p,)} 


+2 py {1-8,(T(p;))}, 
where 


1/2 s—1 


M.M. : : 
gs) = zor] Mabe ay > +m(ver+o! vm] pe |tj—tysals 
fF 


G(x)=1, if |x| =y, 
6(x)=0, if |x| <>». 


Hence it follows that 
: Q,/. M 
[Dprats ceotpca| al ase oe cyt ay 
s 
+25 {1-8,(T(p,))}. 
Starting from inequality (2.21) and repeating, with this modified 


bound, all the previous arguments, we obtain the following inequality 
in place of inequality (2.24): 
1 j-—7\"2 
vi | ) 


>< , Q M2 
_ (s+1)/2 = 
x DF le gaal + etme & [ert alsa, (2.26) 


MiM2 








Dui] =c1s, X, ma) {( MMe ae] Vers 5D 


where 


A= Kultsya Ave, as De tore (T(q))) } il am} 
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But, obviously, 


1 s/2 1 
psie MW 2 Gale {! 7 5(T(4))} U (G+™M,P 


Sore 


s/2 1 
ys MW 2 Gale {! 7 BTa@)}TT (qé+ M,)? 


1 1 1 1 (s/2)—-1 
7 7 yl -U@) arm ly > GF 


1 


1 
<{(CM errr -F Gime: 


q 
(|(@%/8m)— «| <r) 


Further, we have 


q (q?+M,)? (2x8 +My vr 
(@%/am)— | <7) (| @*/am)—4| <7) 


where C is a constant, and 


1 _ 44 _F dq 

(27) Q?+Mi (2x) (¢?+M,)? 
({(4#/an)—x| <7) (\@*/2m)—p| <y) 

gle. _laldg amy 

(2x)? (q?+M,)? = (20)?>M}? * 


(\q2/2m)—4| <y) 


Thus, from (2.26) we obtain 


a = —— I 
|DHI|-<CUs, Ry Mn) | (MG + Ma] Ver FVM)" 








s—1 a W ae 
OF i ly— tas [+ aerne 2 r+ te} (2.27) 
j=1 y 4V 
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The quantity » can be chosen arbitrarily. Here, we shall put 


Me \14 
y= ker a) ’ 


where k is some constant. Then we shall find from (2.27) that if «, + 0 
(V + ~) uniformly in the temperature interval (0 < @ < 4p), then 
the right-hand side of the inequality (2.27) will also tend to zero 
uniformly in this interval as V + oo. In the same way, analogous 
observations can also be made in later arguments of this chapter. 


§ 6. GENERALIZED ASYMPTOTIC RELATIONS FOR THE 
GREEN’S FUNCTIONS 


We turn now to the Green’s functions. First of all, we shall consider 
a so-called T-product of our field functions, which we shall define here 
as follows: 


T{o,(t1, t1) «++ Pots, ts) 
=} HP O(t1— te)... Ots_1—ts)Po,(F1, ti) --- Po, (ts, ts)} (2.28) 
where 


1, t>0 
A(t) = | 0 ree (2.29) 


@ denotes a permutation of the field functions in the field-function 
product, accompanied by the corresponding permutation of t1, te, 
...,¢, in the product of 6-functions; e.g. 


PiG(t1—te) ... Ots_1— ts) Po, (ra, 1) -- - Pos, ts)} 
= 1p Wt), — tin) + Gi bye) Poin» th) «+» Poy (tin t)-" 


N@ =+1, according to whether the permutation @ is even or odd. 
The summation runs over all permutations D. 


t Here the indices j;, fz, ..., j, are obtained from the set of indices 1, 2, ..., 5 
by means of the permutation D. 
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Arranging (1, te, ..., t, in decreasing order from left to right 


hy 


eh, > ey, 


we have 
T{Po,(r1, hh)... Pat's, t,)} = £ po, (tj,, t,,) see Pay (h),» t;,). (2.30) 


It is clear from (2.28) and (2.29) that the T-product introduced here 
will vanish if any two time arguments ¢,, t; coincide. Since we can 
always substitute the product (2.30) into inequality (2.24), we see 
that 


fan, eas Rs) {(T{po,(r1, th)... Pols, ts))r 
—(T{po,(r1, t) sae Po (t's, ts)})r,} dr; aes drs 


<C(s, Kn, Mn) ( ts ranfvers Gh vin} (2.31) 





4v 
s—1 Me 
lf $+1)/2 pesbeeads 
x wilt tyyi| +26 XY er W ), 
and, consequently, for the many-time Green’s functions 


(T{...}), 


we have generalized asymptotic relations of the form 


(T{Go,(t1, t1) «. + Po,(s, ts)})r 


(2.32) 

—(T{o,(r1, hh) Past Po,(Fs, ts)}r, +0 

as V+, 

Finally, we take the two-time Green’s functions 

Gr rt, T |i, 6.5 Pes Peta, oo. Ps) (2.33) 
= 0t—1)([A(), BO) 4x, 7. 
where 

A(t) = Po,(r1, ft)... Po, (Pk h), (2.34) 


Br) = Paya s(Pk+1> 3 eee Pots; T). 
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It is clear that we can substitute the field-function products 
A) Bt), Blt) A) 


into the relations (2.24) and multiply the result by 6(¢—1). We thus 
obtain the inequality: 


| f alr, Gass rs) {Gr—Gr,} ar, nsrerce dr, | 


= C16, Xn May ( My tan| vere Vi] s|r—e1 00-0 





+ 20(t—1)20+ DAK Vv Mal (2.35) 


whence follows the generalized asymptotic relation: 


Gr—Gr, - 0 (2.36) 
as V+ o, 
To summarize, we can see that the following theorem is valid: 


THEOREM III. Jf the conditions I are fulfilled, then on the class £ of 
functions h(r,, ro, ..., t,) we have the following generalized asymptotic 
relations for the correlation functions and Green’s functions: 


(Po, (" 1, ty)... Pa,("s, ts))r — (Pao, (hi, ty)... Po,(¥'ss ts))r, +0 


as V+, 
(T{Po,(¥1, £1)... Pots, ts)})r —(T{Po,(t 1, t1) «+ Pots, ts)}>r, > O 
as Vo, 


and also 


O(t—T) ([Po,(r1, f) --- Poker Ds Pra i(e+1 T) «++ Pols, T))4)r 
—O(t—1) ([po,(t1, £) --- Pox (ks 1), Pou ss(Ye+15 T) +++ Pols, T)+ra-> O 
as V+, 


It clearly follows from this that to prove the existence of the gene- 


ralized limits 
lim (r. é ) 


V —> 00 
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of the averages considered here, it is sufficient to prove the existence 
of the generalized limits 
lim ts ‘ rr, 
(V > o) 
of the corresponding averages for the trial Hamiltonian I°,. We shall 
now investigate a proof of the existence of these limits. 


§ 7. THE EXISTENCE OF GENERALIZED LIMITS 


Apart from the conditions I, we shall now have to impose the 
following additional conditions, which we shall call conditions I’: 


1. The functions 
Af) = A(p,o) (6 =+4) 


are defined and bounded in the whole space E of points p and 
are independent of V. 

2. The discontinuities of these functions form a set of measure 
zero in the space E.* 


These conditions I’ mean that as a result of the passage to the limit 
(V + co), we go over from sums to Riemann integrals. In fact, the 
integrability in the Riemann sense of some bounded function is ensured 
by the fact that the set of its discontinuities is of measure zero. 

We shall take some function A(r, ..., re) from the class 2 and 
write, using (2.13) 


[rs Vo, ...,0s) (Po,(h1, th) aes Po,(ts> ts))r, dr, ... dt, 


1 
= per Y h(er Pr, ++ Ess) (Ap,o, (41) ade Apo {t))r, - 
Pir +++» Dos 


(2.37) 


By virtue of the rules of Bloch and de Dominicis, we have 


(Ap,o,(t1) ease Apo Ats)ra = Y n II (Anjo Ati) Apofti)ra> (2.38) 
t See Shilov and Gurevich (1964). 
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where 7 = +1. The sum here extends over all possible ways of pairing 
and, therefore, the number of terms in it is bounded by the finite 
number /(s). 

Further, [] is a product of s/2 factors. For a given pairing, p, = p, 
if o,=0, and p;=—p, if o,=—o, We denote the set of left 
indices p; in a given product by q,, ..., Qs respectively. Then, for 
a given way of pairing, each of the indices p,, ..., p, will be equal to 
q or —q. For a given way of pairing, these signs are completely deter- 
mined by the signs of 0,, ..., o,. Therefore, for such a pairing, we can 
write: 

mh(erPi, .. +» &sPs) = H(qi, -- +5 Qsi2)s 


where H(q,, ..-, 952) is a continuous function in the space E“”), 
satisfying the inequality 
s/2 1 


A(qi, -.-, Qsj2)| = aT RESP oRU EC 
|H(q qsi2) | Xu] G+M,? 


On the other hand, because of the relations (2.7), which give expres- 
sions for binary correlation averages of the form 


(Apyo,(t)) Apo t, Tas 
the function 
Odi. «+5 dara) = TL (Ayo ts) Areal ti) drys 


will be a function which is defined in E“” and independent of V. 
We have now from (2.37) 


[i 1 ey Hs) (Po,(F1, C1)... Po,(ts, ts))r, at... ats 
(2.39) 
= 2 [jaw 2, Hae ad Or. sau, 


veeee Aele 


where in the first sum } the number of terms is bounded by the 
number /(s), while the second sum 


(...) 


Ms -+ +9 Gai2 


is taken over the quasi-discrete set E(7”). 
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We note that, since the set of discontinuities of the function A(p, ¢) 
has zero measure in E£, the functions (2.7) will also possess this pro- 
perty. Thus, the discontinuities of the functions 


Aq, sees s/2) a1, we Ry qs/2) 


form a set of measure zero in the space E®”, 
On the other hand, since 


|O(q1, sey qsi2) | <= 1, 
we have the bound 


8/2 1 
Hi yg eeeg Qe eeg s v4 Ee Tare, Yt 
|H(qi 9s.) O41 qs/2) | mul G+M,? 
In view of the above, 


1 
pre Ways dO «a 
Gar oss Tale 





1 
a oar | vee [aa + +05 Wiz) O91, «+s Qsi2) 4g... Asia 


as Vso, 


and therefore, using (2.39), we obtain 


[te aeey Ps) (Po, (hi, th) sito Po,(ts5 ts) )r, dr; Ss ave dr, 
1 (2.40) 
pa Laan | Ha. ++ +5 G52) (G1, -» +5 Gsiz) dq ..- dGsia 


as V+o, 


where in the sum )° there are no more than /(s) terms (As) is the 
number of different ways of pairing). We emphasize here that the 
order of the field operators in the averages (...)r, can be chosen 
arbitrarily. 

In addition, the relation (2.40) also holds for any real parameters 
t,, .--, f,, so that, for example, this relation can be multiplied by 
products of functions of the form 6(t,— t,). Hence follows the existence, 
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on the class 2 of functions h, of the generalized limits 
Rll (Po,(r1, t) ae Pots, ts)r,> 
Rs (T{o,(r1, th) enced Po,(rs, ts)})r,, 


ao O(t—t) ([Go,(r1, th)... Poy("k: nh), Pozss(Pk+1> t).. Pots, Tar: 


Thus, we can see that the following theorem is valid: 


THEOREM IV. If the conditions I (Chapter 1, § 1) and I’ (Chapter 
2, § 7) are fulfilled, then, on the class £ of functions h(r,, ..., 1,), the 
following generalized limits for the correlation functions and Green’s 


functions exist: 
im (Po,(r1, th) ce Pots, ts))r 


= yim (Pot 1, t1) «.- Pols, ts)re3 


a (T{Go,(r1, t1) .. + Pots, ts)}r 
= um (T{Go,(r1, th) ves Pots, ts)})r,3 


line O(t—t) ([po,(ri, t) eh Po, ("ks n), Poy (rks 1, t) Bare Pos, t)).>r 


lim O(t—1)([Po,(r1, f) «Polk 4), Pop ys (Pk +197). » Po,(Pes T)] ar, « 


As we have seen, expressions of the form 


lim (...)p, 


V—> co 


can be calculated by applying the rules of Bloch and de Dominicis 
and replacing the quasi-discrete sums by the corresponding integrals. 


CHAPTER 3 


CORRELATION FUNCTIONS FOR 
SYSTEMS WITH FOUR-FERMION 
NEGATIVE INTERACTION 


§ 1. CALCULATION OF THE FREE ENERGY FOR MODEL 
SYSTEMS WITH ATTRACTION 


In this chapter we shall study dynamical systems which correspond 
to attraction of fermions. We shall begin by calculating the free energy 
for model Hamiltonians with four-fermion interaction. This problem, 
as we have shown (Bogolyubov, Jr., 1966b), is of great interest in the 
study of model problems in the theory of superconductivity and serves 
as an example of an exact calculation of the free energy for model 
systems of the BCS type (Wentzel, 1960, and Thirring, 1968).t The 
results and upper bounds obtained here also constitute a proof of the 
Theorem 3.1 formulated in this section. 

We shall start from the Hamiltonian 

H=T-W YY J,Jt. (3.1) 
1lsa<s 
If we take the following Fermi-operator expressions for the operators 
T and J,: 


T=P Tapa, Jam wv Adf) atat,, (3.2) 


t This treatment also solves a number of problems raised by Wentzel (1960) 
concerning the asymptotically exact calculation of the free energy. and can also 
be applied to certain quasi-spin Hamiltonians (cf. Bogolyubov, Jr., and 
Shumovskiy (1970, 1971) and Shumovskil (1971). 
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we obtain the usual BCS Hamiltonian 
1 , + 
= + a,p— —— fe Af aL faA_pAy’. 
H= PTS) aa oy » ALL) af a2 4_pay (3.3) 


In fact it is not necessary for our discussion that the operators T and 
J, have the explicit form (3.2). 
It is sufficient to impose the following general conditions: 


| Ja| <« M,, {TJ.—J,T | = Mo, 

\JtIp—Ipdt a2, tly Iola =, _ 
where M,, Mz and M3 are constants as V - co and the symbol |.. .| 
denotes the norm of the indicated operators. We assume also that 
the free energy per unit volume for the Hamiltonian H = Tis bounded 
by aconstant and that the number of terms s in the sum (3.1) is fixed. 

We thus start from the Hamiltonian (3.1) with the conditions (3.4). 
We take the trial Hamiltonian to have the usual form 
H°=T-W Y (CJt+CV,)+W Y ICI. 
l<a=<s l<a=s (3.5) 
Here the C, are complex constants determined from the condition 
that the function 


fu(C) =— , 6 In Tr e-#%6 (3.6) 


have its absolute minimum value in the domain of all the complex 
variables (Ci, ..., C,). We shall denote this complex set of points 
(Ci, ..., C,) by {£*}. By making use of the minimizing values of C, 
we calculate the free energy per unit volume for the trial Hamiltonian: 


Su = min fy0(C). (3.7) 
{E*} 


We also take the corresponding free energy for the Hamiltonian 
(3.1): 


tu = -} 61n Tr e~#, (3.8) 
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We shall prove that the difference fj.—f, tends to zero as V + ~, 
For this it is more convenient to consider first the auxiliary problem 
with the Hamiltonian 

T=H-V ¥ (vJ.+v2Jt), (3.9) 
1<exs 
where v1, ..., ¥, are arbitrary non-zero complex parameters. In this 
problem, the corresponding trial Hamiltonian has the form 
P=H°-V Y (v,J,+v2Je). (3.10) 


l<a<s 
The complex quantities C = (Ci, ..., C,) occurring here are also 
determined from the condition for the absolute minimum of the func- 
tion 


fr(C) = -< In Tr e-7%9 (3.11) 


We shall obtain an upper bound for the difference fpo—fy, and 
show this difference to be asymptotically small as V + oo. Here, 


fro = min fro(C) 
{E*} 


and f, is the free energy per unit volume for the Hamiltonian I’. 
Although this bound for f-o—fp will be found for |v,| > 0, it turns 
out that it is uniform with respect to », + 0, so that we can then 
pass to the limit v, =O(1 «a <5). We thereby obtain a bound 
for fyo—fy, proving it to be asymptotically small as V+ -. 

We therefore begin by treating the trial Hamiltonian J®. It is not 
difficult to show that the problem of the absolute minimum of the 
function (3.11) has a solution and that this absolute minimum is 
realized for finite values C, = C2 (1 <k <s). This can be seen by 
using the inequalities 


2 (1C.)?+(|Ca|+2My)}—4Mis+y+fr = fro(C) 
= » {IC.? + (|C.| —2My)?"}—4Mis—y+fr 


=P IC,?-4Mis—ytfr, 
C=(Ci,...,Cs), p=2Mid|n| (l<ass). 


(3.12) 
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Thus, the function f7(C) has an absolute minimum at some point 
C= (C2,...,€%. 

Because the function f7(C) is continuously differentiable, at the 
point C = C° we have 





TAO) =0 (l<ass), 
i.e. an equation for C,: 
Tr Jne— 78 
C. = (Jaro = “Tr e-te . 


Taking into account the conditions (3.4) on the operators J,, we 
obtain |C, | < M, = const. 
We turn now to the derivation of inequalities limiting the difference 
in the free energies per unit volume fp—/fy in terms of averages of 
SF 2-3 F U,-ej02-&, (3.13) 
V l<e<s 
For this we note that 1° = J+ and introduce the intermediate 
Hamiltonian 
It = 1%+ tot, 
which for t = 0 coincides with the trial Hamiltonian (3.10) and for 
t = 1 coincides with the original J‘ (3.9). The constants C = (Ci, ..., 
C,) occurring in I“ are assumed to be fixed and to be independent of 
the parameter ¢. 


Let us consider the partition functiont and free energy for the 
intermediate Hamiltonian I“: 


Qr=TreF;  f(Cr, 5.) =—F In Qs Qe = eM 


(3.14) 


Differentiating the equality (3.14) twice with respect to ¢ using operator 
differentiation, we arrive at the formula 


1 
mae v Ot aa Of - = 1 —(P/8)r —(F/0K1 — 7) 
o oe e (a) => vg, {7 {he ote } dz. 
0 


+ See footnote added in proof on p. 136. 
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Taking into account that 


Of, _ 1 Trote-r 


1 
eV Trent = 


we find 





- os AN - C10)" fe-TVN1-") dr — (oA)? 
= 57 r 0, fr ote ote 1-1) dr—(A) 
0 


— (VO Be— (178) — 7) 
— 0, | {Be ‘Be } dr; 


B= A-(A). 


Going over to a matrix representation in which the Hamiltonian J‘ 
is diagonal, we have 


Of 1 —{((El,-E!)|0}<—Et/0 
— St = Wo, je » Brune 


a 2 p~{(E,—-E4)/O}x-EY/O 
_ st | #3 om e > 0. 


Hence it follows, in particular, that 6%f,/6t2<0, and therefore 
6f,/Ot = (A),/V decreases with increase of the parameter ¢. Further- 
more, taking into account that f, does not depend on C, we have 


fil) fo =— | Ff dt =— |S a= 
0 0 


Since this relation is true for all C = (Ci, ..., C,), we also have 


minfpo(C)> fr, ie. fro=/Sfp. 
{E+} 
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We shall integrate the inequality (cf), = (ct) (0 < t = 1). Substi- 
tuting the expression (3.13) for of, we convince ourselves that the 
inequality 


Sr(C)—fr = 2 Ps (Ja—Cz) JE —CD))r 


is true for any C = (Ci, ...,C,). We put here C, = (J,)p (I <a <5) 
and note that 

fro= mn SroC) = fro((J)r). 
Thus, 


fr—fr=flr)—-fr =? (Je—VJa)r) JE —Jn))r 
and, finally, 
O<fre~fr=2 yr (Je—(Je)r) Jt —WJt)r))r, — (3.15)t 


where, as always, 
Sro = min fro(C). (3.16) 
{E*} 


Let us recall our main problem. We want to show that the difference 
Sr—fr is asymptotically small as V + o. It follows from (3.15) that 
we shall have solved our problem if we can demonstrate the asymptotic 
smallness of the average on the right-hand side of (3.15). 

Taking into account the main idea of a paper by the author (Bogo- 
lyubov, Jr., 1966a), we express this right-hand side in terms of 
&f/dv! dv,. Differentiating, we have 


I or Vv 1 Tr (D@e-C/OrD+ @e—(A—-V9}T) dr 


~6 avian, & Tr e-F/8 ’ 
0 


where 
D© = J,—(J) (l sas). 


t For convenience, in the remainder of this section we shall omit the subscript 
in statistical averages over the Hamiltonian I, i.e. we shall write (...) = ¢...)r. 
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Going over to the matrix representation in which I’; is diagonal, 
we find 


1 
1 Of - V (2) —(/0)Em py +(@),,— (1-1/8) En 
a 6 av Ov" Ov, ~ O06 X | vite Dinan e dt 


1 


V a -~(a- 
= ow Y —s (#/0)Em — {(1—1)/8} En, dt 
nm 
0 
—Enld _ o—Eal8 


Gis 


Using Hélder’s inequality, we have the following bound: 


oe | D® |? |e Em? _ e~Ew/o| 


oer \*5 i 2 : 1/3 
«(~32%,) (5: % | Diim|* |En—Em|? (e EmlO 1 @ =) 


We carry out the simple transformations: 
V é “ ee 
5 ©, |DEAI |En— Em Me" Fo” + 0° Fo) 
= S Tr e-79 ((TD©— D©L) (D+@Pr—TD+) + (D+@P—TD+) 


x ([D®— Der) 
= VW(CI.—IeP) (Ie —IeV* + (TFe—TaP)* Tn— TeV) 
=2VM, 


where 
SM = M2+4MiM3s+2M3 Y |». 


lsaxs 


Hence we obtain 


V 7 2 = or 2/3 
5 Smeal lente |< (ara) anys. 
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Furthermore, 


V a)|2— V Dam 2 i 
OD, Dime wate Ds eel ch eB Enl@ 6 w/0) 


+g Epil lee eh, 
a,m 


where 
V @iz Rio V 2)» +(a), — 1/0 2) py +(a) 
@ &,| Dim? e =o HDD e 7 — yD Dt) 


= V(Ja—(Je)) (JE —(J2))). 


Thus, we finally obtain 


((Ja- WJ) JE - (JP) 
arf 2 ony ( Of ie 


~ Out Ov, | V V23 \ Ove Ov, 


Substituting this inequality into (3.15), we find 


0= fro—fr=2 a n> (- arm) 


<acs 


(3.17) 
af 28 
a Sas aes j are) , 
Hence we can see that our problem would be solved if we could show 
that the second derivatives |62//Ov! av,| are bounded by a constant 
as V + co. Unfortunately, we are unable to prove such a statement. 
We must start from the boundedness of the first derivatives | 0f/0», | 
<M, (1 sa <5). 
Because of this, we develop a method in which it will not be neces- 
sary to use the boundedness of the second derivatives and by means 
of which we can demonstrate the asymptotic smallness of the difference 


a = fro—fr. 
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For the following, it will be more convenient to transform to the 
polar variables r,., 9.: 


ry = r, Va, va) Pa = PalPas v3) (1 <a 5) 
in inequality (3.17). Accordingly, 
f(@, vt, Ce | Vs, v3) >f(n, Pi; eee ls, Ps). 


ef -a{t a Z\+3 of al (3.18) 


Then 





Ov Ov, 3894) r. Or,\* Or.) AP r 
We now make use of the inequality 
[a(ra, .. 5 %s3 1, «+ +s Ps)|  la(h1, «- 5 785 P2, +++ Ps) 
—a(Er, ...,&s5 M1, ---,s)| +]a(Gr, ....&s3 My +++. Ns) 
2) eer 


Ore lata l<a<s Ov. 
+ fa(é, sey E53 Nhy vos Ns)| E) 
and take 


= 








|Pa— Ma | (3.19) 


max 











l<a<s 


ntl<t<ntl, ge<te<petdy b=, (lea<s), 
so that 
a(E1, ..-5 Es3 M1, - + +> Ms) 


n+l r,tal +4, +8 


j dr, ... dr, f ie (dps... dpsatray «5153 Pty «6 s) TI] 7 
l<a<cs 


— matt retl 71 ee 


( oie TT (e+ 21)P?—(@a +1)" 

















(3.20) 
We note that 
| o <=2Mi, <= 2Mir, 
Ore 
(3.21) 
Ga 
and | ar. ae ~ | aMire . 
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Therefore, the first two terms in inequality (3.19) can be bounded as 


follows: 
da 
or, 





1G. | 


max 


|r. —€| + Dy 


Oa 
1<a<s Opa 


I<a<s 


<4M,s-21+4M,ls = 12M,sl. 


max 














(3.22) 


Starting from the formulae (3.17) and (3.20), we find a bound for 
the expression a(£1,..., &3 71, ..., 7%). We put (3.18) into the right- 
hand side of the inequality (3.17). We then multiply (3.17) by the 
product rire ... r, and integrate it over all values of the variables 
1}, «+5 ly Pi, -+-, Ps Within the following limits: 


tet] < ry < y+ 2l, Pa = Po = Pat 42, é= 4 (<e<s) 


a 


We then obtain 


Of n- Patrs eet oa os goiter da «odd 1. ADs 


i] 
< ay | o> | tres wee Pst Forts... Pst... 


+... +¥Fyrire ... tsi} dry ... drs dpi ... dps (3.23) 


S23 
+ Tas [- : [ e vee Heb FRB AB ret. 


+... +F28ryry ... 13} dr, ... dr, dp, ... dp. 








R=a(" “Go te on (G2)-0 (i =a <5). 


By considering the separate terms of the part of (3.23) containing 
the factor 0/2V, we see that bounds can be found for these by integrat- 
ing successively over r, and ¢, in each of them (1 < « <5) and using 
the inequalities (3.21) to find bounds for the resulting first derivatives 
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Of/Op, and Of/Or,; for all the terms of this sum, we then find 


() 2M1(6p+2) 2 
WF Je TET LL, (Oot P(t De 8.24) 


By applying Holder’s inequality and using analogous reasoning for all 
the terms in the sum containing the factor ?/V?, we obtain 


SN?! ((6-+2)2M1)?8 : 
yrazis es ~ Tp GES J, {(ra +21)? (ra + 1)?} 62. (3.25) 


Using the formulae (3.19}(3.25), we now obtain a bound for 
Sro—Sr = a(ri, sea Mss Ply ve ey Ps): 
20M, (6+ 2) 


V ee 155 


SH (6g+ 2)? 
+ pax OM FR 


0 <fro—fr < 12Mils+ 








(3.26) 


On the other hand, we note that 6, = //r, and choose R in such 
a way, 
R= |r, sey |sl, 


that R= 7, (8 = 1,...,5)and 6, = 6 = 1/R. Then, using the obvious 
inequality (J+ 2R)*? < [3 4-(2R)**, we find from (3.26) 








20M: 4M;)23 
US aee Vie ADR) a 3A Dae 
(4M,)?3 (3.27) 
pis 5 pas 22/3 R23 | 


We now choose /, which is an arbitrary positive quantity, such that 


SN23(4M 8 
12Mil = rss 
Then 
2/5 
l= at P AM = const. 


ys” ~ "22/833/5 gu 
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Putting this expression for / into the inequality (3.27), we find 


P  20M;s  40MisR 
0 <fro—fr = 24Mis ag + aap “Pispr 
SMNs 


+ Fspas AMA YR2PR® for |m|<R (1<a<s). 


Hence it is clear that the difference f,.—f; vanishes as V > ~. 

We note that in the above bound, we can take the limit », = 0 
(1 <a <s) and finally prove the statement we made earlier about the 
asymptotic smallness of the difference fpo—fy: 


Safwan ees 


pat yap’ 3:28) 


where P is a simple combination of the original constants M1, M2 and 
Ms3. It is also clear that the above bound is uniform as 6 — 0, and, 
therefore, the inequality (3.28) is valid for 6 = 0. 

We have thus proved the following theorem: 


THEOREM 3.1. Let the Hamiltonian of the system be 


H=T-W Y. g,J.Ji (3.29) 


l<axs 
and let the operators T and J, in (3.29) satisfy the following conditions: 
T=T*, \J.| < Mi, 


[TI-J.T|<Ma, \JeJp—Iphe< 2, 3.30) 
|JtIp—JpJz| <™, My, M2, Mg = const. 


In addition, let the free energy per unit volume, calculated for the 
Hamiltonian T, be bounded by a constant: 


|f(T)| <= My = const, (3.31) 
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We construct the trial Hamiltonian 


H(C) = T-W ¥g(C. JE +Ci J.-C) (3.32) 


where C = (Ci, ..., C,), and Ci, ...,C, are complex numbers. Then 
the following inequalities are valid: 


; 1 
O< min A{H(C)}-f(A) < (7) : (3.33)t 


where &(1/V) + 0 (as V + ©) uniformly with respect to 9 in the 
interval (0 < 6 =< 6) where @o is an arbitrary fixed temperature. 


§ 2. FURTHER PROPERTIES OF THE EXPRESSIONS FOR THE 
FREE ENERGY 


Having formulated Theorem 3.1, we shall now study the question 
of the existence of the limit 


lim f(A). 
V-> 00 


We shall assume that, in addition to the conditions of Theorem 3.I, 
the following condition is fulfilled. For any complex Ci, ..., C,, the 
limit 


lim ftH(C)} 
V-—>oco 
exists. We shall denote this limit by 


Soo {H(C)}. 


t (1) We shall denote the free energy per unit volume for some Hamiltonian 
A by f(A), or, if we wish to emphasize its dependence on the volume, by /;(A). 
(2) By min f(C) we shall always mean the absolute minimum of the function f(C) 

(©) 


in the eee of all the points C. 
(3) e(1/V) is given by (cf. (3.28)) 


| 1) _ 24MsP , 26Ms 
) Veves — VeveP ° 


where P is a constant and g > 0 is the smallest of the gz. 
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We put! 
FC) = f{H(C)}—2Y &C.Cz 


and note that 
OFC) _ 








ac —2g.Ji)Ho> 
OF YW(C 
) = —2g.(Ja)HO : 


Therefore, from (3.30) we have 








OF (C OFC 
ORO. |< 28M | ne <2gM;, 
whence 
|FY(C’)—FAC")| <4M, ; YX ga|Cr—-Cy' |. (3.34) 
saxs 


Thus, the set of functions 
{Fy(C)} V+) 
s uniformly continuous. 
Since we have the convergence 


FC) ~ F.(C) = foolH(C)}— 2) 8sC.Cz (V + ~) 


at each point C, we see that this convergence will be uniform on the 
set /M(R) of points C defined by the inequalities: 


IC1| =R,..., ICs| <=R, 


for any fixed value of R. 
Therefore, 


|fAH(C)} fo (H(C)}| = [FAC)— F.C) (3.35) 
=nv(R)+0 V+) 


t We stress that the function F,(C) need not be treated from the standpoint 
of the theory of functions of a complex variable. It is not difficult to see that, in 
essence, the function Fy(C) is a function of real variables and that these can be 
taken as the real and imaginary parts of the variables C,. 
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for 


CéeM(R). 
On the other hand, it follows from (3.34) that 
|F(Ci, Ce, ..., Cs)—Fy(O, C2, ...,Cs)| = 4Migi |Ci}. 
Hence, we have 
S{H(C1, C2, ..., Cs)}—f{HO, Co, ..., Cs)} 


= Fy(Ci, Ce, ..., Cs)—Fy(0, Ce, ..., Cs) +291|Ci|? (3.36) 
=—4g1M, |Ci|+2¢1|Cil?. 


We denote the lower bound of f{H(C)} in the space of the points 
C by 
inf {(H(C)}. 
(Cc) 
Obviously, 
S (HO, Ce, sey Cs)} = ine) 


and, therefore, it follows from (3.36) that 


JC) MEO) = 291 |C1| (|Ci|—2M3). 


Replacing C, by C, (« = 1, 2, ..., s) in the above discussion, we find 
also 


S(H(C)— int F{H(O)} = 2ge|Cu\(|\Cx|-2M), % =1,2...5. 


Hence it is clear that if 
|C.| > 2M, 
for at least one «, then 


S{H(C)} > nt AG). 
Therefore, the lower bound of f{H(C)} on the set M(2M;) is equal 


to the lower bound of this function on the whole space of points C. 
Since f(H(C)} is continuous and the set (2M,) is bounded and closed, 
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this lower bound is attained on (2M), i.e. an absolute minimum 
of the function under consideration exists and is realized at certain 
points :* 

C = CME MM(2M}). 


On the other hand, taking (3.34) into account and passing to the 
limit V ~ o, we find 


| F.(C’)—Fo(C")|<4Mi1 Y) 8a|Ca—Cr’ |. 
l<a<s 


Hence, repeating exactly the above treatment, we see that the func- 
tion 
fu{H(C)} = Fu(C)+2¥ geCeC? 


also has an absolute minimum in the space of all the points C, which 
is realized at certain points 
C = CeM(2M)). 
From (3.35), we have now: 
Sof H(C)}—f{H(C)} = nv (2M1), 
{H(C)} —f-{H(C)} <= (2M). 
But, by definition of the absolute minimum, 


Fo.{H(C)} = f-{H(O)}, 


S{H(C)} = AH(C)}. 
Consequently, 
Sool H(C)} —f{H(CM)} =< nv (2M1), 
S{H(C)} —f.{H(C)} < nv(2M1), 
or 


| f-{(C)}— min F{H(C©)}| < 4y, | 


t Generally speaking, the point C” of the absolute minimum is not unique. 
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where 
by = 72M) +0 (V+). 


Taking Theorem 3.I into account, we finally obtain 
1 
—by < f.(H(©)}—f(H) < (7) bp: 
Thus, we have now proved the following theorem: 


THEOREM 3.II. If the conditions of Theorem 3.1 are fulfilled, and if for 
any complex values of Ci, ..., C, the limit 


fo(H(©)} = lim FHC) 
exists, then: 


1. This limit function has an absolute minimum in the space of all 
the points C, which is realized at certain points 


C = Ce M(2M)). 
2. The inequalities 


~ by = fal) FD) = e( 7) +97 


are valid, where 


1 
(7) +o by >0 as Voom, 


|AH(C)}—f-{HO)}1. 


by = max 
(C € M@M;)) 


§ 3. CONSTRUCTION OF ASYMPTOTIC RELATIONS FOR 
THE FREE ENERGY 


We shall now make a special study of those cases when the operators 
T and J, in the Hamiltonian (3.29) have the form (3.2). 
As can easily be shown, the conditions of Theorem 3.I will be 
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fulfilled in such cases if 
PUIMP A. <0, (2) 
PLM <O1, 0) (3.37) 
FLIM P=. ©) 


Here, « = 1,...,53 0 =+i, Qo, Q1, Qe = are constants. Then, for 
example, in the inequalities (3.30) we can put 


M1 =Q1, Ms3=Q2, M2 =2Qo. 


Here, let the functions 4,(p, o) satisfy, in addition to the inequalities 
(3.37), the following conditions: 


|A.(p,o)|<O, QQ =const. (3.38) 
The set of the discontinuities of the functions 1,(p,o) is a set of 
measure zero in the space E. (3.39) 


We shall show that in this situation the conditions of Theorem 
3.11 are also fulfilled 

Before proceeding to this problem, we note that the inequalities 
(3.37) and (3.38) are not independent. 

In fact, (3.37c) follows from the inequalities (3.37b) and (3.38). 
Also, (3.37b) follows from (3.37a) and (3.38). Thus, all the inequalities 
imposed here on the A, are fulfilled if the inequalities (3.37a) and 
(3.38) are true. 

We note further that (3.37a) and (3.38) hold if the A, satisfy the 
inequalities 

lA, 0) | gg (3.40) 
: (p?+a)>’ 
XX, a = const. 


We turn now to the question of the fulfilment of the conditions of 
Theorem 3.11. Since, in the situation being studied, the conditions of 
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Theorem 3.1 are fulfilled, we need only show that for any fixed 
complex quantities Ci, ..., C,, the limit 
Soo{H(C)} = aim SvV{H(C)} (3.41) 


exists. For this, we write the operator form of the trial Hamiltonian as 
follows: 


H(C) =. T(f)afa,—F x {A(f)apat, + A%f)a_pa;} 
Z (3.42) 
+2V ¥ gCiCu 


where 


A(f) = 2 » BCA Sf). (3.43) 


Going over to the Fermi amplitudes «,, «;, which are related to the 
old a,, a7 by the transformation (1.18), we diagonalize the form (3.42) 
and obtain 


H(C) = »» E( fap a, 
+ Ppy accra DAT}, 3.44) 


where 


Ef) = VT{/)+/AP)E. 


Hence, we find 
Sv{H(OC)} =- 4 In Tr e~ (Ov 


= 28, 8.0.02 w PEN-TH)} 


6 re 
eae —E(f)/@ 
) In(1+e ys 


or, separating the indices p and o in the summation, 
Sv {H(C)} = 2¥g.C.Cr 


by fne-(E-|) ow 


2m 


1 
OF 7D In (1te-Re omy, 
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On the other hand, it is not difficult to see that if we have some bound- 
ed function F(p), defined everywhere on the space E, whose dis- 
continuities form a set of measure zero, then 


1 1 
7 2a FO) > Gaye | FO) (3.46) 


S, 
for any sphere S, with arbitrary fixed radius r. In fact, such a function 
will be Riemann-integrable in the region S,. For the summation 
points 











Ga 27n2 r) 








CONES Ee Be Ek 
we have 
2n \3 2n 
ApxAp,Ap, = (Z) = oe > 
so that 
(2x)8 
F 
7 pie (P) 


will be the Riemann sum for the integral 5) F(p) dp. We note also 
that if ‘ 


Sy |F(p)| < A = const, 
V > 
then 


1 
T py |F) < 4 


Hence, passing to the limit V -> 0, we have 


ane { lF(p)| dp < A. 
S, 


Because of the arbitrariness of the radius r, we see that F(p) is an 
absolutely integrable function in the whole space EF, such that 


1 
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For a given function F(p), let the following inequality be valid: 


1 
yy Fi Ts 
Toe. fP Ml <n 


where (E-—S,) denotes the set of the points of E lying outside the 
sphere S,. 7, does not depend on V, and 


Un > 0, (r > co). 
Then, obviously 
1 1 
ae eee 3.47 
7 LF) One [ro dp. (3.47) 
In fact, we shall fix an arbitrarily small number ¢ > O and choose 
r = Fro such that 


Tre = 


Alm 


In view of (3.46), we can find a number Vo such that, for V = Vo, the 
inequality 
yy FO aus [ro dp| =< 
| V pe, (2x8 ; 2 


fo 





holds. We have, therefore 
1 1 
7 DF ays | Fp) do] 


1 1 
7 FO Gap [FO% 


fo 


= 








1 1 
Oy, Fe ag FOIt Cpe | |F(p)| dp 


—S,,+E 


éE E E 
<5+7+7=6 


2 4 4 
for any V = Vo, and this establishes the validity of (3.47). 
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After these trivial remarks, we turn to the expression (3.45). It is 
clear from (3.43) and the conditions imposed above on the A, that 


| A(p, )| = Ao = const, 
LIA, 0)? = As = const. 


We see also that the discontinuities of the function A(p, o), and conse- 
quently also of the function 


x,0)-(2—n) = (Z-2) +140, 9°-(2 a) 


form a set of measure zero in the space F. Further, for 


p? = 4mu, 
we have 
ee 
2m 4m?’ 
_(P_,\ 1A, _ 2m - 
0 = E(p, 0) (Sr #) < ae = > 4@, o)|?. 
4m 
Hence, 
1 Pp 2m 
T pedis, [8 (Sm -#) | = pe 
(p?=4mp) 


Thus, taking into account the remarks made above, we see that 
1 Pp ) 
__) LE(p, c)—|——— 
7 Re (Sa) 


Further, we have 


In {1 + e{El,, o)}/6¥ < e—E(p. 9/? < const e—77/28, 
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Since this function decreases sufficiently rapidly as p + ~ and its 
discontinuities form a set of measure zero, we also obtain 


1 —E(p, 0)/0 wih: —E(p, 0/8 
yuin{ite P, oy/6} (any In {1+ e~ HP. 9} dp. 


This also proves the validity of the property (3.41). We have here 
- ty] P 
foo {H(C)} = 2V KCCI 5D Gaye E(p, 0)—(5——n) ¢ dp 


= eS —E(p, 9)/6 
0S aie | in tlte P/9\ dp, 


or, more compactly, 


fo {H(C)} = 2¥ Bie sp | {E(f)—T(P)} of 
; (3.48) 
~ aap |! {1+ e- FU) af. 


Here, the integration 


f(...)df implies the operation J) f(...) dp. 


Thus, in the case under investigation, if the conditions (3.37a), (3.38) 
and (3.39) are fulfilled, the conditions of Theorem 3.II are satisfied. 
As was noted above, in the proof of this theorem, the convergence 


SAH(C)}-f.{HO)} +0, V+) (3.49) 


is uniform on any bounded set of points C. 


§ 4. ON THE UNIFORM CONVERGENCE WITH RESPECT TO 0 
OF THE FREE ENERGY FUNCTION AND ON BOUNDS FOR 
THE QUANTITIES 6, 


We shall show that, in the above case, the convergence (3.49) is also 
uniform with respect to 0 in the interval (0 < 6 = 69), where 6o is any 
fixed temperature. It can be seen that this property will be established 
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once we have shown that 


= (fr {H(C)}— foo (H(C)}) | =X=const (0<6<6y) (3.50) 


and this will be our aim now. We have 


E(f)e-2swe 
PAHO) =-7F hn {1 +e Fie} a ae 
Since 
E ene 2 
5 e < 2° 


we can write 


8 Pet asa 
Smo) | < vye seas 7a 


2 
< ( +3) re e— FSYI200 | 


Completely analogously, we find 


¥ e- Funi2e 
4 


(6) | 2 1 7 
aa Foo {H(C)} | < (: + 2) Gs One fe ENV/260 df 
and, therefore, 


SO ALO) 





2\} 1 1 
aes ye —E(f)/2H9 4. —E(f)/269 
<(1+2){ bre tape a}. 
In view of the rapid falling off of e— (1/260 as | p| + oo, the integral 


f e-FU/260 df 


has a finite value, and 


i -E(S 1280 a — E(f)/299 
y>° (ay | oem a. 
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Inequality (3.50) is thus established, and the uniformity of the conver- 
gence (3.49) with respect to @ in the interval (0 < 6 = p) is thereby 
also proved. Therefore, in the case under consideration, in Theorem 
3.II the relation 

by +0 V+) 


holds uniformly with respect to @ in this interval. 
Thus, putting 
1 


(7) +o = by ; 


we can formulate paragraph 2 of Theorem 3.II in the form 
| foo{H(C)} —f(H)| = 5, (3.51) 


6, +0(V +) uniformly with respect to 0 in the interval (0 < 0 =< 6p). 

An explicit expression for e(1/V) has been obtained. (See footnote * 
on page 101.) It would also not be difficult to obtain an explicit expres- 
sion for the bound 4, of the difference 


SAH(C)} fo {H(O)} (3.52) 


if we impose on the A,(f) the appropriate conditions of smoothness 
and falling off as |p| ~-o. In fact, as we have seen, (3.52) is the 
difference between a Riemann sum and the corresponding integral, 
so that here we can make use of well-known techniques from the theory 
of the approximate calculation of three-dimensional integrals. 

Thus, we can show, for example, that if the functions A, of the 
point (p) are continuous and differentiable everywhere, with the 
possible exception of certain sufficiently smooth discontinuity sur- 
faces, and go to zero, together with 0A,/Op, sufficiently rapidly as 
p>, then 
const _ const 

L ypis* 








by = 


In the case when the A, are everywhere continuous, possess deriva- 
tives of second order with respect to (p), and, as p +, go to zero 
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along with their derivatives of up to and including second order, we 
can obtain the stronger bound: 


Fs const 
v= Gals: 


§ 5. PROPERTIES OF PARTIAL DERIVATIVES OF THE FREE 
ENERGY FUNCTION. THEOREM 3.11 


We shall now study the partial derivatives of the function 
Soo {H(C)} = 23) BeCiCa 


(3.53) 
Sofas iees st -EV))\ df 
200) [ew T(f)+26 In (L+e )} df 
with respect to the variables C1, ..., C,, Cy], ..., C3. We have 
U = = {E(f)—T(f)-+26 In (1-+e-200)} 

a 2) OE(f) _ , tanh [E(f)/26] 

~~) Tp eke OC, COE) 


*) A* 4 
E(f) {Fee AN} a(S) 
But, by virtue of (3.38) and the inequality 


bie tanh x 1, 


we see that U is a bounded function of (p) in E: 


1 ; 
Ul =< Y 88 1Co| Dies. 


It is clear also that U is a continuous and differentiable function of C 
in the whole space of points (C). On the other hand, since |tanh x| <1, 
we also have 


IU = ep Y Bow lCol-l AP) af) | 
1 


P fee rene 


= 





aa sd 


FOUR-FERMION NEGATIVE INTERACTION 115 


Hence, it is not difficult to see that U(p) is absolutely integrable in £, 














and 
[ua [ud - | U dp|= | |U| dp 
: Ss, E-S, E-S, 

— 18. i 2 P 

: 2 3, 861C \Aa f)I? ap} |A(f) |? dp 
for 

r => 4mu. 

Consequently, 


fUdp ~ [Uap (r + «) 
S, 


uniformly with respect to C on any bounded set of points C. 
Thus, the expression (3.53) can be differentiated with respect to 
C, (or C2) under the integral sign, and the corresponding derivatives 


Of. {H(C)} _ 8. [tanh [E(f)/26] . 
~ OC, = 2geCa— (2x8 E(/) ay BeCpAa f | Aah) af, 


(3.54) 








FeO) « 29.0, fg ( SALI > gepap| aia 


will be continuous functions of C in the whole space of points (C). 
It is not difficult to see that an analogous treatment is valid for 
partial derivatives of f,,{H(C)} of any order with respect to the 
variables C,, ...,C, Cy, ..-5 C3. 
In fact, on further differentiation of the expression for U, in addition 
to the factor {tanh [E(f)/26]}/E(f), the expressions 


| 1 @ (= I 
E OE E E=EU) 


1 @ (1 @ tanh[E/26] 
E OE (x OE E Viensen 


also appear. These are bounded functions of E (since {tanh [E/26]}/E 
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can, for small Z, be expanded in a Taylor series in even powers of E), 
and, as E + o, fall off like 


1 const 1 const 
Be pee! ae op 
Moreover, on differentiation of U with respect to the variables C, 
there appear further polynomials in C,,C* and 4,(f), which also fail 
to invalidate the above arguments. 

Returning to the expressions (3.54) for the first derivatives, we see 
that, since they are continuous functions of C at the points C = C at 
which the absolute minimum of the function f,.{H(C)} occurs, we 
have 








26.0:— fs | sen {y chan} 2(f) df = 0, 
(3.55) 
280. aes | EGY [HCOMD| MN) A= 0. 


Thus, summarizing the results just obtained, we see that the following 
theorem holds: 


THEOREM 3.III. If in the Hamiltonian (3.29) the operators T and J, 
have the form (3.2) and the functions 2,(f) satisfy the conditions (3.37a), 
(3.38) and (3.39), then: 


i. | fv {H(C)}— fo {H(O)} | by (3.56) 
for 
[C.J 2M1, x2=1,...,8 


where 5, +0 uniformly with respect to 0 in the interval (0< 0=<6o). 


Here, f.{H(C)} is given by the expression (3.53) and possesses 
continuous partial derivatives of all orders with respect to the vari- 
ables Ci, ...,C,, Cj, ..., Cy for all complex values of these variables. 

This function has an absolute minimum in the space of all the 
points (C), which is realized at certain points C = C: 


min f..{H(C)} = f.{H(C)}, 
(C) 


satisfying the equations (3.55). 
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2. The inequality 
fH) —f..{H(C)}| < dy (3.57) 
holds, where 5, = (e(1/V)+5y) + 0 uniformly with respect to 
6 in the interval (0 < 6 < 6p). 


We shall now add a rider to this theorem. 


§ 6. RIDER TO THEOREM 3.11] AND CONSTRUCTION OF 
AN AUXILIARY INEQUALITY 


The point C = C at which the function f,,{H(C)} attains an absolute 
minimum is, in general, not unique. However, in the particular case 
when the absolute minimum is realized at the point C = 0, the 
uniqueness property holds. 

In other words, if 


wy S-{H(C)} = fof H(0)}, 
then 


So{H(C)} > f{HO)} 


for C # O(ie. for C such that at least one of the components C, ~ 0). 
To establish this property of the free energy (3.53) taken for the trial 
Hamiltonian, we shall assume that the opposite is true. 

Then there exists a point C # 0 such that 


f-{H(C)} = foo{HO)}.- (3.58) 
We put 
C=C, 1>0, 


and consider the function 
9(0) = fe {H(V/2C)}. 


Then, because of (3.58), 
$(1) = $(0). (3.59) 
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Making use of the expression (3.53) and differentiating, we find: 


d h 26 5 
HO) = 28 BIC. gue ( MELO 9 ereasinf a 








d°p(t) _ ef/6 sinh sinh E)—EI6 E/6 < 
de — (2n)8 sar + eb els é xn/ a 


Since 


sinh (E/9)-E/9 _ 5 
E3 > 


d°/dt? can go to zero only if, for all f, 
| Lacnty| =0 
identically. But in this case, 


| 
(= Peco } ae 
also, so that 


#O) =2SglC> 0, 1H0. 


But this inequality contradicts (3.59). Consequently, 





FA) 0; (3.60) 


On the other hand, since C = 0 gives the absolute minimum of 
F..{H(C)}, we have 


o(t) = 90), t>0. 


Therefore, (d#(t)/dt),9 cannot be negative: 


Cae 
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Hence it follows from (3.60) that 


aR) 4 for t>0O 
dt 





and, consequently, 


$1) > $(0), 


which again contradicts (3.59). Our rider is thus proved. 

To conclude this section, which contains preliminary results relating 
to the properties of the free energies f(A), fy{H(C)} and f.{H(C)}, 
we shall prove one more inequality, which we shall use frequently 
in the following discussions. 

We shall consider systems defined by a Hamiltonian which depends 
linearly on some parameter T: 


H, =Ilottl. 


We shall formally define the expression 


SA) =- a In Tr e~ #26, 


which we shall call the free energy per unit volume for the model 
system H,. Differentiating this expression, we have 


d 1 Trl ye—#9 | 
Ziv) = “Tress = p Pie (3.61) 

and 

1 

af) 1 Tr {Pye~ Pee A-9} ge 
av? ~ 6 Tr e—H:/9 ? 

0 

where 


Pr, = I\—W1)u, - 
But, as we have shown in §1 of this chapter, 


a*fy (H;) <0 


ar (3.62) 
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in view of which, 


dfy(H:) dfy(Hx) _ {dfv (Hz) 
A |< cde (ot Cans) 


and, therefore, for the difference 
1 
d 
ful orl) ful) = [ FAM 


0 
we obtain the inequality 


ee | <fyp(Tot P1)—frP 0) = =e | 


Thus, on the basis of (3.61), we have established the following 
important inequality : 


t™=0 


+ Piren =Sv(Pot Ll) —-fr(P 0) = Z (Ti)r,. (3.63) 


We shail make use of these inequalities later when we specify concrete 
model and trial systems and choose the source terms in an appropriate 
way. 


§ 7. ON THE DIFFICULTIES OF INTRODUCING QUASI- 
AVERAGES 


We shall now study the question of the determination of quasi- 
averages. Let of be some operator of the type for which the limit 
Theorems II, III and IV were formulated in Chapters 1 and 2, e.g. a 
product of Fermi amplitudes, field functions or similar operators. 

Then the quasi-average 

<A> 
of such an operator will be defined, for the Hamiltonian (3.29) under 
consideration, as the limit 


<ot>n = lim ( Jim (ct)r (3.64) 
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of an ordinary average 


(Aé)r 


taken over a Hamiltonian J’ obtained from H by adding “source 
terms” to it: 


P= H-VY (di + visa) 


= T—2V LIS? = vy (VIE + Vea). oO) 
We now wish to call attention to certain difficulties associated with 
the definition (3.64). Thus, in the definition given, there is no indica- 
tion in which region the parameters » must lie, or how they must tend 
to zero in order to ensure convergence in the definition (3.64). 
We shall show that, even in the simplest cases, if |»| tends to zero 
arbitrarily, the limit lim may not exist. 


»2+0 


We shall take, as an example, the Hamiltonian (1.1): 
H =T—2VgJJ*, 


which we have examined in a number of papers;’ the basic results 
of these were summarized briefly at the beginning of Chapter 1. 
We recall that here T and J are given by the formulae (1.2), and the 
function A(f) satisfies all the conditions imposed in Chapter 1, § 1. 
For I’, we took a Hamiltonian with real positive v: 


P =P, =T-WegsJIt—wW(J+I4, (v=0). — (3.66) 
As we have shown, I" reduces to the form (see (1.11)) 


Pr =T,—Wg(J—C(v)) (Jt —C()), 
r,=T-} 4h) {atat,+a_,a,}+2gVC?, 


Af) = 2f) e+ ze}. 


t Bogolyubov, Jr. (1966 a, c, 1967). 


122 A METHOD FOR STUDYING MODEL HAMILTONIANS 


Here, 
CO)+ 55 >0 


and the quantity C = C(v) realizes the absolute minimum of the 
function f,, {7(C)}: 


min f..{(C)} = fo{P(CO))}.- 
(Cc) 


Moreover (see (1.13)), 
(J-C()) Jt¥-C()))r<ev +0 (WVr+e). (3.67) 


It can be seen that the Hamiltonian I belongs to the class (1.14), and, 
by virtue of the inequality (3.67) and the conditions imposed on A(f), 
the conditions I of Chapter 1 and I’ of Chapter 2, § 7, are fulfilled. 

Because of this, we can make use of the above-mentioned limit 
theorems and establish the existence of limits of the type 


lim (c#)r = lim (c#)r,. 
V— co V—> oo 
We note further that, as has already been pointed out (see (1.7)), 
Ci”) > C0)=C (v>0,»-0). (3.68) 
On the other hand, the expression 
lim (A), 
V —+co 


can be expanded in explicit form using the rules of Bloch and de 
Dominicis, and it can be proved in a completely elementary way that 
the passage to the limit » + 0 (y > 0) can be made and reduces 
simply to replacing C(v) by C in this expression, i.e. to replacing the 
averaging over I”, by averaging over H(C). 

In this way we can establish the existence of the quasi-averages 


y+ 0\ V-+oco 
0) 
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With the definition (3.64), and for the case when 
C #0, (3.70) 


we shall now examine how the situation changes when we go over 
to complex values of v and, in place of the Hamiltonian (3.66), take 


Ty,.8 = T—2VgJJ+ —V(vs* + 9*J). (3.71) 


Here, we shall put » = |»| e” and note that I’, ,. can be reduced to the 
form I’ =I\,) (i.e. to the Hamiltonian (3.66) with |»| in place of ») 
by means of the gauge transformation 


ar > arei@2), af > af e—iol2), 


Thus, we obtain, for example, 
(af (t) at At))r, y= e'? (at (t) at,(t))r,,, 


‘ 3.72 
= i (at (0) atp(t))r,, Soe 


The limit 
| iat Hit (af (0) at ,())r,,, (3.73) 
obviously exists and is given by the formula (1.8) in which C replaces 
C in the expressions for u(f), of) and E(f). Then, in the case (3.70) 
under consideration, the expression (3.73) does not go identically to 
zero. 
Consequently, although 


im (af (0) at ,(t))r, ,. 


always exists for |»| > 0, the limit 


lim hin (aft (t) a*,(t))r, :. (3.74) 


v—>0 


does not, for the trivial reason that the ratio »/|»|does not tend to 
any limit as » + 0. 
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The limit (3.74) exists only when » tends to zero in such a way that 
the ratio »/|»| is finite. 

In the general case (3.65), with the passage to the limit » > 0 the 
situation is found, naturally, to be even more complicated. Apart from 
gauge invariance (due to the gauge group), other groups of trans- 
formations can also occur, e.g. the rotation group. 

We shall now direct our attention to a difficulty which is specific 
fors> 1. 

We take the Hamiltonian 


P= H-V{(J,t Jp) + Jot Jb)}. 


Here, we take »; and v2 to be real and positive. 
We put here 


where the operators J and T have the same form as in the Hamil- 
tonians (1.1), and (3.66). 

In the given case, H will, in this way, be the same Hamiltonian 
(1.1) that we have just considered. 

We take 1 = v2. Then the source terms will drop out completely: 


T=H 


and, since the operator H conserves the number of particles, we have, 
identically, 
(afaty)r = 0. 


It can be seen that in such a situation we cannot define a quasi- 
average correctly at all. 


§ 8. A NEW METHOD OF INTRODUCING QUASI-AVERAGES 


In order to avoid difficulties of the above type, we propose that v be 
taken proportional to C with positive proportionality coefficients: 


Yy=nCe, te>0, «2 =1,2,...,5. (3.75) 
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In such a case, we shall consider the trial Hamiltonian 


r, a T-2WY g.{C.Ji+Cr Ja} 


3.76 
VY r(CJt + CtJ.)+ const. ee 


Here we shall not write out the constant term, since it affects neither 
the calculation of the averages (...),, nor the equations of motion. 
It can be seen that we have obtained a trial Hamiltonian for H with 
the transformed parameters 


+ Ta 
8a 8a 2 > 


Vo 


H,=T-WY (s+ ‘ ere 


In order that I, from (3.76) be a trial Hamiltonian not for H, but 
for the original H, we must replace g, by g,—7,/2 in the expression 
(3.65) for I’ (in which »y is taken in accordance with (3.75)), thereby 
obtaining 


Yq 


P=T-wy (s.- f rest VY r(ItCut IoC). 


It is clear that here, apart from adding “sources”, we have performed 
a “renormalization” of the parameters g,. 

We can add any constant term to this expression for I" since it will 
not affect either the average (...), or the equations of motion. 

As such a constant term, we shall take 


VY rCIC.. 


Then the Hamiltonian I will be represented by the form 
T=T-wy (8-H) Joti VY r(IEC,+IC2) 
. 7 iy (3.77) 
+V¥7rC.C2 = H+VY r(Ja—Ca) (JE — C2). 
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We emphasize that here, as always, C denotes the point at which the 
absolute minimum of the function f, {H(C)} (see (3.48)) is attained. 
For notational convenience, in (3.77) we put 


r, = 2%.2, where t, > 0. 
Thus we shall be concerned with a Hamiltonian having the form 
r= H+2V ¥t.8(IJ—Ce) (Jf-CD =T- vy g.J,Jt 

+2V Y) tee(Ja— Cx) JE —Co)s 


T > 0, a =1,2,...,5. 


(3.78) 


We shall show that, with the above choice of I’, no difficulties will 
now arise in the definition of the quasi-averages 


<A>y = lim lim (ea)r, u>O0, a= 1,...,8. 


0 Vo 
For this, we note first of all that, for 
: 5 ee Cr t; =o, 
from (3.78) we shall have: 
P=T-WY gdJiC.+ JeC+2V > g.C.Czr = H(C) 


(cf. also formula (3.84)). Since 


(J.—C,) JE — C2) = 0, 
we see that 
A(C)-l =0, 
and (3.79) 
r'-H=0 for 0<1,< 1. 


Consequently, the inequalities 
flHO)}=frl)=frlH) (for O<%<1) (3.80) 


are valid. However, 


0 = fy{H(C)}— fA) =| fool H(C)} fr (A) + [F-fH(O)}-fAHO)} |. 
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Therefore, on the basis of Theorem 3.III, we obtain 
0 =f, {H(C)}—fr (A) = 5+ 4r. 
Hence, taking (3.80) into account, we find 


0=<f/()—f(H) < 5+ 6r, 


(3.81) 
0<f/{H(C)}-f(T) « 5y+5r. 
We make use now of inequality (3.63) and substitute into it 
Po=H, PF, =l-H = WY t.8(Je—C.) (Jt —CD. 
Then, from the first of the inequalities (3.81), we obtain 
2Y; t8a((Ja— Cx) (JE —CD)r = 5y + by. (3.82) 


We have thus proved the following theorem: 


THEOREM 3.1V. Let the conditions of Theorem 3.111 be fulfilled, and 
let I be represented by the expression (3.78) in which 


O<1,<1, «=1,2,. 
Then the following inequalities hold: 
0O=<f(l)-f(H) <5y+6y >0 as Vo, 


Y 8.((Ja—C,) JE — Ca)rw tov by (3.83) 


>-0 as Vo, 
where to is the smallest of the quantities t1, T2, ..., T,- 


RIDER TO THEOREM 3.1V 
We shall consider the more general case when, in the expression 


H =T-WY g,J.Jé, 


the operators T and J, are not represented by (3.2) and satisfy only 
the conditions of Theorem 3.II. 
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Then, replacing Theorem 3.III by Theorem 3.II in the arguments 
carried through above, we see that Theorem 3.IV remains true. 

We note, further, that in the case when the operators have the speci- 
fic form (3.2) and the conditions of Theorem 3.III are fulfilled, the 
Theorem 3.IV proved above enables us to transform the Hamiltonian 
I directly to the form (1.14), (1.15), and the conditions I (Chapter 1, 
§ 1) and I’ (Chapter 2, § 7) are found to be fulfilled. 

In fact, we have 


H =T-W Yi g(JiC.+J.C2)+2V » BCC. 
—2V >" Ae C.) (JE -Cd, 
and, therefore, ; 
P= H(C)-2v yall —%.) (Jz— Cz) (Je— CD; (3.84) 
where 


H(C) = T-2V ¥. g(JiC.+JeC2)+W Yg.C.C?. 


It can be seen that this Hamiltonian has the form of the Hamiltonian 
(1.14), (1.15) in which we put 


r, = HC), 
Af) = 23 BAA NC2, 
K=WY 8.0. S89) 
G, = 2¢.(1 1.) a0). Ce Cy. 
By virtue of Theorem 3.IV, the following inequality is fulfilled: 
YG. (Ia Cu) Ia ED) = Fy, (3.86) 


where 
by+4 
et DUTY iy (V +) 
To 


uniformly with respect to the temperature 0 in any interval of the form 
(0 <O« Oo). 
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The validity of paragraph 3 of conditions I is thereby also es- 
tablished. 
The remaining paragraphs of conditions I and I’ follow trivially 
from the inequalities (3.37) and (3.38), the condition (3.39), the fact 
that s is finite in the sums over «, and the fact that the quantities C 
are independent of V. 
We can therefore make use of all the limit theorems proved in 
Chapters | and 2. 
Since ., = H,(C), we write the theorems on the existence of the 
limits 
lim (oA#)p = lim (cA)r, 
V-—+> oo V— co 

in the form 
jim (of)r = lim (cbmc). 


But H(C) is independent of the parameters t in the case under consi- 
deration when 


O<1t,<1, w=1,2,...,5. (3.87) 


Therefore, the expression 
lim (cA)r, 
V—> co 


is also independent of the t lying in the region (3.87). 
Consequently, when all the t1, ..., t, tend to zero while remaining 
positive, we have, trivially 


lim lim (c#)r = lim (A)ae. 
t7>0 Vo V— co 


We can define the quasi-average in this situation by the relations 


<A> = jim (A)r = Vis (A) He) (3.88) 


in which the t can take any values from the range (3.87). 

We emphasize again that the most important point in our arguments 
was the establishment of the inequality (3.82), based on the inequality 
(3.81). 
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It can be seen that from the inequality (3.81) there follows the 


asymptotic relation 
jim f(D) = lim fH). 


(3.89) 


§ 9. THE QUESTION OF THE CHOICE OF SIGN FOR THE 


SOURCE-TERMS 


We note, in passing, that the above asymptotic relation ceases, in 
general, to be true for negative values of t. In fact, we take, for 


example, 
Te. =—~Wz3 o,.>0 (« =1,2,...,5). 


Then 
r= Ee = H-2V ¥@,8.(J.—C,) (Jt — CF). 


We make use of the inequality (3.63), substituting in it 


Io = r, Tr; = WV YY w.8AI.—C,) (Jt — C2). 


Then in (3.63) the Hamiltonian will be 


Potlhi = H 
and 


fM)—-fy(Co) = 2 YO.82((Je— C.)(JZ—C))a- 


But, obviously, 
(afar p)n = (a_,4y) nH = 0, 
and therefore 
Jan = 0, Wt) n =0. 


We have, consequently, 
(Ja—C.) JE-CD)n = Sant |Cal*, 
whence, by virtue of (3.91), we find: 
fH) —fr Po) = 2) @eBe IC.I?, 


(3.90) 


(3.91) 
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and, passing to the limit, 
_ SAH)— lim fp(Po) > 2Y ange |Cu|*. (3.92) 
—> 00 V —> co « 


Thus, if C ~ 0, the equality (3.89) is not true for negative ¢ (3.90). 

It was to take this fact into account that we required that the 
proportionality coefficients r, in our choice (3.75) of parameters », 
characterizing the sources included in the Hamiltonian be positive. 


§ 10. THE CONSTRUCTION OF UPPER-BOUND INEQUALITIES 
IN THE CASE WHEN C = 0 


We give here a separate treatment of the case when C = 0, ice. 
when 


CpeCr= =, S0 (3.93) 
H(C) =T, 
SAA)-f(T) +0 as Vo, 


Thus, the interaction terms 
-wWV Y gJJe 


l<a<s 


In this case, 


and, therefore, 


of the Hamiltonian H are asymptotically (as V + -) ineffective in the 
calculation of the free energy. 
Further, we have (compare with (3.78)) 


T= rr, = H+WY g.teJaJt = T-2Y» &.(1—t)JaJIt (3.94) 


and, in view of our earlier proofs, we can write an upper bound for a 
correlation average constructed on the basis of this Hamiltonian: 


/ Yat) sp OG ae VSs - GSH) 
* r 


x 2% 


We shall show that, in the given case (3.93), we have also 


(y fJuJt <t,-0 V+). (3.96) 
\‘e H 
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For this, we shall take the Hamiltonian 


H, =T-2W(1+o)>)g.J.J¢ (1>o> 0), (3.97) 


and formulate a trial Hamiltonian 
H,(C) = T—2V(1+0) ¥ g(C.J¢+C2 J) 


4+2V (140) F gCiC,. B78) 
We denote by C™ the point C giving the absolute minimum of the 
function f..{H,,(C)}. If for any positive w, however small, 


C™ = 0, (3.99) 
the proof of the relation (3.96) is trivial. We need only replace H by 
H,, in the inequality (3.95) and for t, take 


@ 


= Te 


in the Hamiltonian ’,. Then the Hamiltonian I in (3.95) coincides 
with H. It remains, therefore, to consider the case when (3.99) is not 
true for some positive value of w, however small. 

We note that the value C = C™ must satisfy the equations 


FoF} 0 (l<a<s); 


i.e., from (3.54), 


1 tanh[E.(f)/26] 


Cp = al, ( manleaCnr20 fy a +arecpanh raf) of 





where 
2 2 
Ff) = (F,-#) +L +or yeep 
Hence, it follows that 


. 1 Q: 
Ol gaap| MD = (3.100) 
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We now put 
(It+@)C, = X, (3.101) 


and note that 
SoH AC)} = fo (H(®)}- 277 wp 8a Xl. (3.102) 


Thus, for 
XP) = (1+.0)CO?, 


the expression on the right-hand side of (3.102) attains an absolute 
minimum. 
Therefore 


foo HX) —2 PY. Ba |X|? < foo HO)}. (3.103) 
On the other hand, 
foe HO) = min fu (HX)), 


in view of which, 


Soo{H(XO)} > fu. {H(0)}. 


Hence, 
O=<f.{H(X™)}—f..{HO)} < eae @ US| |X|? (3.104) 
We shall show now that 
X©+0 (w +0). (3.105) 


In fact, we assume the opposite to be the case. Then, since, by virtue 
of (3.100), X™ is bounded: 


[XO] < (140) 2 <3, 


we can always choose a sequence of positive w’ -- 0 such that 
X@) + ¥ 
with 
Xz 0. (3.106) 
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Putting w = w’ in (3.104) and passing to the limit, we find 
foo{H(X)} = fo.{HO)}. (3.107) 


But, as we have seen, if the point C = 0 gives the absolute minimum 
of the function f,{H(C)}, no other point realizing the absolute mini- 
mum of this function exists; in view of this, (3.107) and (3.106) are 
inconsistent and we have arrived at a contradiction. 

Thus, the relation (3.105) is proved. Noting that 


fo{HO)} = f(T), 
from (3.102) and (3.104) we obtain 


2575 Toe 8a |X|? <= f..{H,(C)}—f.(T) < 0; 


i.€., 


0 <f.(T)— fa. {Ha(C)} < w&(o), (3.108) 
where 
E(w) = Tia 28 [X2?+0 as w+0. (3.109) 
We now invoke Theorem 3.III. Since, in the case under consideration, 
H(C) = H(0) = T, we can write 
IfAH)—f.(T)| < dy +0, V+). 
For the Hamiltonian H,,, we also have 
If Ha)—foo{Hu(C™)}| < dy) +0, as Vro. 
Here, 5,(w) denotes 5, for H,,. On the other hand, 
0 <fXH)—f Ha) = fr(A)—f.(T) + f(T) 
Soo {Ha( CON} + foo {HalO)} fH); 


and, therefore, 


FOUR-FERMION NEGATIVE INTERACTION 135 
We shall now make use of the inequality (3.63); we substitute in it 
Iy=H., T= Woy egJdJt, 
Mo +li = H. 
Then, from (3.110) we obtain 


by+ by(w) 1 
(Letett) «AG + 740) (3.111) 


However, this inequality is valid for any value of w in the interval 
(0 < w < 1), and its left-hand side is entirely independent of w. Con- 
sequently, the left-hand side of (3.111) will not exceed the lower 
bound of the right-hand side in the given interval: 


(Laks) < Cy > 


by+ b(w) 1 
eee 


It only remains for us to show that ¢, goes to zero as V + 0, 

We shall fix an arbitrarily small number 0. On the basis of (3.109), 
we can, in the interval (0 < w < 1) under consideration, fix a number 
@ such that 


fy = inf | 


(0<w<1) 


E(wo) < a. 
We see then that 
by+5wo) , 0 
C ve =e + 2 s 
But, since wp is fixed, we have 


by+ by(wo) 


Xoo 0 as Vo, 
We can find a value Vo such that 
dy+d/(wo) _ @ 
== Siig. 4 oy for V= Vo . 
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Thus, 
Cyv<o for V=Vo; 


ie., ¢, + 0 as V +o, and the relation (3.96) is proved. 
It is now clear that H is of the type (1.14), (1.15), with 


r,=T, Af)=0, C,=0, Ge = 2g. 


By virtue of (3.96), paragraph 3 of conditions I (Chapter 1, § 1) is 
fulfilled; the remaining paragraphs of conditions I and I’ (Chapter 
2, § 7) are trivial in the given case. 

We shall, therefore, make use of the limit theorems proved in Chap- 
ters | and 2. We see that, for the operators of with which these theo- 
rems are concerned, we can write 


V—+ co V— co 


Applying (3.88), we see that, in the case under investigation (C = 0), 


we have 
V + 00 V => co 


Thus, the quasi-average and the usual average of the operators con- 
sidered above are asymptotically equal to the corresponding averages 
taken over the Hamiltonian 7. Moreover, the interaction terms in H 
turn out to have no effect here. 


Footnote to p. 92 (added in proof). Mathematical questions concerning the 
existence and analyticity properties of partition functions have been considered 
in a paper by the author (Physica 41, 601 (1969)). The theorem proved there 
has been further generalized by H. D. Maison in Preprint Th. 1299-CERN 
(Geneva), 25 Feb. 1971. 


CHAPTER 4 


MODEL SYSTEMS WITH POSITIVE 
AND NEGATIVE INTERACTION 
COMPONENTS 


IN THIS chapter we consider model systems with positive interaction 
components and also systems containing both positive and nega- 
tive four-fermion interaction components. We construct, under the 
conditions I of Chapter 1, § 1, an asymptotically exact solution for 
such model systems, 

The technique developed here leads to a minimax principle for 
model problems with four-fermion interaction, and enables us to find 
asymptotically exact expressions for the free energies and many-time 
correlation averages 

A second, equivalent, formulation of the minimax principle is also 
proposed. 


§ 1. HAMILTONIAN WITH NEGATIVE COUPLING 
CONSTANTS (REPULSIVE INTERACTION) 


We shall consider the Hamiltonian obtained by replacing g, by —g, 
in the Hamiltonian studied in the previous chapter: 


H=T+W Y. g.J.Jt, (4.1) 
1 


<e=xr 


&>O0, a=l1,...,7. 
137 
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Reasoning by analogy, we shall attempt to construct, for this H, the 
usual form for the corresponding trial Hamiltonian 


HC) = T+2VY 8.C.JE+CtI.—C.Cd). (4.2) 


We have 
H = H(C)+2V ¥ 8.(Ju—Ca) (Jt — Cz) = H(C). 


Now, therefore, in contrast to the preceding chapter 
SH) > fr {H(©)}- (4.3) 


This inequality is fulfilled for any complex C. Consequently, if we 
wish to select the constants C in (4.2) in such a way that f,{H(C)} 
be as close as possible to the expression for f,(H), we must determine 
C from the absolute maximum condition 


C=C, 
SV{HC)} = na Sv {H()}. (4.4) 


We shall assume now that the operators T and J, in the Hamiltonian 
(4.1) under consideration satisfy the conditions 


\fT)| «Mo, |J.l<= Mi, 


(4.5) 
Mo, M, = const. 
We introduce the expression 
FC) = fr{H(C)}+2¥) 8.C.Cz 
(4.6) 


= fil T+2V Y sCatt +Cz1d} 
and note that 


OF y(C OF Y(C 
ee = 2¢,(JD) Ho» iss ) = 28.(JaHC)- 
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Therefore, by virtue of the conditions (4.5), 


OF (C) oF y(C) 
ac. 


<= 2g,Mi, act = 





<2g,Mi. (4.7) 














On the other hand, the problem of the absolute maximum of the 
function 


SAO) (4.8) 


is equivalent to the problem of the absolute minimum of the function 
—fv {H(C)} =—Fy(C)+2¥ g.Cic, . 


Therefore, taking the inequalities (4.7) into account and repeating 
the arguments from the beginning of the previous chapter, we con- 
vince ourselves of the existence of quantities C =C realizing the 
absolute maximum of the function (4.8) in the space of all the points 
C (Ci, ..., C,) (with complex values C,). 

Since the function (4.8) possesses continuous partial derivatives 


of all orders with respect to the variables ...C,...C?..., we see 
that the values 
CSC. ClaCh (ea ht) 
are solutions of the equations 
OfV{HC)} _ OFAH(C)} _ 

as 0, ace 0. (4.9) 

Hence, 
Ci = Juda, C2 = Ud)me)- (4.10) 

Consequently, on the basis of (4.5), 

[Cee Mi: <= Nyse, Fe (4.11) 


We note that in the case of the model system (4.1) the point C = C is 
unique. Moreover, we shall show now that the solution of equations 
(4.9) is unique. For this, we construct the expression 


Sf AH(C+ tz)} = Fy (C+ tz)—2¥ (Cr+ tz) (Cut tzu), ver 


Z = (21, ..., Zp) 
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(z, is an arbitrary complex number.) But 
Fy(C+tz) =fvr(lot+ 1), 
where 
Py = T+2V ¥8.(C.JiE+Cz2J.), 
r, = WY glZaIt +22 a)s 
and, therefore, by virtue of (3.62), 
@ 


Hence, it follows from (4.12) that 


Pf {H(C + tz)} 
YD 


7 =—4 » Bal Zn » 


(4.13) 


(4.14) 


It is not difficult to see that the equations (4.9) have the solution which 


gives the function (4.8) its absolute maximum. 


Let the equations (4.9) also have some other solution C = C. We 


put 
z=C-C, C=C. 


Then, with these values of z, 


dfy {A(C+ tz)} 


at =0 for ¢=0 and t=1, 


whence 
1 
afy{H(C+tz)} , 
(eS dt=0 
0 


and, therefore, by virtue of (4.14), 


Yi 8.2072 = 0; 


ie, z= 0. 


Thus, we have proved that the solution of equations (4.9) is unique 
and gives the absolute maximum of the function (4.8). Finally, we 
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shall make use of the inequality (3.34), in which we put: 
Ty = HC), Pr; = 2V YY ga(Ja—Cu) (JF — C2), 


Pot li = A. 
We obtain: 
2 Y8a((Sa— C.) (JE -—Co))n = f(A) —f{H(O)} 
= 2) a((Ja—C,) (JE -—CD) ney. (4-15) 


After this preliminary discussion of the properties of the form of 
the trial Hamiltonian, we return to our consideration of the Hamil- 
tonian (4.1) itself. 


§ 2. FEATURES OF THE ASYMPTOTIC RELATIONS FOR THE 
FREE ENERGIES IN THE CASE OF SYSTEMS WITH POSITIVE 
INTERACTION 


Here the question arises as to whether a theorem analogous to 
Theorem 3.I is true for a Hamiltonian of the type (4.1). 
In other words, can we assert that 


fiH)—frlH(C)}+0 as V +o (4.16) 
when the operators T and J, occurring in (4.1) satisfy the conditions 
(3.30) and (3.31)? We shall show now that such an assertion is not 
true even in the simplest case r = 1, when 











H=T+2VgJJ* . (4.17) 
For this, we shall consider an operator J of the form 
1 
= + 
J WV py A(fapat,, (4.18) 
where Pp 
Ad(o), Om | = A, 
Af) = (p, 0) = : 
0, f- el > A, 
ie (4.19) 
140) = | DOs. 
—4, o=-+; A=const > 0, 


A = const > 0. 
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For the operator T in (4.17), we take 
T= Ayo = To—2VgJJ* * 


T= 5 (Sn) 7 


This Hamiltonian Ho obviously belongs, by virtue of (4.18), (4.19) 
and (4.20), to the class of Hamiltonians considered in the previous 
chapter, and satisfies the conditions of Theorem 3.III. 

Here, we choose the values of g, J and 6 in such a way that the 
quantity S = 5 which realizes the absolute minimum 


om So{Ho(S)} = fue {HolS)} 


(4.20) 


is non-zero. Then 
foo{Ho(S)} < f.{Ho(0)} = f.-(To), 
and, on the basis of Theorem 3.III, 


jim fv (Ho) = fo {Ho(5)}. 


Consequently, we have 

jim {fv (To) —fV(T)} > 9. (4.21) 
On the other hand, by virtue of (4.18), (4.19) and (4.20), the opera- 
tors T and J in (4.17) satisfy the conditions (3.30) and (3.31). 


Therefore, if the assertion (4.16) were true for the model system 
(4.17), we would be able to write 


FAT)-f{HC)}} +0 V+), (4.22) 


since, in the present case, H = Tp identically. 
We shall now find C. We have 


H(C) = T+2Vg(CI++C*I-C'C). 


As has been pointed out, the value C = C is determined from the 
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equations 
Ofv{H(C)} _ PVAHCO)} _ o. 
aC = 0, aes aaa 0; (4.23) 
i.e., 
C= Due, C= Io: 


It is clear that these equations have the trivial solution C = C* = 0, 
since, by virtue of (4.18), the identities 


(Dao = Jr = 0, 
(J*) Ho = (J*)r = 0 


are true. Since the equations (4.23) have a unique solution, we have 
C=0, H(C)=T. 
Therefore, from (4.22) we obtain 


lim ff (To) —f(T)} = 0, 


which contradicts (4.21). Thus, in contrast to the case of negative 
interaction (3.29), in the positive-interaction case (4.1) under considera- 
tion the asymptotic relation (4.16) can be false even when the con- 
ditions (3.30) and (3.31) are fulfilled. 

The absence of an analogue of Theorem 3.1, which served as the 
basis of our investigations in the preceding chapter, can, generally 
speaking, make the study of the model Hamiltonian (4.1) in the 
general case more difficult. 

However, when the operators T and J, have the special form (3.2), 
we obtain an even simpler situation than in the preceding chapter. 


§ 3. BOUNDS FOR THE FREE ENERGIES AND 
CORRELATION FUNCTIONS 


Let the operators T and J, occurring in (4.1) be represented by the 
formulae (3.2). We shall also assume that the functions 


AAS) = ACP, 7) 
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satisfy the following conditions: 


(2) FL IAlp. )| <Qs = cons 


(b) |A.(p,o)| = Q = const, (4.24) 
(c) The discontinuities of the functions 4,(p, 0) form a 
set of measure zero in the space E. 


We note, in passing, that, in particular, it follows from these condi- 
tions that 


7 |A.(p, 6) |? < QO, = const (4.25) 


where as Q» we could take Oz = OQi. 
We shall consider the equations (4.9) for C and write them in the 
form (4.10): 
Ci = Janey, Cr = Jt)ue- 


Jar = Jt)r = 0 


identically and H(O) = T, we see again that the equations under 
consideration have the solution 


C.=0, Cr=0, 


Since 


and, by virtue of our earlier proof, this solution is unique. 
Thus, 
H(C) = H(0) = T. 


We now make use of the inequalities (4.15) and obtain 
0 <=fy(H) f(T) <2 Y 8a IeIt 7, (4.26) 
2 » 8a Jade) nx =f (H) —fy(T). (4.27) 


We shall write out the right-hand side of inequality (4.26) explicitly. 
We have 


Uedtr = aya EMAL) afatyap apr 
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But the average 

(ap at pa_pdp)r 
can be non-zero only if f’ = f or f’ = —/f. Further, 


(ajar pa_,ay)r <1 
and 


| afar para_z)r| = \(apaatpa_r)r| <1. 
Therefore, 


1 1 
+ Pike y Bg 
Jade r< we | ACS) | <= V Q2, 
whence, from (4.26) and (4.27), we find 


0< fH) fT) <7 SB, (4.28) 


(Lalit) <7 abe (4.29) 
« H « 


Thus, the interaction terms 
V aeIeJt 
a 


contained in H have, asymptotically (V — ©), no effect on the free 
energy: 


fAA)-fUT) +0 as Vo. 


Furthermore, it is not difficult to see that H reduces to the form (1.14), 
(1.15), if in the latter we put 


Pr=H, f,=T, C,=0, A=0, G,=—2g,. 


By virtue of (4.29), paragraph 3 of conditions I of Chapter 1, § 1 
is fulfilled. The remaining paragraphs of conditions I and I’ can be 
verified trivially starting from (4.24). 

We can therefore use the limit theorems of Chapters | and 2. 

It follows from these that, in the given situation, averages 


(A)n 
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of the operators cf to which these theorems apply are asymptotically 
equal to the corresponding averages 


(A)r 


taken over the Hamiltonian T, so that the interaction terms in H are 
also asymptotically ineffective here. 

In the case being considered, the quasi-averages do not differ from 
the asymptotic values of the usual averages, inasmuch as they are 
defined in terms of the Hamiltonian 


P,=T+WY gl +t)JaJt , 


which has the same form as H, except with the transformed coupling 
constants: g, > g,(1+t,); in view of this, averages of the operators 
cA taken over the Hamiltonian I’, are asymptotically equal to the cor- 
sponding averages over T, provided that t, > —1. 


§ 4. EXAMINATION OF AN AUXILIARY PROBLEM 


We shall consider a model Hamiltonian of rather more general 
form, namely: 


H=T+2 > {Ulf +Uglj}+2 > gi JJt (g. > 0), 
1<s<! l<a<r 


(4.30) 
where 


T= » (Sa — nu) aay, 
J, = wh jafaty, p= a > ud fata, (4.31) 


AA—f) =—A/(f), Ma —f) =— pS), 


and it is assumed that the functions 2,(f) and y,(f) satisfy the condi- 
tions (4.24). 

Furthermore, here U, arecomplex parameters which can take any 
finite values. 
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It is clear that a Hamiltonian of the type (4.30) reduces to the 
form (4.1) if in the latter we put 


T =T t+ ue {Upp + Ugle}. (4.32) 
It is easily verified that the conditions (4.5) are fulfilled, with, in this 


case, 
M, = Qi. 


Therefore, the form of the trial Hamiltonian will be 
H(C) = Tot 2V y {UI + Uplp} 


+W Y g.{C.Jt+CtJ.—C.C%). (4.33) 
The quantity C, is determined from the equations 
P{HC)} _ Ofv{H(C)} _ 

aCe th = 6G On (4.34) 


According to our earlier proofs, these equations have a single unique 
solution. This solution gives the absolute maximum of the function 
JS, {H(C)} in the space of all the points C. 

Since the expression (4.32) contains terms J, and / +, averages (J,)r 
taken over the Hamiltonian T (4.32) are no longer necessarily trivially 
zero; thus, the solution of the equations (4.34) can be non-zero. Since 
this solution depends in general on V, we shall denote it by 

C=C, 
In accordance with (4.11), 
|C|<0,, «=1,2,...,7. (4.35) 
For brevity, we put 
H{C™} = Ha, (4.36) 
and write, using (4.15), 


0 <f(H)—fr(Ha) <= 2¥ 8a((Ja— CX”) (JE -— O°) wv, . (4.37) 


We shall now write out the right-hand side of this inequality explicitly. 
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It follows from (4.34) that 
OM = (Sdn, CP = Ta, 
and therefore we have 
(Ja 7 oO) (Jz 3 CO a, = Ja Je\an, 7 (Ja) Ht, (Jn: 
Consequently, 
(J.-C) (JE - CLO uz, 
1 ; 
= 42 p> A) MS) {af at pa-pap)a,— (a a* ya (a-rar)u}. 
(4.38) 


Since the Hamiltonian H, (see (4.36)) is a quadratic form in the Fermi 
operators, we can use the technique of Bloch and de Dominicis to 
expand the expression (aj a*,a_,d,,) 4. Consequently, we obtain 


(afar pa_pay)n,— (Gf a*y)n, Arp) i, 
= (af Op), A= (4_¢ un, X4F Op), SOAs) 
But the first term in the right-hand side of (4.39) can be non-zero 
only if f’ = f, and the second only if f’ =—f. 
In both these cases, the right-hand side does not exceed unity in 


absolute magnitude. 
Thus, from (4.38), we find 


os 1% 1 2 Q2 
(CNIS CON) = 5579 FDP 


Hence, on the basis of (4.36) and (4.37), we see that 


(4.39) 


0<f(H)—fiLH(O)} <2 OY ge. (4.40) 
On the other hand, taking (4.31) and (4.33) into account, we have 
HCO) = To— + Y (A fyapat + Af )a_saj}—2V 1 BCC. 


A(f) =—2 Usd f)-25 BeC2), Cm) 
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and, therefore, 
fr{HO} =-2 Pacey DA: CU), (4.42) 


where 
F(f;C, U) = 3{-E(f)—-T(f)— 26 In (1 + e-2)}, 


FIN=TANHANR TN=F-n. (443) 


By repeating our arguments from Chapter 3, § 2, it is not difficult to 
show that the difference between the sum and the corresponding 
integral is bounded by the inequality 


1 a on 1 7 | —_ 
EPG CO -aae [PLC Od <M +0 4 


as V +o 


for 
[Cel = Qi, |Up| =. 


Here 0,(/2) goes to zero as V + co (M is fixed) uniformly with re- 
spect to 6 in the interval (0 < 6 < 6p). 
Thus, 
If {H(O)}—fo-{H(O)}| < er M0, ne 
\Cal=Q1, |Up| =A, 
where 


fu H(C)} =-2Y, 8CIC.4 Gey | FUf;C, U) df. (4.46) 


From the above-mentioned arguments it follows also that f,,{H(C)} 
is a continuous function of ...C,...Cl... and ...Ug...Ug... 
and possesses continuous partial derivatives, of all orders, with respect 
to these variables. 

In addition, we can verify that a point C = C exists which realizes 


the absolute maximum 


f-{H(C)} = ss fo.{H(C)}, 
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and, for this point, 


IC. | <s Q1, 
and the equations 
Of. {H(C)} _ Of..{H(C)} _ 
= aoe 0, ac? = 0, (4.47) 
c=C 


are satisfied. We note now that, by virtue of (4.45), 


SAH(C)} = f{H(C)} = f..{H(C)}— ev(M), 
fAH(C)} = fa. {H(C)} = fo (H(C}— ov MO), 


in view of which, 
fv {H(C)} — f.{H(C)}| = er(M). 
Hence, taking (4.40) into account, we obtain 
FAH) -fAHC)}| <2 O25 set eM) +0 as V+e, (4.48) 
for |U,| <A. 


§ 5. SOLUTION OF THE QUESTION OF UNIQUENESS 


We shall make a few more remarks on the equations (4.47). 
We shall take the arbitrary complex numbers z), ..., Z, and the 
teal number ¢. We consider the second derivative 


F fol(C+12)) = AY, gattta 
: : a (4.49) 
+ aap | (ge USC uv) df. 


We note that the function F(f; C+ tz, U) is given by the right-hand 
side of formula (4.43), in which, in the expression for E(f), A(/) is 


replaced by 
A(f)+tAi(f) 
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where 


Alf) =-—2) g,234.(f). 
We have, therefore, 


PF(f;C+tz,U) _ I | E 


=e E/o ~~ 
de = — BEI pebDp OG 


x | Af) (Af) + taf) + (AP) 4 Af) AIP)? 
+[47%(f) | Af)? — | Af) (AP) + tai) 


(Af) + tAi(f)) ANY] (e#?— 1) (e#? + of <0, 
whence, on the basis of (4.49), 


a <—4Y g,z%z, <— 4g) ziz,, (4.50) 


where g is the smallest of the numbers gi, ..., g,. Hence it follows 
(asin § 1 of this chapter) that the solution of equations (4.47) is unique. 
These equations have no solution, other than C = C = C(U), giving 
an absolute maximum of the function f,.{H(C)}. 

As has already been remarked, C(U) is defined for any complex 
U =(... U, ...) and is bounded uniformly by the inequalities 


|C.(U)| = Qi. (4.51) 


We shall show now that the C,(U) are continuous functions and 
possess continuous partial derivatives with respect to the variables 
blag actd Ug titan 
For this, we shall make use of a well-known theorem about implicit 
functions. Since, in its usual formulation, this theorem refers to equa- 
tions with real variables, we now go over from our complex quantities 
C, and U, to real quantities. 
We put 
Cy = XatiXages Ug = upg+iupss, 
Cl = Xa—iXa4r, Us = up—inpsi, 
a=],...,7, BH; exh 
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Thus, the expression 
Soo{H(C)} = O(x1, «.-, Xer3 U1, ».-, Uat) (4.52) 
will be a real function of the real variables 
Xi, e+) Kar} U1, ey Uae 


This function will be continuous and will possess continuous partial 
derivatives of all orders for all real values of the variables x and wu. 
Putting 


: : : 
Za =YatiVatr, 2a =\a-Warr, 4 =1,2,...,7, 


in (4.50), we write this inequality in the form 
a . 3 
We d(x14+ ty1, weey Xap tars Ui, .. +5 Uy) =—4g Y Ya + (4.53) 
1l<a=2r 


Noting that 





* — 
Xe = nee > Xatr ifs Ce) = I, ot; 
we have 
ro) 1 @ i 0a 0 1 a rs) 


OC, ~ 2 Om 2 Ou’? Ol 2 Ox) 2 Omar 
Consequently, the equations (4.47) being considered are equivalent 


to the equations 


Bol ay ons ach ty sos Ma) 9, (Q@=1,...,2n. (4.54) 


Therefore, the equations (4.54) have, for any real 11, ... u, the unique 
solution 
Xa = Xa(U1, ..-, Ue), a=1,...,2r. 


By virtue of (4.51), we have 


| Xe(u1, ---, Uet)| = Q1 
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everywhere. Furthermore, taking the inequality (4.50) at t = 0, we 
obtain 
&O(x, u) 2 
2 Fx xr Vee SAB (4.55) 


for arbitrary real ...x,...Ug...y, ..-. But it follows from the fact 
that the quadratic form in y in the left-hand side of inequality (4.55) is 


negative-definite that 
BO(x, Y) 
Det \ eee Ox, Ox! r} <0 (4.56) 
everywhere. 

The conditions for the application of the theorem on implicit func- 
tions to the equations (4.54) are thus fulfilled. Thus, we convince 
ourselves that, everywhere (i.e. for all real ug), the x,(u1, ..., My) 
are continuous functions and possess continuous partial derivatives 
with respect to ... ug... 

Returning to our comple variables C and U, we see that 


C.=C(U) (@=1,...47 


are indeed continuous functions and possess continuous partial deriva- 
tives with respectto...U,... U; ... in the whole space of the points 
(U). 

To conclude, we summarize the results obtained above in the 
following theorem. 


THEOREM 4.1. Let the Hamiltonian H be defined by the formulae (4.30) 
and (4.31) in which the functions 1,(f) and u,(f) satisfy the conditions 
(4.24). Then, 


I | fAH)—falE)| < 2 O25 82+ AD, 


where MM is any number greater than or equal to all the |U,| and 
where, for fixed MM, 


o(M)+0 as V+ 
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uniformly with respect to 6 in the interval (0 < 0 = 6p). 
2. fu(H) = max foo H(C)} 


where H(C) is given by formula (4.33). 

F.AH(C)} is defined by the equality (4.46) for all C and U, 
and is a continuous function which possesses, everywhere, 
continuous partial derivatives of all orders with respect to the 
variables ...C,...Ch.-. Ug...Ug ..-. 

3. The equations 


Pf (HC)} _ 9 _FetHO} _ 
oc. 


2 oC 
have, for complex U,, a unique solution 
Cc =C=C(U) 
This solution realizes the absolute maximum referred to in para- 
graph 2 of this theorem: 
fo(H) = max fue (H(C)} = foo {H(CU))}. 


The functions CU) are everywhere continuous and possess 
continuous partial derivatives with respect to the variables 
eS eee Sean 


§ 6. HAMILTONIANS WITH COUPLING CONSTANTS OF 
DIFFERENT SIGNS. THE MINIMAX PRINCIPLE 


We pass now to the study of a model Hamiltonian with positive 
and negative interaction terms: 


H =T)+2V > SnJaJq—2V Y BSI es 
) 


(lsea=xr r+lea=cris 


Je= ap EM hain MP =-WP. be 


e=1,...,rt+s, 


T= (se-u)ay, 


> 0, 
(4.57) 
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We shall assume that the functions A(f) = A(p, o) satisfy the condi- 
tionst 


(1) |A(p, 0)| =O = const. 
(2) ple | Aa(p, 6)| <O = const. 


(3) The functions A,(p,o) are continuous everywhere, with the 
possible exception of a set of points p of measure zero. 


We note, first of all, that it follows from conditions | and 2 that: 


7 LIAB, 0)| < Qs = const. 


1 (4.58) 
YY Y |A.(p, o) ? < Qs = const. 
V Dp 
We now put 
T=T + Y gdadi. (4.59) 
a=] 
Then, 
rts 
H =T-2V >’ g.J.Jd . (4.60) 
«=1 


Starting from the conditions 1 and 2, it is not difficult to verify that, 
for the Hamiltonian (4.60) with the operator T given by (4.59), the 
conditions of Theorem 3.] are satisfied. 

To verify that the conditions of Theorem 3.II are also fulfilled, we 
need only to prove that, for any fixed complex S, (a = 1, ..., r+s), 


SAHr(S)}—foe{Hr(S)} +9 (V+), (4.61) 


where 


r r+s 
H-(S) = To + 4V Zi &.5,jJt—2V py 8a {S.J + S2I.—S,S2}. (4.62) 


t The conditions (1) and (2) will obviously be fulfilled if, for example, 
|A.(p, 9) | <= 


where A and B are positive constants. 


eae hs 
(p*+ By’ 
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To demonstrate the validity of (4.61), we note that the expression 
r+s 
Hr(S)—const = H7(S)—2V )° g.S.Sz 
a=) 


coincides with the Hamiltonian (4.30) if, in the latter, g, is replaced 
by 2g, and U, by g,S,, and that the conditions (4.24) for the applica- 
bility of Theorem 4.1 are fulfilled. 

Hence it follows, in particular, that 


fol Ha(S))—foo(Hr(S)}| <7 On Y gat oN, 


IM => 8, |Sa\; a#=1,2,....r+s, 


(4.63) 


where 
o(M)-0 V+-) 

for any fixed value of /M. The asymptotic relation (4.61) is thus proved, 
and the conditions of applicability of Theorem 3.1] are found to be 
fulfilled. 

By making use of this theorem, we can convince ourselves of the 
existence of points S = § at which the function f..{H,(S)} attains 
its smallest value: 


fu{Hr(S)} = min f..{H7(S)}. (4.64) 


Further, on the basis of (4.63) and paragraph 2 of Theorem 3.II, we 


see that 
\fu(H)~f{Hx(S))| = By = e( 7) +8, 


in (4.65) 
by = 7 Q2 >» 8. + 0(2GQ1), 


where G is the largest of the constants g1, ..., g, and e(1/V) is, for 
the Hamiltonian (4.60), the quantity corresponding to e in Theorem 
3.1L 
We emphasize that 
by > 0, op +0 as V +0 (4.66) 


uniformly with respect to 6 in the interval (0 < 6 = 6p). 
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We now turn our attention to the rider to Theorem 3.IV. Since the 
conditions of Theorem 3.11 are fulfilled, we can assert, on the basis 
of this rider, that 


0<=f(l)—fA) <= 6y+5,, 
and 


(4.67) 


rts é é 
2 Y gee S)UE-SI)yr, = EE, 


where 


rts 
P= H+2V ‘ave J,—5, Jt —S? > 
X, 8e— Si) ( , (4.68) 


0<7 < I, =), 245 FES: 


% is the smallest of the numbers 7,. 

Having made these comments, we return to the application of 
Theorem 4.1 in the case of the Hamiltonian H;(S). 

Taking paragraphs 2 and 3 of this theorem into account, we see 
that 


fo HAS} = max fo{HC, 9} = fo{Hr(CtS), )}, (469) 
where 
Hy(C, 8) = T + 4V ¥. gx{CaJt+CtJ,—C.C3} 
ae (4.70) 
=2V > 8a {SIT + StJ.—S.Sz}, 
a=1 


and where C = C(S) is determined by the equations 
e {H7(C, S)} = 0 mn H-(C,S)}=0, (4.71) 
act oo AFT» )} =e act f ae ’ Ms at | 


which, for a given S, have a single unique solution. 
It follows from (4.69) that 


foo{Hr(5)} = foo{Hr(C(S), 5)} 
= min fu.{Hr(C(S), S)} = min max f..{H7(C, 5)}, (4.72) 
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and therefore, by virtue of (4.65), we obtain 
lim f(A) = min max f..{H7(C, S)}. (4.73) 
(V2) (Ss) (©) 


We have found here an expression for the asymptotic calculation of 
the free energy by means of the “minimax principle” developed in 
papers by the author (Bogolyubov, Jr., 1968 f, g, 1969). 

We note further that, on the basis of Theorem 4.I, the expression 


foo {Hr(C, S)} (4.74) 


is a continuous function with continuous partial derivatives with 
respect to the variables ...C, ... Cy... Sg...S3 .-.,asisC,(S) with 
respect to the variables ....S,...S3 ---- 
Consequently, f..{Hr(C(S), Ss) is Continuous and possesses conti- 
nuous partial derivatives with respect to the variables ....S,...53.... 
Therefore, if the point S = S realizes the absolute minimum of 


this function, we shall have 


fn {Hx(CS). S)}_ Ae {Hx(CCS), S)} 


a5; : as =0 (475) 


for S = §. But, in view of the fact that C = C(S) satisfies the equa- 
tions (4.71), we can write 


Of.{Hr(C(S), S)} eee =) _8CS) 
C=C(S) 





OS, 1 aC, OS, 
Fa {= {Hr(C, S)} 
OSs c=C(S) 
Thus, 
C, = C, = C,(S), Sp = Se (4.76) 


is a solution of equations (4.71) and of 


Of. {Hx(C, S)} _ Of..{Hr(C, S)} _ 
ar ana er pom (4.77) 


But, by definition of the absolute minimum, 


fo{Hr(C, 5)} < f.{Hr(C(S), S)}. 
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Now let 

Cac. -S=S (4.78) 
be some other solution of the system of equations (4.71), (4.77). Then 
it follows immediately from (4.71) that 


C= C(S). 


fa{Hr(C, S)} <f.{Hr(C, S)}. 
As can be seen, (4.76) is a solution of equations (4.71), (4.77) that 
gives the function (4.74) its smallest value on the set of the points 
(C, S) satisfying the equations (4.71) and (4.77). 
Conversely, if (4.78) satisfies (4.71) and (4.77) and gives the function 
(4.74) its smallest value on this set, then 


foo{Hr(C(S), S)} = f..{Hr(C, $)} = f..{Hr(C, 5)} 
= min fo {Hr(C(S), S)}. 


Therefore 


Thus, the minimax principle for determining C and S, i.e. the 
determination of C = C(S) from the condition that the function (4.74) 
be an absolute maximum in the space of all C for given S, followed 
by the determination of S = § from the condition that f., {H,(C(S), 
S)} be an absolute minimum in the space of all S, is equivalent to 
finding a solution of the equations (4.71) and (4.77) that gives the 
function (4.74) a smaller value than do any other solutions of these 
equations. 

We shall now write the equations (4.71) and (4.77) explicitly. 

Starting from the expression (4.70) for H;(C,S), we find, in the 
usual way, 


folHr(C, S)} =—4 D3 g.CiC,+2 . g.S*S. 
_t ( ; (4.79) 
+ esd {T(f)-E(/)- 26 In a + e—E(fW/0)\ df, 


TN =F B= TNHANE, 
A(f) = 2¥! eed fSI-4 Y. BAI NCE 
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and, therefore, the equations under consideration can be represented 


in the form 
* E(f) 
C.= ms aay | ANAL) mn (FP) 


C= A*(f) A(f) tanh (4p. ) a 


mF | 
‘ Ef) 
Ges rer 5 [ 4002109 am (GP ) a, 


; EY) 
= apap | AU) AD tanh (p>) a 


Hence it follows that, for any solution of these equations, 


C,=S. (@=1,....9, 


and 

foo {Hr(C, S)} = f.{H(C, S)}, 
where 

H(C, S) = H(Cy, «.., Crs Spay «+65 Spats) 
= T+2V os 8x {C.Ji+C2J,—C.C3} 
r+s 
—2V 8 {S.It + S2Ja—SiSz}, 
a=r+)] 

and 


r r+s 
Soo{H(C, S)} =-2 Y BCCe+2 y 8S. bs 
a=l a=r+1 


(4.80) 


(4.81) 


(4.82) 


(4.83) 


(4.84) 


+ aRF | {T(f)—E(f)— 26 In (1+ e- 2") df, 


1) =F - BA) =VPAHADE, 
Af) = 2S) seh fSI-2 ¥ BANE 
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It is clear also that the equations (4.80) are equivalent to the equa- 
tions: 


A@FoLH(C, S)} _ 4 OfeAH(C, St _ 9, 


aC, > ace for «= 1, 2, ex agl’s 
(4.85) 

GoAMC SN) 2g. : SoettC 9. poe gis r+1 s 

A £) ost f) a ete Dy 


with the additional equalities (4.81). On the basis of (4.82), we can 
now state that the method formulated earlier for the determination 
of the quantities C,, ..., C, and Si, ..., S,,, is equivalent to the 
following: 

(a) The quantities 


Ci = Ci, uC, = Cy Sr41 = Ortls ++ +5 Sr+s = Orts 


are also solutions of the system of equations (4.85) that give the func- 
tion 
foo H(C, S)} (4.86) 


its smallest value on the set of all the solutions 5, = C, (a = 1, ...,r) 
of these equations. 

(b) Repeating the arguments given above for f.,{H7(C, S)}, we 
see that such a formulation of the determination of 


Ci, sae Crs S415 2p 419, Sr+s 
is equivalent to the following minimax formulation. 
First, we find 
C = C(S) = C(S41, ---, Sr+s) 
from the condition that the function (4.86) be an absolute maximum 
in the space of all the points C (C-, ..., C,) for given S, and then 


find 
S+41 ste S.% 


from the condition that the function f,, {H(C(S), S)} be an absolute 
minimum in the space of all the points S (S,,,, ..., S,,,). It can be 
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seen that H(C, S) is the usual form of the trial Hamiltonian for the 
Hamiltonian H (4.57). Unlike the form H,(C, S), which depends on 
2r+s complex parameters, H(C, S) depends on only r+s complex 
parameters. Returning to the relation (4.73), we can write it also in the 


form 
on SAX) = manne max f..{H(C, S)}. (4.87) 


V — oo) 


To conclude, we turn to the problem of constructing “Hamiltonians 
with sources” for the determination of quasi-averages. 
As in Chapter 3, we shall take 


lr, = H+ > g,1.(J,—C,) Jt -C*) 
(4.88) 


s 


+2V > g.t.(Ja—S.) (JE —S2), 
a=r+1 
where t, are fixed parameters for which 


O<1<1, aw=l1,...,r4+s. 


It is clear, in view of (4.81), that I”, (4.88) coincides with the Hamilto- 
nian (4.68) and, therefore, the inequalities (4.67) are valid for it. 
Furthermore, we have 


lr’, = HC, S)+W ¥ g(1+t,) Je-C,) JE -C2) 
es sis (4.89) 
—2V > - g.(1 —T.) (J.—5.) (J —S2). 
a=r+ 


Thus, I", reduces to the form (1.14), (1.15) with 


Vr, = HC, 5), 
C.=C., G,=—2¢,fl+t,), «=1,...,7, 
C, = §,, G, = 2¢,(1—t,), a=rt+l,...,rts. 


It follows from the inequality (4.67) that 


r+s 
(¥ iG. eect) ey, 
2x rr 
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where 
y= > vt by +0 (V+) 
To 


uniformly with respect to @ in the interval (0 < @ « 6). Clearly, 
paragraph 3 of the conditions I is fulfilled. The validity of the 
remaining paragraphs of conditions I and I’ follows trivially from the 
conditions 1, 2 and 3 of this Section. 

We can thus apply the limit theorems of Chapters 1 and 2 to the 
above case and obtain for the quasi-averages the same results as were 
obtained in Chapter 3. 
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